CN103628866B - 获取地层的地层压力系数的方法 - Google Patents

获取地层的地层压力系数的方法 Download PDF

Info

Publication number
CN103628866B
CN103628866B CN201310597853.3A CN201310597853A CN103628866B CN 103628866 B CN103628866 B CN 103628866B CN 201310597853 A CN201310597853 A CN 201310597853A CN 103628866 B CN103628866 B CN 103628866B
Authority
CN
China
Prior art keywords
pressure coefficient
stratum
reservoir pressure
test point
coefficient
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
Application number
CN201310597853.3A
Other languages
English (en)
Other versions
CN103628866A (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 National Petroleum Corp
BGP Inc
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201310597853.3A priority Critical patent/CN103628866B/zh
Publication of CN103628866A publication Critical patent/CN103628866A/zh
Application granted granted Critical
Publication of CN103628866B publication Critical patent/CN103628866B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种获取地层的地层压力系数的方法，该方法包括以下步骤:(a)获取作为检测点的多个钻井的测井数据；(b)根据在步骤(a)中获取的测井数据，通过两种不同的计算方法计算得到各个钻井的所述地层的第一地层压力系数和所述地层的第二地层压力系数；(c)计算得出第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系；(d)获取工区平面内各检测点的地震数据，计算得到各检测点的第一地层压力系数和所述地层的第二地层压力系数，根据第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系计算得到各检测点的地层压力系数。本发明可以为预测页岩气产能提供技术支持。

Description

获取地层的地层压力系数的方法
技术领域
本发明涉及地层压力技术领域，更具体地讲，涉及一种获取地层的地层压力系数的方法。
背景技术
目前，我国页岩气勘探开发还处于起步阶段，具有风险大、成本高的特点。页岩气通常具有低渗透率和低孔隙度的特征，除少数裂缝发育带可能具有较高的自然产能外，其他地带的页岩气产能需要进行采集数据进行估算。
目前页岩气地震预测识别尚处于摸索阶段，常用的方法是利用常规的叠后速度反演和伽玛反演对优质页岩厚度进行预测，而优质页岩厚度并不能反应页岩地层气储量，对页岩开发有利区的选取的准确性不高，因此，还需要获取其他参数进行综合分析预测。
发明内容
本发明的目的在于提供一种获取地层的地层压力系数的方法，以解决目前页岩气产能预测缺乏技术支持的问题。
本发明的一方面提供一种获取地层的地层压力系数的方法，包括以下步骤:(a)获取作为检测点的多个钻井的测井数据，所述测井数据包括钻井处的所述地层的地震数据及所述地层的实测地层压力系数；(b)根据在步骤(a)中获取的测井数据，通过两种不同的计算方法计算得到各个钻井的所述地层的第一地层压力系数和所述地层的第二地层压力系数；(c)根据所述第一地层压力系数和所述第二地层压力系数与实测地层压力系数相关性大小，计算得出第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系；(d)获取工区平面内各检测点的所述地层的地震数据，计算得到各检测点的第一地层压力系数和第二地层压力系数，再根据步骤(c)中计算得到的第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系计算得到各检测点的地层压力系数。
优选地，所述地震数据，包括孔隙度接近于零的地理位置的地层速度、刚性接近于零的地理位置的地层速度、检测点的地层速度、平均地层压力系数、检测点的地层埋深。
优选地，第一地层压力系数通过以下公式计算得到：
P f = L n ( V i / V m a x ) L n ( V min / V m a x ) P o v ,
其中，Vmax为孔隙度接近于零的地理位置的地层速度，Vmin为刚性接近于零的地理位置的地层速度，Vi为检测点的地层速度，Pov为平均地层压力系数，Pf为第一地层压力系数。
优选地，第二地层压力系数通过已获取的地层埋深以及地层压力系数与地层埋深之间的函数关系计算获得。
优选地，所述函数关系通过已获取的多个钻井的测井数据：实测地层压力系数、检测点的地层埋深求取。
优选地，还包括步骤(e)，根据步骤(d)中计算得到的工区平面内各检测点的地层压力系数，绘制工区平面内的所述地层的地层压力系数分布图。
本发明获取地层的地层压力系数的方法，适用于获取地壳中的各个地层的地层压力系数。当本发明应用于页岩层时，能够获取工区内页岩层中的各检测点的地层压力系数，由于页岩地层压力系数为判断气藏量的重要的参数之一，因此可以为页岩气勘探开发有利区选取和页岩气井井位部署提供技术支持。
附图说明
通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中：
图1为根据本发明的实施例的获取地层的地层压力系数的方法的流程图。
图2为根据本发明实施例的地层压力系数分布图。
具体实施方式
现在，将参照附图更充分地描述不同的示例实施例，其中，一些示例性实施例在附图中示出。
图1示出根据本发明的实施例的获取地层的地层压力系数的方法的流程图。
如图1所示，根据本发明的实施例的获取地层的地层压力系数的方法包括以下步骤。
在步骤101，获取作为检测点的多个钻井的测井数据，所述测井数据包括钻井处的所述地层的地震数据和所述地层的实测地层压力系数P。所述多个钻井可以是区内已有的钻井，当然也可以是获取测井数据而专门开发的钻井。所述钻井的数量可以根据区内已有钻井的数量来选定。
所述地震数据可以包括孔隙度接近于零的地理位置的地层速度Vmax、刚性接近于零的地理位置的地层速度Vmin、检测点的地层速度Vi、平均地层压力系数Pov、检测点的地层埋深Zi。所述检测点为需要计算的地层压力系数的地层所在的地理位置。所述的多个钻井的测井数据可以通过原始的测井资料计算获取，其计算获取方法为现有技术。例如检测点的地层速度Vi的获取方法为：利用工区的声波测井数据和井旁地震道的地震子波来合成地震记录，再用合成记录对井旁地震道进行层位标定，建立每一口井的时深关系，即为初始阻抗模型，以此为基础，进行基于模型的速度反演运算，得到检测点的地层速度Vi。
在步骤102，根据在步骤101获取的测井数据，通过两种不同的计算方法计算得到各个钻井的所述地层的第一地层压力系数Pf和所述地层的第二地层压力系数P2。
第一地层压力系数Pf可以通过以下公式计算得到：
P f = L n ( V i / V m a x ) L n ( V min / V m a x ) P o v .
即将步骤101中获取的各钻井处的所述地层的地震数据代入到公式中，通过计算得到各钻井的所述地层的第一地层压力系数Pf。
第二地层压力系数P2可以通过已获取地震数据中的检测点的地层埋深Zi以及地层压力系数与地层埋深之间的函数关系计算获得。地层压力系数与地层埋深之间的函数关系可以是本领域的技术人员通过已有经验获取；优选通过已获取的多个钻井的测井数据：实测地层压力系数P、检测点的地层埋深Zi求取，求取方法可以是拟合。
在步骤103，根据各个钻井的所述第一地层压力系数Pf和所述第二地层压力系数P2与实测地层压力系数P相关性大小，计算得出第一地层压力系数Pf和第二地层压力系数P2与实测地层压力系数P的权重关系。在实际应用中，可以将各个钻井的所述第一地层压力系数Pf和所述第二地层压力系数P2与实测地层压力系数P分别输入到计算机计算软件中，计算机计算软件将输出上述得出第一地层压力系数Pf和第二地层压力系数P2与实测地层压力系数P的权重关系。
在步骤104，获取工区平面内各检测点的所述地层的地震数据，计算得到各检测点的第一地层压力系数和第二地层压力系数，再根据步骤103中计算得到的第一地层压力系数Pf和第二地层压力系数P2与实测地层压力系数P的权重关系计算得到各检测点的地层压力系数。
所述的工区平面内各检测点的所述地层的地震数据包括孔隙度接近于零的地层速度Vmax、刚性接近于零的地理位置的地层速度Vmin、检测点的地层速度Vi、平均地层压力系数Pov、检测点的地层埋深Z，这些数据可以通过以本领域公知的方法进行获取，例如，可以通过野外勘探采集得到原始的数据，然后，经过地震资料处理，得到纯波带地震偏移数据体，再开展构造精细解释和地震速度反演工作，获得检测点的地层速度Vi和检测点的地层埋深Zi。
根据获取到的地震数据计算各检测点的第一地层压力系数Pf和第二地层压力系数P2的计算方法与步骤102中的第一地层压力系数Pf和第二地层压力系数P2计算方法相同。
再根据计算得到的各检测点的第一地层压力系数Pf和第二地层压力系数P2及步骤103中计算得到的第一地层压力系数Pf和第二地层压力系数P2与实测地层压力系数P的权重关系计算得到各检测点的地层压力系数。
根据需要，可以利用计算获得的工区内预定地层的各位置点的地层压力系数，编制工区平面内的预定地层的地层压力系数分布图。图2示出了根据本发明实施例的页岩气地层的地层压力系数分布图，由于页岩地层压力系数为判断气藏量的重要的参数之一，因此根据工区内地层压力系数分布图可以为页岩气勘探开发有利区选取和页岩气井井位部署提供技术支持。
在编制工区平面内的预定地层的地层压力系数分布图时，步骤104中的各检测点可以包括步骤101中的多个钻井的位置，即将多个钻井处的预定地层的地层压力系数按照图1中的方法重新进行计算获取；步骤104中的各检测点也可以不包括步骤101中的多个钻井的位置，即将多个钻井处的预定地层的实测地层压力系数作为多个钻井位置处的预定地层的地层压力系数。
在应用地层压力系数分布图的过程中，应该考虑断层影响，在大断层附近，不适宜进行页岩气钻井部署。可以通过构造精细解释和根据断层断距的大小，划分区内断层级别。
可以理解，根据本发明实施例获取地层的地层压力系数的方法中的地层不限于页岩气地层，该方法可以应用于获取地壳中的其他地层的地层压力系数。
除非另有定义，否则这里使用的所有术语(包括技术和科学术语)具有与示例性实施例所属领域的普通技术人员通常理解的含义相同的含义。还应该理解，除非这里明确定义，否则术语(诸如在常用词典中定义的)应被解释为具有与所述术语在相关领域的上下文中的含义一致的含义，而不应理想化或过于形式地被理解。
尽管已经参照其示例性实施例具体显示和描述了本发明，但是本领域的技术人员应该理解，在不脱离权利要求所限定的本发明的精神和范围的情况下，可以对其进行形式和细节上的各种改变。

Claims (5)

1.一种获取地层的地层压力系数的方法，其特征在于，包括以下步骤:
(a)获取作为检测点的多个钻井的测井数据，所述测井数据包括钻井处的所述地层的地震数据及所述地层的实测地层压力系数；
(b)根据在步骤(a)中获取的测井数据，通过两种不同的计算方法计算得到各个钻井的所述地层的第一地层压力系数和所述地层的第二地层压力系数；
(c)根据所述第一地层压力系数和所述第二地层压力系数与实测地层压力系数相关性大小，计算得出第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系；
(d)获取工区平面内各检测点的所述地层的地震数据，计算得到各检测点的第一地层压力系数和第二地层压力系数，再根据步骤(c)中计算得到的第一地层压力系数和第二地层压力系数与实测地层压力系数的权重关系计算得到各检测点的地层压力系数，其中，工区平面内的检测点包括步骤(a)中的作为检测点的多个钻井的位置，或者不包括步骤(a)中的作为检测点的多个钻井的位置，
其中，还包括步骤(e)，根据步骤(d)中计算得到的工区平面内各检测点的地层压力系数，绘制工区平面内的所述地层的地层压力系数分布图。
2.根据权利要求1所述的获取地层的地层压力系数的方法，其特征在于，所述地震数据，包括孔隙度接近于零的地理位置的地层速度、刚性接近于零的地理位置的地层速度、检测点的地层速度、平均地层压力系数、检测点的地层埋深。
3.根据权利要求2所述的获取地层的地层压力系数的方法，其特征在于，第一地层压力系数通过以下公式计算得到：
P f = L n ( V i / V max ) L n ( V min / V max ) P o v ,
其中，Vmax为孔隙度接近于零的地理位置的地层速度，Vmin为刚性接近于零的地理位置的地层速度，Vi为检测点的地层速度，Pov为平均地层压力系数，Pf为第一地层压力系数。
4.根据权利要求2所述的获取地层的地层压力系数的方法，其特征在于，第二地层压力系数通过已获取的地层埋深以及地层压力系数与地层埋深之间的函数关系计算获得。
5.根据权利要求4所述的获取地层的地层压力系数的方法，其特征在于，所述函数关系通过已获取的多个钻井的测井数据：实测地层压力系数、检测点的地层埋深求取。
CN201310597853.3A 2013-11-22 2013-11-22 获取地层的地层压力系数的方法 Active CN103628866B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310597853.3A CN103628866B (zh) 2013-11-22 2013-11-22 获取地层的地层压力系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310597853.3A CN103628866B (zh) 2013-11-22 2013-11-22 获取地层的地层压力系数的方法

Publications (2)

Publication Number Publication Date
CN103628866A CN103628866A (zh) 2014-03-12
CN103628866B true CN103628866B (zh) 2017-02-08

Family

ID=50210281

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310597853.3A Active CN103628866B (zh) 2013-11-22 2013-11-22 获取地层的地层压力系数的方法

Country Status (1)

Country Link
CN (1) CN103628866B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104484573B (zh) 2014-12-30 2017-09-15 中国石油天然气股份有限公司 一种确定地层刚性系数的方法
CN104483703B (zh) * 2014-12-31 2017-05-10 中国石油天然气股份有限公司 一种井震联合古地层压力预测方法
CN105298483B (zh) * 2015-10-22 2018-03-09 中国石油天然气股份有限公司 一种获得注水过程中储层综合伤害的方法及装置
CN109339771B (zh) * 2017-08-01 2022-04-12 中国石油化工股份有限公司 一种页岩油气层孔隙压力预测方法及系统
CN114135276A (zh) * 2020-08-12 2022-03-04 中国石油天然气集团有限公司 地层压力的确定方法、确定装置、处理器和确定系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101906966A (zh) * 2010-07-16 2010-12-08 中国石油天然气股份有限公司 储层产能的预测方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103382838A (zh) * 2013-07-25 2013-11-06 中国石油大学(北京) 一种基于压裂地质体可压性的储层分析方法及装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101906966A (zh) * 2010-07-16 2010-12-08 中国石油天然气股份有限公司 储层产能的预测方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Fillippone地层压力预测方法的改进及应用;马海;《石油钻探技术》;20121130;第40卷(第6期);第56-61页 *
三维地层压力预测方法及应用研究;陈学国;《石油钻探技术》;20050531;第33卷(第1期);第13-15页 *
地震资料预测地层压力方法综述;王越之;《钻采工艺》;19971231;第20卷(第2期);第5-8页 *
用地震层速度预测地层压力的方法;郭燕珩;《内蒙古石油化工》;20081231(第5期);第66-67页 *
辽西凹陷北洼下第三系异常地层压力分析;刘震等;《石油学报》;19930131;第14卷(第1期);第14-23页 *

Also Published As

Publication number Publication date
CN103628866A (zh) 2014-03-12

Similar Documents

Publication Publication Date Title
EP2960680B1 (en) Fracturing and reactivated fracture volumes
EP3036398B1 (en) Formation stability modeling
US8768672B2 (en) Method for predicting time-lapse seismic timeshifts by computer simulation
EP2649271B1 (en) A method to improve reservoir simulation and recovery from fractured reservoirs
US8548782B2 (en) Method for modeling deformation in subsurface strata
US10007015B2 (en) Methods, systems and devices for predicting reservoir properties
CA2913827C (en) Methods, systems and devices for predicting reservoir properties
EP3018502A2 (en) Modeling fluid-conducting fractures in reservoir simulation grids
CN103628866B (zh) 获取地层的地层压力系数的方法
CN101421640A (zh) 用于预测井位的碳氢化合物产量的方法和设备
CN103914620B (zh) 一种计算断层破裂带裂缝张开空间分布的方法
CN105588883B (zh) 三维岩石力学参数获取方法和系统
CN103038670A (zh) 对地球物理数据进行处理
US11578564B2 (en) Systems and methods for predicting shear failure of a rock formation
CN106443770A (zh) 一种页岩气地质甜点的预测方法
Chen et al. Study on the application of a comprehensive technique for geological prediction in tunneling
CN104502996A (zh) 一种密度曲线校正方法及系统
US20160377752A1 (en) Method of Digitally Identifying Structural Traps
CN104820237B (zh) 计算地层的方差体的方法
CN106842316A (zh) 裂缝确定方法和装置
CN104749618B (zh) 泥页岩缓倾角裂缝叠后概率定量表征方法
MX2015004890A (es) Aparatos, metodos y sistemas de doble busqueda de datos.
CN106650036B (zh) 水平井钻井方法和装置
WO2010047858A1 (en) Method for predicting time-lapse seismic timeshifts by computer simulation
RU2563862C2 (ru) Атрибут напряжения в горных породах

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180212

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200916

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Co-patentee after: BGP Inc., China National Petroleum Corp.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.