1.一种确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,包括:获取水平井多段簇主裂缝干扰应力场,根据岩石力学动态参数及静态参数,确定原始地应力场,根据页岩气藏水平井生产规律,确定孔隙压力降干扰应力场;
将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加,获取现地应力场分布。
2.根据权利要求1所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述获取水平井多段簇主裂缝干扰应力场,包括:采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值,结合二维裂缝干扰应力场计算方法,获取水平井多段簇主裂缝干扰应力场。
3.根据权利要求2所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述二维裂缝干扰应力场计算方法,包括:σy诱导=μs(σx诱导+σz诱导)
其中,σx诱导、σy诱导及σz诱导分别为x,y,z方向上的干扰应力;pn为裂缝壁面上净压力;r为裂缝中心距离二维空间中干扰应力场研究目标点的距离;r1为裂缝底部距离二维空间中干扰应力场研究目标点的距离;r2为裂缝顶部距离二维空间中干扰应力场研究目标点的距离;θ为二维空间中干扰应力场研究目标点偏离裂缝中心的角度;θ1为二维空间中干扰应力场研究目标点偏离裂缝底部的角度;θ2为二维空间中干扰应力场研究目标点偏离裂缝顶部的角度;c为缝高的一半;μs为静态泊松比。
4.根据权利要求2所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述采用每段裂缝的缝内净压力的加权平均值作为缝内净压力值,包括:Li=i*Lf/n
其中,pn0为最大缝宽处裂缝壁面上净压力;Li为距离井筒的裂缝长度;γ为裂缝内净压力分布指数;G为岩石剪切模量;w为缝长对应的缝宽;Lf为等效缝长;n为缝长分割段数;μs为静态泊松比;pni为第i段裂缝的缝内净压力;pn为裂缝壁面上净压力;Es为静态杨氏模量。
5.根据权利要求1所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述根据岩石力学动态参数及静态参数,确定原始地应力场,包括:其中,σz为深度H处的垂向应力;ρr(h)为随深度变化的上覆岩体密度;H为压裂层位深度;g为重力加速度;σH和σh分别为最大、最小水平主应力;μs为岩石静态泊松比;KH和Kh分别为最大、最小水平地应力方向的构造应力系数;α为Biot系数;pp为孔隙压力。
6.根据权利要求5所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述根据页岩气藏水平井生产规律,确定孔隙压力降干扰应力场,包括:其中,Pe为原始地层压力;Pwf为井底流压;A和B为二项式产能方程系数;q为产量;qn为归一化拟产量;Gp为累积产气量;t为物质平衡时间;Zi为原始地层孔隙压力下的天然气偏差系数,采用Dranchuk-Abu-Kassem方法获取;Z为当前地层孔隙压力下的天然气偏差系数;G为原始地质储量;Δpp为孔隙压力降低值;P为当前底层压力。
7.根据权利要求1所述的确定页岩气藏水平井重复压裂前现地应力场的方法,其特征在于,所述将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加,获取现地应力场分布,包括:其中,σz为深度H处的垂向应力;g为重力加速度;σH和σh分别为最大、最小水平主应力;
σ′H和σ′h分别为现地应力场的最大、最小水平主应力;σx诱导(i)和σy诱导(i)分别为第i段裂缝在x,y方向上的干扰应力;μs为静态泊松比;α为Biot系数;Δpp为孔隙压力降低值;ρs为静态岩体密度。
8.一种确定页岩气藏水平井重复压裂前现地应力场的装置,其特征在于,包括:初始应力场获取模块,用于获取水平井多段簇主裂缝干扰应力场,根据岩石力学动态参数及静态参数,确定原始地应力场,根据页岩气藏水平井生产规律,确定孔隙压力降干扰应力场;
现地应力场获取模块,用于将所述水平井多段簇主裂缝干扰应力场、原始地应力场及孔隙压力降干扰应力场进行叠加,获取现地应力场分布。
9.一种电子设备,其特征在于,包括:
至少一个处理器、至少一个存储器、通信接口和总线;其中,
所述处理器、存储器、通信接口通过所述总线完成相互间的通信;
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令,以执行如权利要求1至7任一项所述的方法。
10.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至7中任一项所述的方法。