CN111722285A - 一种基于微地震数据压后页岩气藏建模方法 - Google Patents

一种基于微地震数据压后页岩气藏建模方法 Download PDF

Info

Publication number
CN111722285A
CN111722285A CN201910220274.4A CN201910220274A CN111722285A CN 111722285 A CN111722285 A CN 111722285A CN 201910220274 A CN201910220274 A CN 201910220274A CN 111722285 A CN111722285 A CN 111722285A
Authority
CN
China
Prior art keywords
shale gas
gas reservoir
modeling
data
grid
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
Application number
CN201910220274.4A
Other languages
English (en)
Other versions
CN111722285B (zh
Inventor
戴城
计秉玉
王卫红
刘华
胡小虎
郭艳东
王妍妍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910220274.4A priority Critical patent/CN111722285B/zh
Publication of CN111722285A publication Critical patent/CN111722285A/zh
Application granted granted Critical
Publication of CN111722285B publication Critical patent/CN111722285B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出了一种基于微地震数据压后页岩气藏建模方法,根据水力裂缝和天然裂缝网络的几何特点,分别采用霍夫变换和离散裂缝粗化技术重构多尺度裂缝系统,建立压后气藏模型。根据本发明的方法建立的模型充分利用了微地震分布信息,重构的裂缝网络更加准确,在此基础上可以更加精确的进行页岩气藏数值模拟,从而更加准确的进行页岩气井产能预测,对开发方案进行设计、调整。

Description

一种基于微地震数据压后页岩气藏建模方法
技术领域
本发明涉及石油开采技术领域,尤其涉及一种基于微地震数据建立压后页岩气藏地质模型的方法。
背景技术
由于页岩气藏储层致密低渗的特点,为提高页岩气产能,必须在页岩气开发过程中采用分段压裂水平井技术对页岩储层进行改造。分段压裂技术形成不同尺度的人工水力裂缝沟通天然裂缝网络会显著提高页岩储层的渗透率,对页岩储层内流体流动影响很大,将直接影响气藏开发效果。数值模拟技术是预测气井产能的重要工具。而建立准确的压后页岩气藏模型是利用数值模拟预测产能的基础。
微地震监测是页岩气藏压裂改造领域重要新兴技术,它通过接收、解释水力压裂过程中释放的微地震信号为评估水力压裂改造效果提供依据。现有根据微地震数据建立压后气藏的方法主要可以分为两类:一是根据微地震事件发生的位置及强度信息建立等效双重介质模型;二是连接微地震事件点,建立离散裂缝模型(discrete fracture model,DFM)。以上两种方法都存在一定的不完善之处。第一种方法建立的等效双重介质模型均化了压裂改造区裂缝信息,不能考虑水力裂缝的真实几何信息。水力裂缝是渗流发生的主要通道,对流动影响巨大。忽略水力裂缝的几何信息会一定程度上造成模拟结果的不精确。微地震信号通常只反映水力裂缝开裂导致孔隙压力扩散诱发周围天然裂缝的剪切形变。现有第二类方法直接将微地震事件点相连不足建立离散裂缝模型,不足以反应水力裂缝的真实几何形态。
因此,需要一种新的建模方法以准确构建压后气藏模型。
发明内容
为了解决上述技术问题,本发明提供了一种在压前页岩气地质模型基础上,根据微地震数据,建立压后页岩气藏的方法。该方法能够根据水力裂缝和天然裂缝网络的几何特点,分别采用霍夫变换和离散裂缝粗化技术重构多尺度裂缝系统,建立压后气藏模型,为后续数值模拟奠定基础。
本发明的基于微地震数据压后页岩气藏建模方法,其特征在于,包括以下步骤:
步骤1.数据准备和处理;
步骤2.利用霍夫变换重构水力裂缝;
步骤3.根据微地震的强度,计算裂缝的开度和渗透率;
步骤4.将重构的水力裂缝嵌入至压前地质模型中,建立压后页岩气藏模型。
进一步的,在所述步骤1中,所述数据包括压前地震模型和微地震数据。
进一步的,所述步骤2包括:
步骤2.1.定义方位角
Figure BDA0002003383420000021
倾角θ和距离r张成的新空间为霍夫空间,微地震事件点(x,y,z)变换至霍夫空间,参考方程为
Figure BDA0002003383420000022
步骤2.2.将所有N个微地震事件位置为(x0,y0,z0),(x1,y1,z1)……(xN,yN,zN),强度为M1,M2,…,MN变换至霍夫空间组成方程组
Figure BDA0002003383420000023
步骤2.3.求解微地震点满足的曲面参数方程组成的超定方程组。
进一步的,在所述步骤3中,假设步骤三中某条由n个强度分别为Mn1,Mn2,…,Mnn的微地震事件重构得到的裂缝f,其开度为
Figure BDA0002003383420000024
其渗透率为
Figure BDA0002003383420000025
进一步的,所述步骤4包括:
步骤4.1.筛选所有与裂缝片有接触的基质网格;
步骤4.2.生成裂缝平面与基质网格的交面;
步骤4.3.根据裂缝外形剪裁交面;
步骤4.4.计算基质网格-交面的传导率;
步骤4.5.搜索裂缝内部交面的连接,计算交面之间的传导率;
步骤4.6.计算跨裂缝的传导率。
进一步的,在所述2.3中,求解所述超定方程租的方法采用投票函数方法。
进一步的,在所述步骤2.3中,将整个霍夫空间离散为NX×NY×NZ个网格组成的离散霍夫空间,遍历离散霍夫空间内每一个网格,统计步骤2.2中方程所代表曲面和网格相交的个数,记为该网格的投票数V,当V大于微地震事件总数的十分之一时,网格的中心点坐标
Figure BDA0002003383420000031
则为重构裂缝的方位角,倾角和离原点的距离。
进一步的,在所述步骤4.2中,所述交面为多边形,多边形的边数≥3。
与现有技术相比,本发明的基于微地震数据压后页岩气藏建模方法能够基于已经获得的压前地质模型和微地震数据建立压后地质模型,该模型充分利用了微地震分布信息,重构的裂缝网络更加准确。在此基础上可以更加精确的进行页岩气藏数值模拟,从而更加准确的进行页岩气井产能预测,对开发方案进行设计、调整。并且本发明可以直接采用Petrel地质建模软件构建的压前地质模型,无需对压前地质模型进行任何非结构化网格剖分处理,使用更为方便。裂缝基于微地震信息重构之后,直接嵌入至压前地质模型中,无需进行粗化,完全保留了裂缝的分布,走向,开度等细节,使得基于本发明构建的压后地质模型的数值模拟更加精确。
上述技术特征可以各种技术上可行的方式组合以产生新的实施方案,只要能够实现本发明的目的。
附图说明
在下文中将基于仅为非限定性的实施例并参考附图来对本发明进行更详细的描述。其中:
图1显示了根据本发明的基于微地震数据压后页岩气藏建模方法的流程图;
图2显示了本发明方法中的霍夫变换示意图;
图3显示了本发明方法中将重构裂缝嵌入压前地质模型流程示意图;
图4显示了本发明实施例中微地震分布示意图;
图5显示了本发明实施例中利用本方法建模后模型示意图。
具体实施方式
以下将结合说明书附图和具体实施例对本发明做进一步详细说明。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
如图1所示,本发明的基于微地震数据压后页岩气藏建模方法,包括以下步骤,
步骤1.数据准备和处理;
所需数据包括压前地质模型和微地震数据,其中,压前地质模型由两组数据定义,角点网格的柱线坐标COORD和网格的上、下表面与柱线交点的深度ZCORN。当压前地质模型的网格数为NX×NY×NZ时(NX、NY和NZ分别为x,y,z方向的网格数),COORD的数量为(NX+1)×(NY+1)×6,ZCORN数量为NX×NY×NZ×8。COORD和ZCORN数据可以通过地质建模软件Petrel导出。
微地震数据包含四列内容,前三列为微地震发生位置的空间坐标(x,y,z),第四列为微地震的强度。微地震数据可由微地震监测设备直接获得。
步骤2.利用霍夫变换重构水力裂缝,具体包括:
定义方位角
Figure BDA0002003383420000041
倾角θ和距离r张成的新空间为霍夫空间,如图2所示,微地震时间点(x,y,z)变换至霍夫空中变为曲面,参考方程为,
Figure BDA0002003383420000042
步骤2.2.将所有N个微地震事件位置为(x0,y0,z0),(x1,y1,z1)……(xN,yN,zN),强度为M1,M2,…,MN变换至霍夫空间组成方程组
Figure BDA0002003383420000043
步骤2.3.求解微地震点满足的曲面参数方程组成的超定方程组。
具体的,采用投票函数方法,将整个霍夫空间离散为NX×NY×NZ个网格组成的离散霍夫空间(NX、NY和NZ分别为x,y,z方向的网格数)。遍历离散霍夫空间内每一个网格,统计步骤2.2中方程所代表曲面和该网格相交的个数,记为该网格的投票数V,当V大于微地震事件总数的十分之一时,该网格的中心点坐标
Figure BDA0002003383420000044
则为重构裂缝的方位角,倾角和离原点的距离。
步骤3.根据微地震的强度,计算裂缝的开度和渗透率;
假设步骤三中某条由n个微地震事件(强度分别为Mn1,Mn2,…,Mnn)重构得到的裂缝f,其开度为
Figure BDA0002003383420000051
其渗透率为
Figure BDA0002003383420000052
步骤4.将重构的水力裂缝嵌入至压前地质模型中,建立压后页岩气藏模型。
包括:
步骤4.1.筛选所有与裂缝片有接触的基质网格;
步骤4.2.生成裂缝平面与基质网格的交面,交面有可能是三角形、四边形、五边形、六边形等多边形;
步骤4.3.根据裂缝外形剪裁交面;
步骤4.4.计算基质网格-交面的传导率;
步骤4.5.搜索裂缝内部多边形的连接,计算多边形之间的传导率;
步骤4.6.搜索裂缝与裂缝的连接,求出裂缝与裂缝的交线,然后在两个裂缝中分别找到交线穿过的多边形,计算共享交线的多边形之间的传导率。
本发明的方法建立的压后页岩气藏建立压后页岩气藏地质模型,一方面可以直接采用Petrel地质建模软件构建的压前地质模型,无需对压前地质模型进行任何非结构化网格剖分处理,使用更为方便,另一方面充分利用了微地震分布信息重构裂缝,重构的裂缝直接嵌入至压前地质模型中,无需进行粗化,完全保留了裂缝的分布,走向,开度等细节,这使得基于本发明构建的压后地质模型的数值模拟更加精确。
按照本发明提供的基于基于微地震数据压后页岩气藏建模方法,以F气田A平台为例,进行了现场实施。
该区块共有三口井,即A-1HF,A-2HF和A-3HF。三口井井长分别为1468m,1262m,1248m,井距为600m。该井平台共覆盖3.24平方公里面积。储层为五峰-龙马溪组页岩储层,埋深为2550m,储层温度为85摄氏度。储层共划分为9小层,1-5小层为开发储层,厚度为40m,束缚水饱和度40%。对A平台进行了微地震监测来探测裂缝的扩展。图4展示的是A平台微地震的分布情况。在平面上,微地震事件分布在井周围250m范围之内。A-2HF和A-3HF微地震事件存在重合区域。在纵向方向上,微地震事件分布均匀,在各个小层中均有分布。图5展示的是利用本发明的方法建立的A平台压后地质模型。
从图4和图5的比较可以看出,利用本发明的方法建立的A平台压后地质模型完全保留了裂缝的分布,走向,开度等细节,压后地质模型的数值模拟更加精确。
至此,虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。

Claims (8)

1.一种基于微地震数据压后页岩气藏建模方法,其特征在于,包括以下步骤:
步骤1.数据准备和处理;
步骤2.利用霍夫变换重构水力裂缝;
步骤3.根据微地震的强度,计算裂缝的开度和渗透率;
步骤4.将重构的水力裂缝嵌入至压前地质模型中,建立压后页岩气藏模型。
2.根据权利要求1所述的基于微地震数据压后页岩气藏建模方法,其特征在于,在所述步骤1中,所述数据包括压前地震模型和微地震数据。
3.根据权利要求2所述的基于微地震数据压后页岩气藏建模方法,其特征在于,所述步骤2包括:
步骤2.1.定义方位角
Figure FDA0002003383410000011
倾角θ和距离r张成的新空间为霍夫空间,微地震事件点(x,y,z)变换至霍夫空间变为曲面,参考方程为
Figure FDA0002003383410000012
步骤2.2.将所有N个微地震事件位置为(x0,y0,z0),(x1,y1,z1)……(xN,yN,zN),强度为M1,M2,…,MN变换至霍夫空间组成方程组
Figure FDA0002003383410000013
步骤2.3.求解微地震点满足的曲面参数方程组成的超定方程组。
4.根据权利要求3所述的基于微地震数据压后页岩气藏建模方法,其特征在于,在所述步骤3中,假设步骤三中某条由n个强度分别为Mn1,Mn2,…,Mnn的微地震事件重构得到的裂缝f,其开度为
Figure FDA0002003383410000014
其渗透率为
Figure FDA0002003383410000015
5.根据权利要求1所述的基于微地震数据压后页岩气藏建模方法,其特征在于,所述步骤4包括:
步骤4.1.筛选所有与裂缝片有接触的基质网格;
步骤4.2.生成裂缝平面与基质网格的交面;
步骤4.3.根据裂缝外形剪裁交面;
步骤4.4.计算基质网格-交面的传导率;
步骤4.5.计算各交面之间的传导率;
步骤4.6.计算跨裂缝的传导率。
6.根据权利要求3所述的基于微地震数据压后页岩气藏建模方法,其特征在于,在所述2.3中,求解所述超定方程租的方法采用投票函数方法。
7.根据权利要求6所述的基于微地震数据压后页岩气藏建模方法,其特征在于,在所述步骤2.3中,将整个霍夫空间离散为NX×NY×NZ个网格组成的离散霍夫空间,NX,NY和NZ分别为x,y,z方向的网格数,遍历离散霍夫空间内每一个网格,统计步骤2.2中方程所代表曲面和网格相交的个数,记为该网格的投票数V,当V大于微地震事件总数的十分之一时,网格的中心点坐标
Figure FDA0002003383410000021
为重构裂缝的方位角,倾角和离原点的距离。
8.根据权利要求5所述的基于微地震数据压后页岩气藏建模方法,其特征在于,在所述步骤4.2中,所述交面为多边形,多边形的边数≥3。
CN201910220274.4A 2019-03-22 2019-03-22 一种基于微地震数据压后页岩气藏建模方法 Active CN111722285B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910220274.4A CN111722285B (zh) 2019-03-22 2019-03-22 一种基于微地震数据压后页岩气藏建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910220274.4A CN111722285B (zh) 2019-03-22 2019-03-22 一种基于微地震数据压后页岩气藏建模方法

Publications (2)

Publication Number Publication Date
CN111722285A true CN111722285A (zh) 2020-09-29
CN111722285B CN111722285B (zh) 2022-05-03

Family

ID=72562187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910220274.4A Active CN111722285B (zh) 2019-03-22 2019-03-22 一种基于微地震数据压后页岩气藏建模方法

Country Status (1)

Country Link
CN (1) CN111722285B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112270093A (zh) * 2020-10-28 2021-01-26 成都北方石油勘探开发技术有限公司 一种基于微地震矩张量反演的储层天然裂缝建模方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010117575A2 (en) * 2009-04-07 2010-10-14 Virginia Commonwealth University Accurate pelvic fracture detection for x-ray and ct images
CN103698817A (zh) * 2013-12-27 2014-04-02 上海川海信息科技有限公司 基于图像的地下洞室掌子面快速安全预测方法
CN104459775A (zh) * 2014-11-28 2015-03-25 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于微地震监测数据的页岩气藏裂缝建模方法
WO2017206159A1 (en) * 2016-06-03 2017-12-07 Schlumberger Technology Corporation Fracture network extraction by microseismic events clustering analysis
CN108240214A (zh) * 2018-01-18 2018-07-03 西南石油大学 一种页岩气藏水平井分段压裂压后产能预测方法
CN108460744A (zh) * 2018-03-21 2018-08-28 武汉理工大学 一种水泥刻槽路面图像降噪增强与裂缝特征提取方法
CN109490949A (zh) * 2018-11-29 2019-03-19 四川圣诺油气工程技术服务有限公司 一种基于微地震解释结果的压后产能评价方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010117575A2 (en) * 2009-04-07 2010-10-14 Virginia Commonwealth University Accurate pelvic fracture detection for x-ray and ct images
CN103698817A (zh) * 2013-12-27 2014-04-02 上海川海信息科技有限公司 基于图像的地下洞室掌子面快速安全预测方法
CN104459775A (zh) * 2014-11-28 2015-03-25 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于微地震监测数据的页岩气藏裂缝建模方法
WO2017206159A1 (en) * 2016-06-03 2017-12-07 Schlumberger Technology Corporation Fracture network extraction by microseismic events clustering analysis
CN108240214A (zh) * 2018-01-18 2018-07-03 西南石油大学 一种页岩气藏水平井分段压裂压后产能预测方法
CN108460744A (zh) * 2018-03-21 2018-08-28 武汉理工大学 一种水泥刻槽路面图像降噪增强与裂缝特征提取方法
CN109490949A (zh) * 2018-11-29 2019-03-19 四川圣诺油气工程技术服务有限公司 一种基于微地震解释结果的压后产能评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHENG DAI ET AL.: "Analysis of the Influencing Factors on the Well Performance in Shale Gas Reservoir", 《GEOFLUIDS》 *
SONG DU ET AL.: "Field study:embedded discrete fracture modeling with artificial intelligence in permian basin for shale formaton", 《SOCIETY OF PETROLEUM ENGINEERS》 *
严侠等: "基于模拟有限差分的嵌入式离散裂缝数学模型", 《中国科学》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112270093A (zh) * 2020-10-28 2021-01-26 成都北方石油勘探开发技术有限公司 一种基于微地震矩张量反演的储层天然裂缝建模方法
CN112270093B (zh) * 2020-10-28 2022-11-08 成都北方石油勘探开发技术有限公司 一种基于微地震矩张量反演的储层天然裂缝建模方法

Also Published As

Publication number Publication date
CN111722285B (zh) 2022-05-03

Similar Documents

Publication Publication Date Title
CN104992468B (zh) 缝洞型碳酸盐岩油气藏三维地质建模方法
CN114663627B (zh) 一种基于三维点云数据库的矿山数字模型建立方法
CN110147561B (zh) 一种含天然裂缝致密油气储层体积压裂缝网预测方法
CA2885293C (en) Analyzing fracture stratigraphy
Qu et al. Incorporation of deformation band fault damage zones in reservoir models
CN110595887B (zh) 一种浅埋硬岩大型地下洞室成洞条件量化评价方法
CN105701319B (zh) 一种水平井规则开发井网下的沉积微相建模方法
CN102609982B (zh) 空间地质数据非结构化模式的拓扑发现方法
CN107742020A (zh) 页岩气储层压裂改造体积预测方法
EP2818897B1 (en) Method and apparatus for estimating a parameter of a subsurface volume
CN101726750A (zh) 地裂缝评价的“三图法”
CN112035939B (zh) 一种双侧壁导坑隧道的岩土体参数随机场建模方法
CN108241170B (zh) 碳酸盐岩储层缝洞雕刻体积与动态储量拟合方法及装置
CN115937467B (zh) 随机三维裂隙在升尺度模型网格中的划分方法及系统
CN105719340A (zh) 一种油气储层中人工压裂裂缝的三维建模方法
CN113177322A (zh) 一种压裂单井控制储量计算方法
CN111722285B (zh) 一种基于微地震数据压后页岩气藏建模方法
Du et al. Generalization of dual-porosity-system representation and reservoir simulation of hydraulic fracturing-stimulated shale gas reservoirs
CN106372295B (zh) 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法
CN110705168A (zh) 构造应力场的模拟方法
CN109558614B (zh) 页岩气藏多尺度裂缝内气体流动的模拟方法及系统
KR101131517B1 (ko) 물리검층 자료의 통계학적 해석결과를 활용한 오일샌드 지역 지하정보의 처리시스템 및 이를 이용한 오일샌드 지역 지하정보의 처리방법
Liu et al. Study of roof water inrush forecasting based on EM-FAHP two-factor model
CA2943189C (en) Bin constraints for generating a histogram of microseismic data
US11353622B2 (en) Systems and methods for hydrocarbon reservoir three dimensional unstructured grid generation and development

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