CN107066718A - 一种四维动态地应力模拟方法 - Google Patents
一种四维动态地应力模拟方法 Download PDFInfo
- Publication number
- CN107066718A CN107066718A CN201710217338.6A CN201710217338A CN107066718A CN 107066718 A CN107066718 A CN 107066718A CN 201710217338 A CN201710217338 A CN 201710217338A CN 107066718 A CN107066718 A CN 107066718A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- model
- stress
- dynamic
- field
- 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.)
- Granted
Links
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种四维动态地应力模拟方法,它包括以下步骤:S1、建立三维地质模型;S2、建立三维油藏模型并利用生产/注入动态参数预测不同时期的三维孔隙压力场和温度场;S3、建立三维地应力模型;S4、形成初始三维地应力场;S5、根据三维初始地应力场模型,以动态三维孔隙压力场和温度场作为边界条件,建立四维动态地应力模型;S6、进行渗流‑应力耦合求解,分析计算动态地应力及孔弹性参数并用于指导钻完井工程实践。本发明的有益效果是:能够准确反应油气开发过程中动态地应力及孔弹性参数变化情况,为老井重复压裂和新井钻完井作业提供指导。
Description
技术领域
本发明涉及油气资源开发领域,特别是一种四维动态地应力模拟方法。
背景技术
在油气资源开发过程中,地应力模型以三维静态模型为主,着眼于描述当前的地应力分布状态。然而,对于油气开发过程中的压实和沉降、盐岩地层压力分布、盐丘或断层附近的孔隙压力衰竭或应力反转以及各向异性特征明显的储层等,常规的三维静态地应力模型无法根据过往的生产或注入动态情况准确描述当前的地应力状态或预测未来一定时间内地应力的变化情况。因此,根据生产或注入动态情况建立高精度的四维动态地应力模型能够更为准确描述及预测地应力状态及变化情况,并用于解决油气开发过程中的各项钻完井问题。
发明内容
本发明能够模拟分析及预测油气资源在不同生产或注入时期内不同位置处地应力及储层和岩石力学参数的变化情况,提供一种四维动态地应力模拟方法。
本发明的目的通过以下技术方案来实现:一种四维动态地应力模拟方法,它包括以下步骤:
S1、建立三维地质模型,该三维地质模型中至少应包含三维离散的地层厚度、储层物性参数及岩石力学参数;
S2、根据三维地质模型建立三维油藏模型,并利用生产或注入动态参数预测不同时期的三维孔隙压力场和温度场,同时建立非储层段模型并计算其属性参数;
S3、根据三维地质模型,建立带有储层物性及岩石力学属性的三维地应力模型;
S4、利用应力平衡法及单井地应力计算结果横向插值校正,形成初始三维地应力场;
S5、以油藏数值模拟结果导出的三维孔隙压力场和温度场作为不同计算时间步的边界条件,以应力初始化后的三维地应力模型为初始模型,建立四维动态地应力模型;
S6、对四维动态地应力模型进行渗流-应力耦合迭代计算,得到动态地应力及孔弹性参数计算结果,分析地应力、地层位移、体积应变率、孔隙比、渗透率以及孔隙压力动态孔弹性参数的变化情况,以用于指导钻完井工程实践。
本发明具有以下优点:(1)本发明通过考虑生产或注入动态产生,建立了渗流-应力耦合的四维动态地应力模型,克服了三维静态地应力模型无法准确反应油气开发过程中地应力及储层参数变化的问题。(2)本发明能够准确描述油气储层孔弹性参数的非均质性和各向异性,更大程度上地反映储层的真实状态。(3)本发明能够广泛适用于各个类型的油藏,为我国油气资源开发甚至是其它相关地质资源的安全高效开发提供技术支持。
附图说明
图1为本发明的流程图;
图2为三维地质模型的实例图;
图3为三维孔隙压力场在三维油藏模型中的实例图;
图4为井眼轨迹附近网格加密实例图;
图5为初始应力场最小有效应力实例图;
图6为耦合迭代计算结果储层段最小有效应力实例图;
图7为本发明计算结果中最小水平主应力随孔隙压力变化情况;
图8为本发明计算结果中储层垂向沉降随孔隙压力变化情况;
图9为本发明计算结果中体积应变率随孔隙压力变化情况;
图10为本发明计算结果中孔隙比随孔隙压力变化情况;
图11为本发明计算结果中孔隙比随最小水平主应力变化情况;
图12为本发明计算结果中体积应变率随最小水平主应力变化情况。
具体实施方式
下面结合附图对本发明做进一步的描述,本发明的保护范围不局限于以下所述:
一种四维动态地应力模拟方法,它包括以下步骤:
S1、建立三维地质模型,该三维地质模型中至少应包含三维离散的地层厚度、储层物性参数(孔隙度、渗透率、饱和度、沉积相)及岩石力学参数(杨氏模量、泊松比)。
具体建立三维地质模型的步骤为:先根据地震资料建立三维地质模型几何形状,并根据其解释结果中的地震层位、断层、地震相、岩石类型、岩石属性等完善三维地质模型,然后结合测井资料、岩心资料等单井或单点资料分析储层物性参数和岩石力学参数并在横向上进行插值,最后生成包含岩石物性及岩石力学属性参数的三维地质模型,三维地质模型如图2所示。
S2、根据三维地质模型建立三维油藏模型,并利用生产或注入动态参数预测不同时期的三维孔隙压力场和温度场,同时建立非储层段模型并计算其属性参数,其具体包括以下三个步骤:
S2(I)将带有储层物性及岩石力学参数的三维地质模型导入到油藏模拟器中,建立有限差分网格的三维油藏模型;
S2(II)如图3所示,结合不同位置处单井在一定时间的生产/注入数据在三维油藏模型中进行渗流分析并预测不同时间的三维孔隙压力场和温度场;
S2(III)根据地震资料、钻井及录井资料对非储层段地层进行地层岩石物性及岩石力学参数插值计算。
S3、根据三维地质模型,建立带有储层物性及岩石力学属性的三维地应力模型。其具体建立过程为:将步骤S2中建立的三维油藏网格模型转换为三维地应力模型,并将油藏模型网格中的地层物性及岩石力学参数赋值到三维地应力模型对应网格中,如图4所示,以此建立带有储层物性及岩石力学属性的三维地应力模型。
S4、利用应力平衡法及单井地应力计算结果横向插值校正,形成初始三维地应力场,其具体包括以下两个步骤:
S4(I)在三维地应力模型中进行应力平衡计算,并将计算结果进行导入,其具体方法为:对模型施加重力载荷、孔隙压力等,根据地应力平衡得到有效垂向地应力,再选取适用的地应力模型,计算形成初始三维地应力场;
S4(II)利用初始应力提取法进行初始三维地应力场校正:根据测井资料、地应力测试结果等资料分析多口井的单井地应力纵向剖面,然后通过在横向上进行插值对初始三维地应力场进行校正,校正结果如图5所示。
S5、以油藏数值模拟结果导出的三维孔隙压力场和温度场作为不同计算时间步的边界条件,以应力初始化后的三维地应力模型为初始模型,建立四维动态地应力模型。其具体建立过程为:在不同时间步下对三维地应力模型施加载荷,载荷包括重力、静水压力、孔隙压力、温度等,并将步骤S2中预测得到的不同时间三维孔隙压力场和温度场作为外载荷,分别作为每一计算时间步的初始条件和边界条件,并以应力初始化后的三维地应力模型为初始模型,建立四维动态地应力模型。
S6、对四维动态地应力模型进行渗流-应力耦合迭代计算,得到动态地应力及孔弹性参数计算结果,计算结果如图6所示,通过计算结果分别分析地应力、地层位移、体积应变率、孔隙比、渗透率以及孔隙压力动态孔弹性参数的变化情况,这些变化情况以用于指导钻完井工程实践。
如图7~ 12所示,分别为对比分析不同储层岩石力学参数下地应力、地层位移、体积应变率、孔隙比、渗透率以及孔隙压力动态孔弹性参数的变化情况用于证明岩石力学特征描述在动态地应力分析中的重要性;如图7所示,地应力随孔隙压力变化情况为井筒稳定性和水力压裂提供基础参数;如图8 所示,储层垂向位移随孔隙压力的变化情况用于分析地层沉降和断层滑移;如图9所示,体积应变率随孔隙压力的变化情况用于分析地层沉降和断层滑移;如图10所示,孔隙比随孔隙压力的变化情况用于分析储层压实;如图11孔隙比与地应力的关系用于分析储层应力敏感性;如图12所示,应变率与地应力的关系用于分析储层应力敏感性。
因此该模拟方法能够模拟分析及预测油气资源在不同生产或注入时期内不同位置处地应力及储层和岩石力学参数的变化情况,能够准确描述油气储层孔弹性参数的非均质性和各向异性,更大程度上地反映储层的真实状态,进一步的克服了三维静态地应力模型无法准确反应油气开发过程中地应力及储层参数变化的问题。
以上所述仅是本发明的模型建立实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (1)
1.一种四维动态地应力模拟方法,其特征在于:它包括以下步骤:
S1、建立三维地质模型,该三维地质模型中至少应包含三维离散的地层厚度、储层物性参数及岩石力学参数;
S2、根据三维地质模型建立三维油藏模型,并利用生产或注入动态参数预测不同时期的三维孔隙压力场和温度场,同时建立非储层段模型并计算其属性参数;
S3、根据三维地质模型,建立带有储层物性及岩石力学属性的三维地应力模型;
S4、利用应力平衡法及单井地应力计算结果横向插值校正,形成初始三维地应力场;
S5、以油藏数值模拟结果导出的三维孔隙压力场和温度场作为不同计算时间步的边界条件,以应力初始化后的三维地应力模型为初始模型,建立四维动态地应力模型;
S6、对四维动态地应力模型进行渗流-应力耦合迭代计算,得到动态地应力及孔弹性参数计算结果,分析地应力、地层位移、体积应变率、孔隙比、渗透率以及孔隙压力动态孔弹性参数的变化情况,以用于指导钻完井工程实践。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217338.6A CN107066718B (zh) | 2017-04-05 | 2017-04-05 | 一种四维动态地应力模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217338.6A CN107066718B (zh) | 2017-04-05 | 2017-04-05 | 一种四维动态地应力模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107066718A true CN107066718A (zh) | 2017-08-18 |
CN107066718B CN107066718B (zh) | 2020-05-15 |
Family
ID=59601434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710217338.6A Active CN107066718B (zh) | 2017-04-05 | 2017-04-05 | 一种四维动态地应力模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107066718B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108227032A (zh) * | 2018-01-04 | 2018-06-29 | 中国石油天然气股份有限公司 | 一种地应力的确定方法及装置 |
CN108629463A (zh) * | 2018-05-23 | 2018-10-09 | 中国石油大学(北京) | 地应力预测方法和装置 |
CN109538185A (zh) * | 2018-10-30 | 2019-03-29 | 中国海洋石油集团有限公司 | 一种耦合温度场条件下的多层套管井井筒完整性分析模型 |
CN110454127A (zh) * | 2019-07-04 | 2019-11-15 | 成都理工大学 | 一种非常规油气藏加密井有利加密时间窗口确定方法 |
CN110705000A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 非常规储层加密井压裂动态微地震事件屏障区确定方法 |
CN112302601A (zh) * | 2019-07-25 | 2021-02-02 | 中国石油天然气集团有限公司 | 一种断层激活控制方法和装置 |
CN112906178A (zh) * | 2020-06-22 | 2021-06-04 | 中国海洋石油集团有限公司 | 一种利用点源函数计算水平注水井压力场方法 |
US11391854B2 (en) | 2019-07-04 | 2022-07-19 | Chengdu University Of Technology | Optimization design method for volumetric fracturing construction parameters of infilled well of unconventional oil and gas reservoir |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN104632157A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 低渗透油藏均衡驱替方法 |
CN105089582A (zh) * | 2015-05-28 | 2015-11-25 | 中国石油天然气股份有限公司 | 基于井下流量控制设备的油藏数值模拟方法及装置 |
CN105134189A (zh) * | 2015-08-24 | 2015-12-09 | 西南石油大学 | 一种LogGMIR测井地应力识别优质储层的方法 |
CN105201484A (zh) * | 2015-10-29 | 2015-12-30 | 西南石油大学 | 一种直井分层压裂层段优选及施工参数优化设计方法 |
-
2017
- 2017-04-05 CN CN201710217338.6A patent/CN107066718B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN104632157A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 低渗透油藏均衡驱替方法 |
CN105089582A (zh) * | 2015-05-28 | 2015-11-25 | 中国石油天然气股份有限公司 | 基于井下流量控制设备的油藏数值模拟方法及装置 |
CN105134189A (zh) * | 2015-08-24 | 2015-12-09 | 西南石油大学 | 一种LogGMIR测井地应力识别优质储层的方法 |
CN105201484A (zh) * | 2015-10-29 | 2015-12-30 | 西南石油大学 | 一种直井分层压裂层段优选及施工参数优化设计方法 |
Non-Patent Citations (3)
Title |
---|
WANG, LIANG.ETC: ""Mineral and pore structure characteristics of gas shale in Longmaxi formation: a case study of Jiaoshiba gas field in the southern Sichuan Basin, China"", 《ARABIAN JOURNAL OF GEOSCIENCES》 * |
孙金等: ""注水开发油藏温度对地应力的影响研究"", 《中国海上油气》 * |
移峥峰: ""页岩气应力/解吸/滑脱联合作用规律和多级渗流模型"", 《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅰ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108227032A (zh) * | 2018-01-04 | 2018-06-29 | 中国石油天然气股份有限公司 | 一种地应力的确定方法及装置 |
CN108227032B (zh) * | 2018-01-04 | 2020-03-10 | 中国石油天然气股份有限公司 | 一种地应力的确定方法及装置 |
CN108629463A (zh) * | 2018-05-23 | 2018-10-09 | 中国石油大学(北京) | 地应力预测方法和装置 |
CN108629463B (zh) * | 2018-05-23 | 2022-08-09 | 中国石油大学(北京) | 地应力变化预测方法和装置 |
CN109538185A (zh) * | 2018-10-30 | 2019-03-29 | 中国海洋石油集团有限公司 | 一种耦合温度场条件下的多层套管井井筒完整性分析模型 |
CN110454127A (zh) * | 2019-07-04 | 2019-11-15 | 成都理工大学 | 一种非常规油气藏加密井有利加密时间窗口确定方法 |
CN110705000A (zh) * | 2019-07-04 | 2020-01-17 | 成都理工大学 | 非常规储层加密井压裂动态微地震事件屏障区确定方法 |
US11391854B2 (en) | 2019-07-04 | 2022-07-19 | Chengdu University Of Technology | Optimization design method for volumetric fracturing construction parameters of infilled well of unconventional oil and gas reservoir |
CN110705000B (zh) * | 2019-07-04 | 2022-09-09 | 成都理工大学 | 非常规储层加密井压裂动态微地震事件屏障区确定方法 |
CN112302601A (zh) * | 2019-07-25 | 2021-02-02 | 中国石油天然气集团有限公司 | 一种断层激活控制方法和装置 |
CN112906178A (zh) * | 2020-06-22 | 2021-06-04 | 中国海洋石油集团有限公司 | 一种利用点源函数计算水平注水井压力场方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107066718B (zh) | 2020-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107066718A (zh) | 一种四维动态地应力模拟方法 | |
CN106991236B (zh) | 一种基于四维地应力动态变化的重复压裂选井选层方法 | |
CN110704888B (zh) | 一种非常规油气藏加密井体积压裂施工参数优化设计方法 | |
CN110472276B (zh) | 一种裂缝型油气储层加密井改造优化方法 | |
CN104992468A (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
Allis et al. | Update on subsidence at the Wairakei–Tauhara geothermal system, New Zealand | |
CN110705000B (zh) | 非常规储层加密井压裂动态微地震事件屏障区确定方法 | |
WO2010047859A1 (en) | Method for modeling deformation in subsurface strata | |
CN106014399B (zh) | 一种非均质地层高精度三维地应力模型建立方法 | |
CN106368694B (zh) | 一种针对复杂领域储层孔隙演化恢复与物性预测方法 | |
CN113011048B (zh) | 一种致密砾岩油藏水平井重复压裂模拟方法 | |
CN113591338B (zh) | 一种荷载和地下水开采引发地面沉降三维变参数全耦合模拟计算方法 | |
Xu et al. | Influence of tectonic uplift-erosion on formation pressure | |
CN113356843B (zh) | 针对地层的井壁稳定性分析方法、装置、介质及设备 | |
Cai et al. | Investigation of geomechanical response of fault in carbonate reservoir and its application to well placement optimization in YM2 Oilfield in Tarim Basin | |
Luo et al. | Discrete element-based regression analysis of initial ground stress and application to an extra-long tunnel in China | |
CN106154326A (zh) | 一种纵弯褶皱裂缝密度评价的方法及装置 | |
Zhang et al. | Present-day in situ stress prediction in Bozi 3 deep sandstone reservoir, Kuqa Depression: implications for gas development | |
Han et al. | In-situ and induced stresses in the development of unconventional resources | |
Diomampo et al. | Integrated dynamic reservoir modeling for multilayered tight gas sand development | |
Fan et al. | Stability Assessment and Optimization Design of Lakeside Open‐Pit Slope considering Fluid‐Solid Coupling Effect | |
Liu et al. | Numerical modeling on anisotropy of seepage and stress fields of stratified rock slope | |
Feng et al. | Wellbore Stability of a Deep‐Water Shallow Hydrate Reservoir Based on Strain Softening Characteristics | |
Jiang et al. | Establishment of geology-engineering integration model for shale gas of Wei AH platform in Weiyuan Area | |
Lamberghini et al. | Performance Based Reservoir Characterization in a Tight Gas Reservoir System-Case Study from Lajas and Punta Rosada Formations in the Neuquen Basin, Argentina |
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 |