横观各向同性页岩井周应力分布模型及代码分享

本文摘要：（由ai生成）

本文介绍了用MATLAB编写的求解横贯各向同性页岩地层井壁坍塌压力的程序。程序接收多种输入数据，如地应力、孔隙压力、弹性模量、泊松比等，同时考虑地应力方位、井眼轨迹和层理产状。程序通过坐标转换处理数据，并基于参数计算井壁坍塌压力。完整代码可通过赞赏获取，本文仅展示部分。该程序为页岩地层井壁稳定性分析提供计算工具，有助于降低钻井风险，提高工程安全性。

如需横贯各向同性页岩地层井壁坍塌压力求解模型MATLAB程序代码，可赞赏后获取。部分代码如下所示，

%% %输入数据

clc;close all;clear;

SH=81.11;          %最大水平主应力，MPa

Sh=58.66;          %最小水平主应力，MPa

Sv=74.09;          %垂向水平主应力，MPa

Pp=39.47;          %孔隙压力，MPa

%%%各向异性程度 k=Ev/Eh

Ev=25.47;        %层理面内弹性模量，GPa(1GPa=1000MPa)

Eh=2*Ev;         %层理面法向上弹性模量，GPa(1GPa=1000MPa)（此程程序中将其设为变量）

vv=0.193;        %层理面法向上泊松比，小数

vh=0.5*vv;         %层理面内泊松比，小数

%% 地应力方位-井眼轨迹

Alpha_s=85;  %N85°E,最大主地应力方向

Belta_s=0;

Belta_b=0;%井眼倾角

%% 层理的产状

Alpha_bp=85;%层理倾向-在大地坐标系的角度

Belta_bp=23.5;%层理倾角-在大地坐标系的角度

%% 本体强度

S_o=21.4;%MPa

u_o=tand(28.3);%角度

%% 坐标转换

Sigma_p=[SH 0 0;0 Sh 0;0 0 Sv];%主地应力

%主地应力向大地坐标转换

R1=[cosd(Alpha_s)*cosd(Belta_s) sind(Alpha_s)*cosd(Belta_s) sind(Belta_s);

    -sind(Alpha_s)  cosd(Alpha_s) 0;

   -cosd(Alpha_s)*sind(Belta_s) -sind(Alpha_s)*sind(Belta_s) cosd(Belta_s)];

Sigma_e=R1'*Sigma_p*R1;%大地坐标下的地应力