CN112257304B - 页岩地层直井水力裂缝高度预测方法 - Google Patents
页岩地层直井水力裂缝高度预测方法 Download PDFInfo
- Publication number
- CN112257304B CN112257304B CN202010968797.XA CN202010968797A CN112257304B CN 112257304 B CN112257304 B CN 112257304B CN 202010968797 A CN202010968797 A CN 202010968797A CN 112257304 B CN112257304 B CN 112257304B
- Authority
- CN
- China
- Prior art keywords
- shale
- hydraulic fracture
- height
- vertical
- shale formation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 81
- 238000005070 sampling Methods 0.000 claims abstract description 14
- 230000035699 permeability Effects 0.000 claims description 29
- 239000011435 rock Substances 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000002790 cross-validation Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 239000003079 shale oil Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000009933 burial Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geochemistry & Mineralogy (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种页岩地层直井水力裂缝高度预测方法。页岩地层直井水力裂缝高度预测方法包括:S1、建立页岩地层水力裂缝高度的影响指标体系,影响指标体系包括多个影响指标;S2、建立页岩地层水力压裂模型;S3、采用扩展有限单元法计算页岩地层水力裂缝高度;S4、在各影响指标的取值范围内开展拉丁超立方抽样,并计算各影响指标的组合条件下的页岩地层水力裂缝高度;S5、基于多元自适应回归样条法建立各影响指标关于页岩地层水力裂缝高度的预测模型,预测页岩地层水力裂缝高度。本发明解决了现有技术中页岩地层直井水力裂缝高度预测方法建模繁琐、数值计算机时多的问题。
Description
技术领域
本发明涉及能源开采技术领域,具体而言,涉及一种页岩地层直井水力裂缝高度预测方法。
背景技术
页岩储层作为我国重要的能源组成部分,具有纵向层数多、单层厚度薄、平面分布广等特征。受页岩油丰度低、储层渗流能力差等地质特征制约,页岩储层单井一般无自然产能或自然产能低于工业油流下限,且储层能量衰减快、补给困难。以地质工程一体化理念为指导,针对页岩油储层薄互层发育的特点,通过实现“甜点区”的科学***评价,通过直井分段体积压裂实现“人造高渗区、重构渗流场”等技术体系,对破解页岩油经济开采困局、实现工业效益开采具有重要意义。页岩储层中的富含有机质岩石在成藏过程中,在埋深、压力、温度、应力等复杂外部环境下,呈现出明显的层理结构特征。在层理发育程度与其力学性质影响下,页岩的抗压、抗拉、硬度等力学特征参数等存在明显的各向异性,层理面往往为地层中的力学薄弱面,受到外力时易先于页岩本体而破坏。在页岩储层开展水力压裂施工过程中,水力裂缝受层理影响显著,裂缝形态与常规地层存在较大差异,对储层有效改造体积影响巨大,因此,急需开展人工裂缝缝高预测,优化压裂工艺及施工参数,达到提高单井产量的目的。
并且,在现有技术中还提供了一种利用测井资料预测判断水力压裂裂缝高度的方法,具体包括利用阵列声波资料以及过套管交叉偶极横波测井资料可准确地检测和预测水力压裂裂缝的缝高,从而定量地对压裂施工作业进行指导和评估。相比于常规测井方法,该论文中的两种方法具有简单、直观、快速、不受井时限制以及便于推广等优势,但在没有加测阵列声波资料的井,该方法的适用性受到限制。
中国专利CN108280275A公开了一种致密砂岩水力压裂缝高预测方法,该方法包括建立由隔层-储层-隔层组成的三维水力压裂裂缝扩展有限元模型;基于水力压裂扩展有限元模型分析致密砂岩气藏水力压裂裂缝纵向扩展的影响因素,进而分析得出裂缝内净压力与储隔层地应力差比值和临界隔储层厚度比的对应数据关系,用以指导现场控制缝高压裂时施工参数的制定。但该方法需要通过建立页岩地层模型和流固耦合数值方程,进而进行扩展有限元数值模拟,以确定层理发育的页岩地层水力裂缝高度,建模过程较繁琐,数值计算所需机时较多。
针对现有技术存在的建模繁琐、数值计算机时较多等问题,需要寻找一种水力裂缝高度的预测方法,使得该方法能够高效而准确地对各工况下页岩地层水力裂缝高度进行预测。
因此,现有技术中存在页岩地层直井水力裂缝高度预测方法建模繁琐、数值计算机时多的问题。
发明内容
本发明的主要目的在于提供一种页岩地层直井水力裂缝高度预测方法,以解决现有技术中页岩地层直井水力裂缝高度预测方法建模繁琐、数值计算机时多的问题。
为了实现上述目的,根据本发明的一个方面,提供了一种页岩地层直井水力裂缝高度预测方法,包括:S1、建立页岩地层水力裂缝高度的影响指标体系,影响指标体系包括多个影响指标;S2、建立页岩地层水力压裂模型;S3、采用扩展有限单元法计算页岩地层水力裂缝高度;S4、在各影响指标的取值范围内开展拉丁超立方抽样,并计算各影响指标的组合条件下的页岩地层水力裂缝高度;S5、基于多元自适应回归样条法建立各影响指标关于页岩地层水力裂缝高度的预测模型,预测页岩地层水力裂缝高度。
进一步地,在步骤S2中,根据目标地层埋深和应力条件,采用横观各向同性材料描述页岩地层的岩石力学特性及渗透特性,建立包含射孔及步骤S1中页岩地层水力裂缝高度的影响指标体系的力学模型。
进一步地,步骤S2包括:S21、根据目标井深位置和应力状态建立力学模型,并确定射孔方位;S22、根据应力边界条件,开展地应力平衡,计算初始井周应力场;S23、采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,定义页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率以及页岩孔隙比。
进一步地,步骤S3包括:S31、根据设计的射孔方位,在地层中设置起裂方位;S32、在射孔位置施加设计的压裂液排量,模拟水力压裂过程;S33、裂缝扩展至设计压裂时间,裂缝停止扩展,统计水力裂缝高度。
进一步地,步骤S4包括:S41、统计各影响指标的取值范围;S42、在步骤S41所统计各影响指标取值范围内开展拉丁超立方抽样,获得N组影响指标的组合条件,其中100≤N≤5000;S43、使用扩展有限单元法,计算步骤S42中各影响指标的组合条件下的页岩地层水力裂缝高度。
进一步地,步骤S5包括:S51、以步骤S42中各影响指标组合为输入矩阵,以步骤S43中各影响指标组合条件下的水力裂缝高度为输出矩阵,建立多元自适应回归样条预测模型;S52、采用交叉验证法确定步骤S51所建立的预测模型所需的基函数数量M;S53、建立基函数数量M条件下的多元自适应回归样条预测模型,预测页岩地层水力裂缝高度。
进一步地,步骤S1中页岩地层水力裂缝高度的影响指标体系的各影响指标根据地层应力状态、页岩岩石力学参数及施工参数确定。
进一步地,步骤S1中页岩地层水力裂缝高度的影响指标体系的各影响指标包括:射孔方位、垂向-水平主应力差、页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率,页岩孔隙比,压裂液排量及压裂时间。
进一步地,射孔方位的取值范围为0°-90°、垂向-水平主应力差的取值范围为8-20MPa、页岩垂向杨氏模量的取值范围为18-20GPa,页岩垂向泊松比的取值范围为0.09-0.11,页岩垂向渗透率的取值范围为0.003-0.005mD。
进一步地,页岩横向杨氏模量的取值范围为10-12GPa,页岩横向泊松比的取值范围为0.09-0.11,页岩横向渗透率的取值范围为0.004-0.006mD,页岩孔隙比的取值范围为0.1-0.12,压裂液排量的取值范围为6-9m3/min,压裂时间的取值范围为180-240min。
应用本发明的技术方案,本申请中的页岩地层直井水力裂缝高度预测方法包括:S1、建立页岩地层水力裂缝高度的影响指标体系,影响指标体系包括多个影响指标;S2、建立页岩地层水力压裂模型;S3、采用扩展有限单元法计算页岩地层水力裂缝高度;S4、在各影响指标的取值范围内开展拉丁超立方抽样,并计算各影响指标的组合条件下的页岩地层水力裂缝高度;S5、基于多元自适应回归样条法建立各影响指标关于页岩地层水力裂缝高度的预测模型,预测页岩地层水力裂缝高度。
本发明采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,避免了传统方法对层理单独建模,导致单元数量过于庞大的弊端。并且,本发明在模型中利用扩展有限元方法判定三维裂缝的扩展方向,并统计水力裂缝高度,采用拉丁超立方抽样结合多元自适应回归样条法建立各影响指标关于水力裂缝高度的预测模型,避免了传统响应面法需要假设基函数阶数的问题,实现页岩地层水力裂缝高度的高效、准确预测。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1示出了根据本发明的一个具体实施例的页岩地层直井水力裂缝高度预测方法流程图;
图2示出了本发明的一个具体实施例中页岩地层直井水力裂缝高度预测方法的力学模型几何及边界条件图;
图3示出了本发明的一个具体实施例中页岩地层直井水力裂缝高度预测方法的压裂200min后模型中的孔隙流体压力分布图;
图4示出了本发明的一个具体实施例中页岩地层直井水力裂缝高度预测方法的GCV随基函数数量M变化趋势图;
图5示出了本发明的一个具体实施例中页岩地层直井水力裂缝高度预测方法的模型预测结果与扩展有限元计算结果对比图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
需要指出的是,除非另有指明,本申请使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
在本发明中,在未作相反说明的情况下,使用的方位词如“上、下、顶、底”通常是针对附图所示的方向而言的,或者是针对部件本身在竖直、垂直或重力方向上而言的;同样地,为便于理解和描述,“内、外”是指相对于各部件本身的轮廓的内、外,但上述方位词并不用于限制本发明。
为了解决现有技术中页岩地层直井水力裂缝高度预测方法建模繁琐、数值计算机时多的问题,本申请提供了一种页岩地层直井水力裂缝高度预测方法。
本申请中的页岩地层直井水力裂缝高度预测方法包括:S1、建立页岩地层水力裂缝高度的影响指标体系,影响指标体系包括多个影响指标;S2、建立页岩地层水力压裂模型;S3、采用扩展有限单元法计算页岩地层水力裂缝高度;S4、在各影响指标的取值范围内开展拉丁超立方抽样,并计算各影响指标的组合条件下的页岩地层水力裂缝高度;S5、基于多元自适应回归样条法建立各影响指标关于页岩地层水力裂缝高度的预测模型,预测页岩地层水力裂缝高度。
本发明采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,避免了传统方法对层理单独建模,导致单元数量过于庞大的弊端。并且,本发明在模型中利用扩展有限元方法判定三维裂缝的扩展方向,并统计水力裂缝高度,采用拉丁超立方抽样结合多元自适应回归样条法建立各影响指标关于水力裂缝高度的预测模型,避免了传统响应面法需要假设基函数阶数的问题,实现页岩地层水力裂缝高度的高效、准确预测。
如图1所示,为本申请的一个具体实施例中页岩地层直井水力裂缝高度预测方法的流程图,主要包括:建立页岩地层水力裂缝高度的影响指标体系、采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,建立页岩地层水利压裂模型、采用扩展有限单元法计算页岩地层水力裂缝高度、在各参数取值范围内展开超拉丁立方抽样,并计算各影响指标组合条件下的水力裂缝高度、基于多元自适应回归样条法建立各影响指标关于水力裂缝高度的预测模型。
具体地,在步骤S2中,根据目标地层埋深和应力条件,采用横观各向同性材料描述页岩地层的岩石力学特性及渗透特性,建立包含射孔及步骤S1中页岩地层水力裂缝高度的影响指标体系的力学模型。
如图2所示,步骤S2包括:S21、根据目标井深位置和应力状态建立力学模型,并确定射孔方位;S22、根据应力边界条件,开展地应力平衡,计算初始井周应力场;S23、采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,定义页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率以及页岩孔隙比。并且,各个参数的取值见下表1。
表1页岩地层水力裂缝高度影响指标取值
具体地,步骤S3包括:S31、根据设计的射孔方位,在地层中设置起裂方位;S32、在射孔位置施加设计的压裂液排量,模拟水力压裂过程;S33、裂缝扩展至设计压裂时间,裂缝停止扩展,统计此时水力裂缝高度H,其中模拟孔隙流体水利压力分布如图3所示。
具体地,步骤S4包括:S41、统计各影响指标的取值范围,各影响指标的取值范围如表2所示。
表2.水力裂缝高度影响指标取值范围
S42、在步骤S41所统计各影响指标取值范围内开展拉丁超立方抽样,获得N组影响指标的组合条件,如表3所示,其中100≤N≤5000;S43、使用扩展有限单元法,计算步骤S42中各影响指标的组合条件下的页岩地层水力裂缝高度,如表3中H列所示。
表3.拉丁超立方抽样组合及相应裂缝高度
具体地,步骤S5包括:S51、以步骤S42中各影响指标组合为输入矩阵,以步骤S43中各影响指标组合条件下的水力裂缝高度(表3中H列)为输出矩阵,建立多元自适应回归样条预测模型,在MATLAB中实现语句为:
params=aresparams2('maxFuncs',40,'maxInteractions',2);
[model,~,resultsEval]=aresbuild(X,Y,params)
figure;
hold on;grid on;box on;
h(1)=plot(resultsEval.MSE,'Color',[0 0.447 0.741]);
h(2)=plot(resultsEval.GCV,'Color',[0.741 0 0.447]);
numBF=numel(model.coefs);
h(3)=plot([numBF numBF],get(gca,'ylim'),'--k');
xlabel('Number of basis functions');
ylabel('MSE,GCV');
legend(h,'MSE','GCV','Selected model');
areseq(model,5);
S52、采用交叉验证法(GCV)确定步骤S51所建立的预测模型所需的基函数数量M,结果如图4所示,当M=24时,GCV达到稳定值,多元自适应回归样条预测模型最优基函数数量M=24;S53、建立基函数数量M条件下的多元自适应回归样条预测模型,预测页岩地层水力裂缝高度。
预测模型可表示为:
BF1=C(x4|+1,0.039254,0.06474,0.5299)
BF2=C(x2|+1,0.34774,0.52205,0.67057)
BF3=C(x2|-1,0.34774,0.52205,0.67057)
BF4=C(x1|-1,0.21355,0.4267,0.44679)
BF5=C(x5|+1,0.34577,0.68422,0.84085)
BF6=C(x5|-1,0.34577,0.68422,0.84085)
BF7=C(x3|+1,0.35546,0.55109,0.6691)
BF8=BF2*C(x1|+1,0.49312,0.51935,0.75663)
BF9=BF2*C(x1|-1,0.49312,0.51935,0.75663)
BF10=BF3*C(x1|-1,0.44679,0.46688,0.49312)
BF11=C(x3|+1,0.081473,0.15983,0.35546)
BF12=C(x1|+1,0.21355,0.4267,0.44679)*C(x2|+1,0.67057,0.81909,0.90792)
BF13=C(x1|+1,0.21355,0.4267,0.44679)*C(x2|-1,0.67057,0.81909,0.90792)
BF14=C(x3|+1,0.6691,0.78712,0.88823)
BF15=C(x3|-1,0.6691,0.78712,0.88823)
BF16=C(x1|+1,0.21355,0.4267,0.44679)*C(x2|+1,0.088224,0.17343,0.34774)
BF17=C(x1|+1,0.21355,0.4267,0.44679)*C(x2|-1,0.088224,0.17343,0.34774)
y=-0.92075+9.8169*BF1+7.5374*BF2-13.917*BF3-15.493*BF4+4.6498*BF5-5.1358*BF6+14.101*BF7-58.297*BF8-12.713*BF9+26.732*BF10+16.505*BF11-45.16*BF12+10.476*BF13-13.936*BF14+22.166*BF15+17.011*BF16-46.639*BF17
采用上述多元自适应回归样条预测模型可实现页岩地层水力裂缝高度的高效预测。开展10种工况下的页岩地层水力裂缝高度预测,如表4所示,并与扩展有限元计算结果对比,如图5所示。
表4页岩地层水力裂缝高度预测结果
由表4和图5可知,页岩地层水力裂缝高度预测结果与计算结果基本相符,验证了本发明中预测方法的准确性和可靠性。
具体地,步骤S1中页岩地层水力裂缝高度的影响指标体系的各影响指标根据地层应力状态、页岩岩石力学参数及施工参数确定。
具体地,步骤S1中页岩地层水力裂缝高度的影响指标体系的各影响指标包括:射孔方位、垂向-水平主应力差、页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率,页岩孔隙比,压裂液排量及压裂时间。
具体地,射孔方位的取值范围为0°-90°、垂向-水平主应力差的取值范围为8-20MPa、页岩垂向杨氏模量的取值范围为18-20GPa,页岩垂向泊松比的取值范围为0.09-0.11,页岩垂向渗透率的取值范围为0.003-0.005mD。页岩横向杨氏模量的取值范围为10-12GPa,页岩横向泊松比的取值范围为0.09-0.11,页岩横向渗透率的取值范围为0.004-0.006mD,页岩孔隙比的取值范围为0.1-0.12,压裂液排量的取值范围为6-9m3/min,压裂时间的取值范围为180-240min。
从以上的描述中,可以看出,本发明上述的实施例实现了如下技术效果:
1.本发明采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,避免了传统方法对层理单独建模,导致单元数量过于庞大的弊端。
2.本发明在模型中利用扩展有限元方法判定三维裂缝的扩展方向,并统计水力裂缝高度,采用拉丁超立方抽样结合多元自适应回归样条法建立各影响指标关于水力裂缝高度的预测模型,避免了传统响应面法需要假设基函数阶数的问题,实现页岩地层水力裂缝高度的高效、准确预测。
显然,上述所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、工作、器件、组件和/或它们的组合。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施方式能够以除了在这里图示或描述的那些以外的顺序实施。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种页岩地层直井水力裂缝高度预测方法,其特征在于,包括:
S1、建立页岩地层水力裂缝高度的影响指标体系,所述影响指标体系包括多个影响指标;
S2、建立页岩地层水力压裂模型;
S3、采用扩展有限单元法计算所述页岩地层水力裂缝高度;
S4、在各所述影响指标的取值范围内开展拉丁超立方抽样,并计算各所述影响指标的组合条件下的所述页岩地层水力裂缝高度;
S5、基于多元自适应回归样条法建立开展拉丁超立方抽样后的各影响指标的组合关于所述页岩地层水力裂缝高度的预测模型,预测所述页岩地层水力裂缝高度。
2.根据权利要求1所述的页岩地层直井水力裂缝高度预测方法,其特征在于,在所述步骤S2中,根据目标地层埋深和应力条件,采用横观各向同性材料描述页岩地层的岩石力学特性及渗透特性,建立包含射孔及所述步骤S1中所述页岩地层水力裂缝高度的影响指标体系的力学模型。
3.根据权利要求2所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S2包括:
S21、根据目标井深位置和应力状态建立力学模型,并确定射孔方位;
S22、根据应力边界条件,开展地应力平衡,计算初始井周应力场;
S23、采用横观各向同性材料描述页岩的岩石力学特性及渗透特性,定义页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率以及页岩孔隙比。
4.根据权利要求1所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S3包括:
S31、根据设计的射孔方位,在地层中设置起裂方位;
S32、在射孔位置施加设计的压裂液排量,模拟水力压裂过程;
S33、裂缝扩展至设计压裂时间,裂缝停止扩展,统计水力裂缝高度。
5.根据权利要求1所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S4包括:
S41、统计各所述影响指标的取值范围;
S42、在所述步骤S41所统计各影响指标取值范围内开展拉丁超立方抽样,获得N组所述影响指标的组合条件,其中100≤N≤5000;
S43、使用扩展有限单元法,计算所述步骤S42中各所述影响指标的组合条件下的所述页岩地层水力裂缝高度。
6.根据权利要求5所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S5包括:
S51、以所述步骤S42中各所述影响指标组合为输入矩阵,以所述步骤S43中各所述影响指标组合条件下的水力裂缝高度为输出矩阵,建立多元自适应回归样条预测模型;
S52、采用交叉验证法确定所述步骤S51所建立的预测模型所需的基函数数量M;
S53、建立所述基函数数量M条件下的多元自适应回归样条预测模型,预测所述页岩地层水力裂缝高度。
7.根据权利要求1至6中任一项所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S1中所述页岩地层水力裂缝高度的所述影响指标体系的各所述影响指标根据地层应力状态、页岩岩石力学参数及施工参数确定。
8.根据权利要求1至6中任一项所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述步骤S1中所述页岩地层水力裂缝高度的所述影响指标体系的各所述影响指标包括:射孔方位、垂向-水平主应力差、页岩垂向杨氏模量,页岩垂向泊松比,页岩垂向渗透率,页岩横向杨氏模量,页岩横向泊松比,页岩横向渗透率,页岩孔隙比,压裂液排量及压裂时间。
9.根据权利要求8所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述射孔方位的取值范围为0°-90°、所述垂向-水平主应力差的取值范围为8-20MPa、所述页岩垂向杨氏模量的取值范围为18-20GPa,所述页岩垂向泊松比的取值范围为0.09-0.11,所述页岩垂向渗透率的取值范围为0.003-0.005mD。
10.根据权利要求8所述的页岩地层直井水力裂缝高度预测方法,其特征在于,所述页岩横向杨氏模量的取值范围为10-12GPa,所述页岩横向泊松比的取值范围为0.09-0.11,所述页岩横向渗透率的取值范围为0.004-0.006mD,所述页岩孔隙比的取值范围为0.1-0.12,所述压裂液排量的取值范围为6-9m3/min,所述压裂时间的取值范围为180-240min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010968797.XA CN112257304B (zh) | 2020-09-15 | 2020-09-15 | 页岩地层直井水力裂缝高度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010968797.XA CN112257304B (zh) | 2020-09-15 | 2020-09-15 | 页岩地层直井水力裂缝高度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112257304A CN112257304A (zh) | 2021-01-22 |
CN112257304B true CN112257304B (zh) | 2021-11-30 |
Family
ID=74232582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010968797.XA Active CN112257304B (zh) | 2020-09-15 | 2020-09-15 | 页岩地层直井水力裂缝高度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112257304B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115126459B (zh) * | 2021-03-26 | 2024-05-28 | 中国石油天然气股份有限公司 | 水力裂缝高度的处理方法和装置 |
CN113216947B (zh) * | 2021-05-17 | 2023-01-13 | 中国石油大学(华东) | 基于监测井分布式光纤应变监测的水平井压裂过程裂缝高度确定方法 |
CN115577645B (zh) * | 2022-12-08 | 2023-04-18 | 中国石油大学(华东) | 一种燃爆压裂破裂范围预测模型的构建方法及预测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248494B (zh) * | 2016-08-29 | 2019-03-01 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种用于页岩气井储层脆性综合评价的方法 |
CN108468538B (zh) * | 2018-03-23 | 2020-04-21 | 西南石油大学 | 一种页岩水压裂缝扩展预测方法 |
CN110348032B (zh) * | 2018-04-08 | 2021-09-14 | 中国石油化工股份有限公司 | 层理发育页岩地层水力裂缝高度数值模拟方法 |
-
2020
- 2020-09-15 CN CN202010968797.XA patent/CN112257304B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112257304A (zh) | 2021-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112257304B (zh) | 页岩地层直井水力裂缝高度预测方法 | |
CN108009705A (zh) | 一种基于支持向量机技术的页岩储层可压性评价方法 | |
CN105201479B (zh) | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 | |
CN106869896B (zh) | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 | |
CN109902918B (zh) | 一种页岩压裂自支撑裂缝缝网渗透率计算方法 | |
CN105760564B (zh) | 分析油层套管失效的方法和装置 | |
CN105046006A (zh) | 一种页岩气藏水平井多段压裂产能预测方法及装置 | |
CN104863560A (zh) | 用于页岩气开采的宽网压裂方法 | |
CN102182453B (zh) | 井壁坍塌分析方法 | |
CN112727534B (zh) | 基于真三向应力和渗透率动态变化的瓦斯钻孔布置方法 | |
CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
CN112682016B (zh) | 油气田薄互层储层穿层压裂参数的确定方法及压裂方法 | |
CN110501758B (zh) | 一种砂砾岩储层纵向连续脆性指数预测方法 | |
WO2017031615A1 (zh) | 一种构建双溶腔盐穴储库地面沉降预测模型的方法 | |
CN107290799B (zh) | 一种岩石可压性的确定方法 | |
Liu et al. | Well type and pattern optimization method based on fine numerical simulation in coal-bed methane reservoir | |
CN108984877A (zh) | 一种致密储层复杂压裂裂缝形成条件的评估方法 | |
CN109902422B (zh) | 一种井眼轨迹选取方法及装置 | |
CN113255092A (zh) | 体积压裂水平井特征点拟合反演方法、电子设备及介质 | |
CN107451671B (zh) | 用于预测页岩地层压裂后初始产能的方法及*** | |
CN112100707A (zh) | 一种穿层压裂图版构建方法 | |
CN112012710A (zh) | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 | |
El Sgher et al. | The Impact of Stress Shadowing and Geomechanical Effects on Gas Production From Marcellus Shale | |
CN115324557A (zh) | 基于多因素分析预测压裂诱发套管变形风险程度的方法 | |
CN115270533A (zh) | 重复压裂设计方法、装置、存储介质及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |