CN108717202B - 一种考虑流体温度的页岩气异常地层压力预测方法 - Google Patents
一种考虑流体温度的页岩气异常地层压力预测方法 Download PDFInfo
- Publication number
- CN108717202B CN108717202B CN201810552746.1A CN201810552746A CN108717202B CN 108717202 B CN108717202 B CN 108717202B CN 201810552746 A CN201810552746 A CN 201810552746A CN 108717202 B CN108717202 B CN 108717202B
- Authority
- CN
- China
- Prior art keywords
- formation pressure
- fluid temperature
- shale gas
- well
- depth
- 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
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 67
- 239000012530 fluid Substances 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 18
- 230000002547 anomalous effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 238000005553 drilling Methods 0.000 abstract description 5
- 230000009466 transformation Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 47
- 230000000694 effects Effects 0.000 description 10
- 239000011435 rock Substances 0.000 description 10
- 239000011148 porous material Substances 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000005056 compaction Methods 0.000 description 3
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 229910052901 montmorillonite Inorganic materials 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000002441 X-ray diffraction Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 244000309464 bull Species 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 229910052900 illite Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种考虑流体温度的页岩气异常地层压力预测方法,包括以下步骤:S1在页岩气地区内选择至少四口试验井,由试验井的纵波速度和横波速度计算出试验井的泊松比参数;S2测量每一试验井的地层压力系数、流体温度和深度;S3根据所有试验井的地层压力系数、流体温度和深度确定公式pρ=log eaT+bδ+cH+d中的待定系数a、b、c和d,其中,pρ表示地层压力系数,T表示流体温度,σ表示泊松比参数,H表示深度;S4由预测井的纵波速度和横波速度计算出泊松比参数,测量预测井内流体温度和深度,根据公式pρ=log eaT+bδ+cH+d计算出预测井的地层压力系数。本发明的有益效果:实现考虑流体温度页岩气地层压力高精度预测,为钻井过程中的井壁稳定性、钻井安全和最终储层压裂改造提供指导。
Description
技术领域
本发明涉及页岩气勘探开发领域,尤其涉及一种考虑流体温度的页岩气异常地层压力预测方法。
背景技术
孔隙流体对岩石弹性参数的作用体现在力学效应和物质效应两个方面。孔隙流体的力学效应是岩石发生应变时的应力背景,即孔隙流体压力对骨架应力场的作用;孔隙流体的物质效应则是岩石发生应变时的相对扰动,由于岩石发生应变,骨架与流体之间存在相互作用力,流-固应力场成为一个整体,影响岩石的等效弹性参数。孔隙流体压力,也称地层压力,压力异常反映了地层的烃源条件、封盖条件以及生烃增压、加热增压、烃类液-气转化增压以及粘土矿物转化增压等各种地质过程。目前,常用的异常地层压力预测方法有Eaton法、Fillippone法和等效深度法等,这些方法主要适用于常规油气藏,不适用于页岩气异常地层压力预测。
发明内容
有鉴于此,本发明的实施例提供了一种考虑流体温度的页岩气异常地层压力预测方法。
本发明的实施例提供一种考虑流体温度的页岩气异常地层压力预测方法,包括以下步骤:
S1在页岩气地区内选择至少四口试验井,由每一所述试验井的纵波速度和横波速度计算出每一所述试验井的泊松比参数;
S2测量每一所述试验井的地层压力系数、流体温度和深度;
S3根据所有试验井的地层压力系数、流体温度和深度确定公式pρ=log eaT+bδ+cH+d中的待定系数a、b、c和d,其中,pρ表示地层压力系数,T表示流体温度,σ表示泊松比参数,H表示深度;
S4由所述页岩气地区内预测井的纵波速度和横波速度计算出所述预测井的泊松比参数,测量所述预测井内流体温度和深度,根据公式pρ=log eaT+bδ+cH+d计算出预测井的地层压力系数。
进一步地,步骤S3还包括:
S3.1再次测量所述试验井内流体温度和地层压力系数的实际值,通过公式pρ=log eaT+bδ+cH+d计算出所述试验井的地层压力系数的预测值,并通过Eaton法、Fillippone法和等效深度法中的任意一种方法计算出所述试验井的地层压力系数作为参考值,当预测值比参考值更接近实际值时,判断待定系数a、b、c和d合格,否则重复步骤S2和S3直至预测值比参考值更接近实际值。
进一步地,步骤S3中公式pρ=log eaT+bδ+cH+d中待定系数a、b、c和d确定方法为,由每一所述试验井的地层压力系数、流体温度和深度,根据公式pρ=log eaT+bδ+cH+d使用软件MATLAB进行拟合,分别确定a、b、c和d的值。
进一步地,步骤S1中在页岩气地区内选择至少四口试验井,且每一所述试验井内流体温度不同。
本发明的实施例提供的技术方案带来的有益效果是:本发明考虑流体温度的页岩气异常地层压力预测方法,考虑流体温度因素对页岩气异常地层压力的作用,可以实现页岩气的高精度地层压力的预测,用于页岩气钻前异常地层压力预测,为钻井过程中的井壁稳定性、钻井安全和最终储层压裂改造提供指导。
附图说明
图1是本发明一种考虑流体温度的页岩气异常地层压力预测方法的流程图图;
图2是海湾沿岸某页岩气井中页岩X-衍射反映的蒙脱石成岩轨迹图;
图3是图2中的页岩气井的中新世泥页岩物性交会图;
图4是当比重G和地层压力一定时某页岩气地区的岩石数字岩心模拟结果图;
图5是当比重G和流体弹性参数一定时图4中的页岩气地区的岩石数字岩心模拟结果图;
图6是使用本发明的预测方法对某页岩气井的预测结果平面图;
图7是使用Fillippone法对图6中的页岩气井的预测结果平面图;
图8是使用本发明的预测方法对图6中的页岩气井的预测结果剖面图;
图9是使用Fillippone法对图6中的页岩气井的预测结果剖面图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
在页岩气中,温度与地层史对页岩气异常地层压力起着关键的作用,页岩气上覆地层载荷增加,不仅会降低孔隙度,而且也会排出孔隙水,在这个过程中,温度对地层压实性质的作用,尤其是对页岩中泥质成分的影响非常大。当泥质遭受压实时,它们往往发生埋藏变质作用,例如,蒙脱石到伊利石的转变。这个现象将引起泥质组分的重新排列,岩石骨架应力重新分布,这个结果取决于沉积地层的时间—温度史。这表明地层压实是温度与时间的函数,不能采用单一压实曲线的思想解释岩石压实。
图2是海湾沿岸一页岩气井中页岩X-衍射反映的蒙脱石成岩轨迹,图3是该井的中新世泥页岩物性交会图,表明:声波时差一定,成岩作用将导致体密度增大,埋藏变质过程是动态的,温度与地质史控制这个过程,如果忽略温度及时间的作用,仅使用与体密度成单一关系的速度模型得到的地层压力有效性不高。
图4和图5是某页岩气地区的岩石数字岩心模拟结果,图中结果表明:当该页岩气地区流体的比重G和地层压力一定时,随着流体温度增加,岩石体密度增加,体积模量减小;当比重G和流体弹性参数一定时,地层压力随着流体温度的增加而增加。这说明:页岩气异常地层压力与流体温度和流体的弹性参数有关,当流体确定后,流体的温度和弹性参数变化,流体的压力发生变化,这里的弹性参数即为泊松比参数。
请参考图1,本发明的实施例提供了一种考虑流体温度的页岩气异常地层压力预测方法,包括以下步骤:
S1在页岩气地区内选择至少四口试验井,且每一所述试验井内流体温度不同,由每一所述试验井的纵波速度和横波速度计算出每一所述试验井的泊松比参数,泊松比参数σ的计算公式为其中vp表示纵波速度,vs表示横波速度;
S2测量每一所述试验井的地层压力系数、流体温度和深度;
S3由每一所述试验井的地层压力系数、流体温度和深度,根据公式pρ=log eaT +bδ+cH+d使用软件MATLAB进行拟合,分别确定a、b、c和d的值,其中,pρ表示地层压力系数,T表示流体温度,σ表示泊松比参数,H表示深度,再次测量所述试验井内流体温度和地层压力系数的实际值,通过公式pρ=log eaT+bδ+cH+d计算出所述试验井的地层压力系数的预测值,并通过Eaton法、Fillippone法和等效深度法中的任意一种方法计算出所述试验井的地层压力系数作为参考值,当预测值比参考值更接近实际值时,判断待定系数a、b、c和d合格,否则重复步骤S2和S3直至预测值比参考值更接近实际值;
S4由所述页岩气地区内预测井的纵波速度和横波速度计算出所述预测井的泊松比参数,测量所述预测井内流体温度和深度,根据公式pρ=log eaT+bδ+cH+d计算出预测井的地层压力系数。
下面举例对上述方法进行验证,以某页岩气井为测试对象,分别采用本发明的地层压力系数预测方法和Fillippone法对该页岩气井进行地层压力系数进行计算。
图6表示使用本发明的预测方法对该页岩气井的预测结果平面图,图7表示使用Fillippone法对该页岩气井的预测结果平面图,对比可知:图7中平面变化较大,不光滑,有些“牛眼”出现,而图6中平面较光滑,变化平缓,因此使用本发明预测方法的预测结果比较合理。
图8表示使用本发明的预测方法对该页岩气井的预测结果剖面图,图9表示使用Fillippone法对该页岩气井的预测结果剖面图,对比可知:图9中剖面上压力有突变的情况,出现预测异常,不符合实际情况,图8中剖面不存在突变,预测结果合理,更为可信。
因此本发明考虑流体温度的页岩气异常地层压力预测方法有效可行,考虑流体温度因素对页岩气异常地层压力的作用,可以实现页岩气的高精度地层压力的预测。
在本文中,所涉及的前、后、上、下等方位词是以附图中零部件位于图中以及零部件相互之间的位置来定义的,只是为了表达技术方案的清楚及方便。应当理解,所述方位词的使用不应限制本申请请求保护的范围。
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种考虑流体温度的页岩气异常地层压力预测方法,其特征在于,包括以下步骤:
S1在页岩气地区内选择至少四口试验井,由每一所述试验井的纵波速度和横波速度计算出每一所述试验井的泊松比参数;
S2测量每一所述试验井的地层压力系数、流体温度和深度;
S3根据所有试验井的地层压力系数、流体温度和深度确定公式pρ=log eaT+bδ+cH+d中的待定系数a、b、c和d,其中,pρ表示地层压力系数,T表示流体温度,σ表示泊松比参数,H表示深度;
S4由所述页岩气地区内预测井的纵波速度和横波速度计算出所述预测井的泊松比参数,测量所述预测井内流体温度和深度,根据公式pρ=log eaT+bδ+cH+d计算出预测井的地层压力系数。
2.如权利要求1所述的一种考虑流体温度的页岩气异常地层压力预测方法,其特征在于,步骤S3还包括:
S3.1再次测量所述试验井内流体温度和地层压力系数的实际值,通过公式pρ=log eaT +bδ+cH+d计算出所述试验井的地层压力系数的预测值,并通过Eaton法、Fillippone法和等效深度法中的任意一种方法计算出所述试验井的地层压力系数作为参考值,当预测值比参考值更接近实际值时,判断待定系数a、b、c和d合格,否则重复步骤S2和S3直至预测值比参考值更接近实际值。
4.如权利要求1所述的一种考虑流体温度的页岩气异常地层压力预测方法,其特征在于:步骤S3中公式pρ=log eaT+bδ+cH+d中待定系数a、b、c和d确定方法为,由每一所述试验井的地层压力系数、流体温度和深度,根据公式pρ=log eaT+bδ+cH+d使用软件MATLAB进行拟合,分别确定a、b、c和d的值。
5.如权利要求1所述的一种考虑流体温度的页岩气异常地层压力预测方法,其特征在于:步骤S1中在页岩气地区内选择至少四口试验井,且每一所述试验井内流体温度不同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810552746.1A CN108717202B (zh) | 2018-05-31 | 2018-05-31 | 一种考虑流体温度的页岩气异常地层压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810552746.1A CN108717202B (zh) | 2018-05-31 | 2018-05-31 | 一种考虑流体温度的页岩气异常地层压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108717202A CN108717202A (zh) | 2018-10-30 |
CN108717202B true CN108717202B (zh) | 2020-01-17 |
Family
ID=63912663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810552746.1A Active CN108717202B (zh) | 2018-05-31 | 2018-05-31 | 一种考虑流体温度的页岩气异常地层压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108717202B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110471129A (zh) * | 2019-09-20 | 2019-11-19 | 西南石油大学 | 一种深层页岩高温高压下的各向异性岩石物理建模方法 |
CN110671095B (zh) * | 2019-09-23 | 2021-03-26 | 中国地质大学(武汉) | 一种地质地层压力智能随钻软测量方法 |
CN110821480B (zh) * | 2019-10-14 | 2023-03-17 | 大庆油田有限责任公司 | 一种快速预测待钻井地层压力系数的图版方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107121703A (zh) * | 2017-06-02 | 2017-09-01 | 成都理工大学 | 一种基于三维地震资料的页岩气地层的地应力预测方法 |
CN107703543A (zh) * | 2017-09-11 | 2018-02-16 | 武汉拓盟能源科技有限公司 | 一种基于广义泊松比的地层孔隙压力预测方法及系统 |
CN107728205A (zh) * | 2016-08-12 | 2018-02-23 | 中国石油化工股份有限公司 | 一种地层压力预测方法 |
CN107843927A (zh) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | 基于井震联合速度的页岩地层压力预测方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105927218B (zh) * | 2016-05-20 | 2019-05-07 | 中国石油大学(北京) | 一种陆相泥页岩储层含气量预测方法及装置 |
CN106251232B (zh) * | 2016-08-12 | 2019-09-06 | 中国石油天然气股份有限公司 | 确定页岩含气量的方法和装置 |
CN107797139B (zh) * | 2016-08-31 | 2019-07-23 | 中国石油化工股份有限公司 | 页岩储层游离气含气量地震预测方法及系统 |
-
2018
- 2018-05-31 CN CN201810552746.1A patent/CN108717202B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728205A (zh) * | 2016-08-12 | 2018-02-23 | 中国石油化工股份有限公司 | 一种地层压力预测方法 |
CN107843927A (zh) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | 基于井震联合速度的页岩地层压力预测方法及装置 |
CN107121703A (zh) * | 2017-06-02 | 2017-09-01 | 成都理工大学 | 一种基于三维地震资料的页岩气地层的地应力预测方法 |
CN107703543A (zh) * | 2017-09-11 | 2018-02-16 | 武汉拓盟能源科技有限公司 | 一种基于广义泊松比的地层孔隙压力预测方法及系统 |
Non-Patent Citations (3)
Title |
---|
"Deepwater geohazard prediction using prestack inversion of large offset P-wave data and rock model";N.C. Dutta;《The Leading Eage》;20020228;第21卷(第2期);第193页右栏最后1段-第188页左栏第1段 * |
"弹性参数联合预测地层压力的方法及应用";王斌 等;《2013年物探技术研讨会文集》;20140531;第613-616页 * |
"碳酸盐岩地层压力预测方法研究";朱伟 等;《CPS/SEG北京2009国际地球物理会议论文集》;20100531;第1170(ID)页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108717202A (zh) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eaton | The effect of overburden stress on geopressure prediction from well logs | |
US8423337B2 (en) | Method for multi-scale geomechanical model analysis by computer simulation | |
CN108717202B (zh) | 一种考虑流体温度的页岩气异常地层压力预测方法 | |
US20120239363A1 (en) | Production estimation in subterranean formations | |
CN109856674B (zh) | 工程甜点评测方法 | |
Khamehchi et al. | Sand production prediction using ratio of shear modulus to bulk compressibility (case study) | |
WO2016118802A1 (en) | Control system and method of flowback operations for shale reservoirs | |
El Ouahed et al. | Application of artificial intelligence to characterize naturally fractured zones in Hassi Messaoud Oil Field, Algeria | |
CN101025084A (zh) | 一种随钻预测钻头底下地层孔隙压力的方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
Bishop | Calculated compaction states of thick abnormally pressured shales | |
Langenbruch et al. | Quantitative analysis of rock stress heterogeneity: Implications for the seismogenesis of fluid-injection-induced seismicity | |
Ranjbar-Karami et al. | A modified fuzzy inference system for estimation of the static rock elastic properties: A case study from the Kangan and Dalan gas reservoirs, South Pars gas field, the Persian Gulf | |
CN108304959B (zh) | 提高地层流体压力预测精度的方法 | |
CN107703543B (zh) | 一种基于广义泊松比的地层孔隙压力预测方法及系统 | |
Ugwu et al. | Integrated approach to geopressure detection in the X-field, Onshore Niger Delta | |
CN101936157A (zh) | 一种利用测井资料检测高压盐水层孔隙压力的方法 | |
CN113221347A (zh) | 一种井壁稳定性钻井优化方法、装置及设备 | |
EP3289180A1 (en) | Determination of horizontal constraints in subsoil | |
Prankada et al. | Analysis of wellbore stability by pore pressure prediction using seismic velocity | |
Ayodele et al. | Pore pressure prediction of some selected wells from the Southern Pletmos Basin, offshore South Africa | |
CN105301642A (zh) | 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置 | |
Michelena et al. | Seismic, geologic, geomechanics, and dynamic constraints in flow models of unconventional fractured reservoirs: Example from a south Texas field | |
Gong et al. | Challenges in pore pressure prediction for unconventional petroleum systems | |
CN112651098A (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 |