CN107817520A - 海相泥页岩地层的压力系数预测方法及系统 - Google Patents
海相泥页岩地层的压力系数预测方法及系统 Download PDFInfo
- Publication number
- CN107817520A CN107817520A CN201710860322.7A CN201710860322A CN107817520A CN 107817520 A CN107817520 A CN 107817520A CN 201710860322 A CN201710860322 A CN 201710860322A CN 107817520 A CN107817520 A CN 107817520A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- velocity
- formula
- optimization
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 208000035126 Facies Diseases 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 claims abstract description 64
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 44
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000000491 multivariate analysis Methods 0.000 claims description 5
- 238000004321 preservation Methods 0.000 abstract description 9
- 238000011160 research Methods 0.000 abstract description 9
- 235000021185 dessert Nutrition 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 5
- 238000005056 compaction Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开了一种海相泥页岩地层的压力系数预测方法及系统。该方法包括：准备钻井数据、测井数据与地震资料；针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；基于钻井数据，分析钻井密度与纵波速度的关系，获得Gardner公式，对Fillippone公式进行优化，获取优化Fillippone公式，拟合优化Fillippone公式中的系数；基于横波速度对优化Fillippone公式进行再次优化，获取再优化Fillippone公式，拟合再优化Fillippone公式中的系数；基于再优化Fillippone公式，预测地层压力系数。本发明通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区。
Description
技术领域
本发明涉及页岩气勘探领域，更具体地，涉及一种海相泥页岩地层的压力系数预测方法及系统。
背景技术
近年来页岩气勘探开发的实践表明良好的保存条件是页岩气高产富集的关键，而压力系数是保存条件的综合判别指标，页岩气产量与压力系数呈正相关关系，因而压力系数的预测及研究对页岩气勘探获得成功至关重要。目前利用地震信息进行压力预测主要有两种方法，大致可分为图解法和公式计算法，图解法包括等效深度图解法、比值法或差值法和量板法；公式计算法包括等效深度公式计算法、Eaton法、Fillippone法、刘震云法、Stone法、Martinez法及有效应力法等。有效应力法、Eaton法及Fillippone法是目前应用最广且技术相对成熟的方法，勘探实践表明利用上述方法的海相泥页岩地层压力系数的预测精度均较低，尤其是复杂构造区，认为主要存在如下几个问题：
(1)川东南地区志留系海相地层以泥岩沉积为主，速度自上而下降低，不符合正常压实理论，是在实际应用中经常难以准确地建立正常压实趋势线，因而等效深度法及Eaton法的预测精度均较低；
(2)有效应力法原理是通过弹性参数反推地层的有效应力，需要已知应变ΔH/H的大小，但应变的大小无法准确获得，受地层埋藏深度控制，尤其是构造起伏较大的构造复杂区，压力系数预测精度较低。
(3)Fillippone法及其改进方法基于的原理是不均衡压实及有机质生烃产生高孔隙压力，形成欠压实，地震波波速比正常压实的波速要小，实现相对简单，陈超及屈大鹏分别在文章《川东南焦石坝地区海相泥页岩含气量预测方法探讨》及《川东南地区基于海相泥页岩地层的压力系数预测_以丁山区块为例》提到了一种基于改进Fillippone法的海相页岩地层压力系数预测方法，但是仅仅单一考虑了纵波速度，构造稳定区取得了一定效果，但推广应用到构造复杂区，压力系数预测精度较低，该预测方法普适性较差。因此，有必要开发一种海相泥页岩地层的压力系数预测方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种海相泥页岩地层的压力系数预测方法及系统，其能够通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区。
根据本发明的一方面，提出了一种海相泥页岩地层的压力系数预测方法。所述方法可以包括：准备钻井数据、测井数据与地震资料；针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；基于所述钻井数据，分析钻井密度与所述纵波速度的关系，获得Gardner公式，基于所述关系和所述纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、所述纵波速度，拟合优化Fillippone公式中的系数；基于所述横波速度对所述优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据所述Gardner公式、所述钻井数据、所述纵波速度与所述横波速度，拟合再优化Fillippone公式中的系数；基于所述再优化Fillippone公式，预测地层压力系数。
优选地，针对所述超压页岩地层，以偶极子声波测井或全波列声波测井，分析对所述超压页岩地层敏感的测井参数，选择敏感测井曲线为所述纵波速度与所述横波速度。
优选地，所述Fillippone公式为：
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
优选地，所述优化Fillippone公式为：
其中，Pc 1为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
优选地，所述Gardner公式为：
其中，为上覆地层的平均密度，Vave为上覆地层平均速度，γ为经验系数，c1为地层平均速度的指数。
优选地，所述再优化Fillippone公式为：
其中，Pc 2为再优化地层压力系数，Vp为目的层纵波速度，Vs为横波速度，Vave为上覆地层平均速度，a2、b2、c2、d，e为再优化系数，其中a2、b2及d为纵波速度及横波速度的拟合经验系数，c2为Gardner公式中的地层平均速度的指数，e为多元统计及回归计算中的经验常数。
根据本发明的另一方面，提出了一种海相泥页岩地层的压力系数预测系统，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：准备钻井数据、测井数据与地震资料；针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；基于所述钻井数据，分析钻井密度与所述纵波速度的关系，获得Gardner公式，基于所述关系和所述纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、所述纵波速度，拟合优化Fillippone公式中的系数；基于所述横波速度对所述优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据所述Gardner公式、所述钻井数据、所述纵波速度与所述横波速度，拟合再优化Fillippone公式中的系数；基于所述再优化Fillippone公式，预测地层压力系数。
优选地，针对所述超压页岩地层，以偶极子声波测井或全波列声波测井，分析对所述超压页岩地层敏感的测井参数，选择敏感测井曲线为所述纵波速度与所述横波速度。
优选地，所述Fillippone公式为：
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
优选地，所述优化Fillippone公式为：
其中，Pc 1为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
本发明的有益效果在于：通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区，对于提高页岩气勘探成功率，推动南方海相页岩气勘探开发进程具有重要的战略意义。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的海相泥页岩地层的压力系数预测方法的步骤的流程图。
图2a和图2b分别示出了根据本发明的一个实施例的纵波速度与横波速度的测井响应特征的示意图。
图3示出了根据本发明的一个实施例的优化Fillippone公式、再优化Fillippone公式预测结果与实测结果的对比示意图。
图4示出了根据本发明的一个实施例的丁山地区过丁页1井-丁页3井-丁页2井五峰组-龙马溪组一段地层纵波速度反演剖面的示意图。
图5示出了根据本发明的一个实施例的丁山地区过丁页1井-丁页3井-丁页2井五峰组-龙马溪组一段地层横波速度反演剖面的示意图。
图6示出了根据本发明的一个实施例的丁山地区五峰组-龙马溪组一段页岩地层压力系数预测的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的海相泥页岩地层的压力系数预测方法的步骤的流程图。
在该实施例中，根据本发明的海相泥页岩地层的压力系数预测方法可以包括：
步骤101，准备钻井数据、测井数据与地震资料。
步骤102，针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；在一个示例中，针对超压页岩地层，以偶极子声波测井或全波列声波测井，分析对超压页岩地层敏感的测井参数，选择敏感测井曲线，获取纵波速度与横波速度。
步骤103，基于钻井数据，分析钻井密度与纵波速度的关系，获得Gardner公式，基于关系和纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、纵波速度，拟合优化Fillippone公式中的系数。
在一个示例中，Fillippone公式为：
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
在一个示例中，优化Fillippone公式为：
其中，Pc 1为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Fillippone公式中Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
在一个示例中，Gardner公式为：
其中，为上覆地层的平均密度，Vave为上覆地层平均速度，γ为经验系数，包含纵波速度，c1为地层平均速度的指数，实际数据拟合得到γ＝1.8，c1＝0.045。
步骤104，基于横波速度对优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据Gardner公式、所述钻井数据、纵波速度与横波速度，拟合再优化Fillippone公式中的系数，其中，所述钻井数据为实测地层压力系数。
在一个示例中，再优化Fillippone公式为：
其中，Pc 2为再优化地层压力系数，Vp为目的层纵波速度，Vs为横波速度，Vave为上覆地层平均速度，a2、b2、c2、d，e为再优化系数，其中a2、b2及d为纵波速度及横波速度的拟合经验系数，c2为Gardner公式中的地层平均速度的指数，e为多元统计及回归计算中的经验常数。
步骤105，基于再优化Fillippone公式，预测地层压力系数。
具体地，Fillippone法是由美国加利福尼亚联合石油公司的W.R.Fillippone提出的，他于1978年和1982年通过对墨西哥湾等地区的钻井、测井、地震等多方面资料的综合分析研究得出的不依赖于正常压实趋势线的计算公式，并在实际应用中取得了较好的效果，计算公式为公式(5)：
其中，Pp为地层压力，Pov为上覆地层压力，h为深度，g为重力加速度。将公式(5)以及静水柱压力计算公式(Pw＝ρwgh，ρw为地层水的相对密度)代入压力系数Pc的定义公式(Pc＝Pp/Pw)，得到Fillippone公式为公式(1)。
准备钻井数据、测井数据与地震资料，钻井数据主要为页岩气探井压力系数实测结果，测井数据包括偶极子声波测井或全波列声波测井资料，至少包括纵波时差、横波时差、密度测井曲线等，地震资料包括常规叠后或叠前偏移处理的成果数据、速度谱数据、构造解释层位及断层数据。
根据页岩地层超压机理研究表明，超压富有机质泥页岩地层具有较高的孔隙度，表现为“欠压实”特征。针对超压页岩地层，以偶极子声波测井或全波列声波测井资料，分析对页岩超压页岩储层敏感的测井参数，优选出敏感测井曲线，即为纵波速度曲线与横波速度曲线，进而获得纵波速度与横波速度。纵波速度与横波速度是通过叠前同时反演技术获取，利用优化后的高质量CRP叠前道集，选择Aki-Richards近似方程为反演方程，通过叠前同时反演技术获得目的层纵波速度及横波速度反演结果。
在测井响应特征分析的基础上，为便于地震预测实现，将公式(1)进行优化，基于海相页岩地层岩性岩相相对一致的特征，首先将目的地层最大、最小速度优化为单一系数a1、b1，同时基于Gardner公式，即公式(3)，由于Fillippone公式中Vmax及Vmin难以准确获取，且针对海相泥页岩地层岩相及岩性稳定的地质特征，对Vmax及Vmin参数优化为统一的经验系数a1及b1，将上覆地层的平均密度优化为平均速度的指数式，系数为c1，从而简化后针对海相泥页岩地层的优化Fillippone公式为公式(2)，根据地层实测压力系数、纵波速度通过多元统计回归算法，拟合优化Fillippone公式中的系数。其中，平均速度是通过模型层析速度建模方法获取的，首先通过CVI约束层速度反演，解决DIX反演速度横向不连续及纵向不稳定问题，能获得较为平滑的速度体，以此为初始模型，再结合构造解释模型及测井速度，采用反射波层析成像算法修改射线节点，反复迭代，最终得到稳定可靠的地层平均速度体Vave。
结合实测地层压力系数结果，利用公式(2)进行海相页岩地层压力系数计算的精度偏低，鉴于目的层纵波速度影响因素的复杂多样性，包括受地层埋深起伏的影响，进一步改进公式(2)，引入横波速度，增加参与压力系数计算的敏感信息，提高计算精度，获得再优化Fillippone公式为公式(4)。
图2a和图2b分别示出了根据本发明的一个实施例的纵波速度与横波速度的测井响应特征的示意图。
根据焦石坝区块及丁山区块10口重点探井测井曲线的统计分析，针对超压页岩地层，以偶极子声波测井或全波列声波测井资料，分析对页岩超压页岩储层敏感的测井参数，优选出敏感测井曲线，即为纵波速度曲线与横波速度曲线，进而获得纵波速度与横波速度的测井响应特征，如图2a与图2b所示，图中地层压力系数随着目的层纵波速度及横波速度增大，压力系数降低，压力系数大于1.2的超压地层为明显的低纵波速度及低横波速度信息。同时，从压力系数与纵波速度的交会分析发现，单纯的依据纵波速度的无法实现压力系数准确预测，尤其是丁页1、丁页3及焦页5，预测误差均大于0.3，因而优化Fillippone公式无法在构造复杂区(丁山地区及焦石坝焦页5井区)推广应用，因此，结合实测地层压力系数结果，利用公式(2)进行海相页岩地层压力系数计算的精度偏低。鉴于目的层纵波速度影响因素的复杂多样性，包括受地层埋深起伏的影响，进一步改进公式(2)，引入横波速度，增加参与压力系数计算的敏感信息，提高计算精度，获得再优化Fillippone公式为公式(4)。根据公式(3)，拟合公式(4)的系数c1，然后依据重点探井压力测试结果及纵横波速度信息，通过多元统计及回归的方法，拟合公式(4)中的系数a1、b1、d及e。基于公式(4)，预测地层压力系数。
本方法通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区，对于提高页岩气勘探成功率，推动南方海相页岩气勘探开发进程具有重要的战略意义。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
准备钻井数据、测井数据与地震资料，钻井数据主要为焦石坝区块及丁山区块10口重点探井(焦页1、焦页2、焦页3、焦页4、焦页5、焦页6、焦页7、丁页1、丁页2及丁页3井)，五峰组-龙马溪组地层压力系数实测结果；测井数据主要为纵波速度及横波速度；地震资料包括焦石坝区块及丁山区块常规叠后或叠前偏移处理的成果数据、速度谱数据、构造解释层位及断层数据。
在测井响应特征分析的基础上，为便于地震预测实现，将公式(1)进行优化，基于海相页岩地层岩性岩相相对一致的特征，首先将目的地层最大、最小速度优化为单一系数a1、b1，同时基于Gardner公式，即公式(3)，将上覆地层的平均密度优化为平均速度的指数式，系数为c1，从而简化后针对海相泥页岩地层的优化Fillippone公式为公式(2)，根据地层压力系数、纵波速度，拟合优化Fillippone公式中的系数为a1＝1.91583、b1＝0.00023、c1＝0.045。其中，平均速度是通过模型层析速度建模方法获取的，首先通过CVI约束层速度反演，解决DIX反演速度横向不连续及纵向不稳定问题，能获得较为平滑的速度体，以此为初始模型，再结合构造解释模型及测井速度，采用反射波层析成像算法修改射线节点，反复迭代，最终得到稳定可靠的地层平均速度体Vave。
结合实测地层压力系数结果，利用公式(2)进行海相页岩地层压力系数计算的精度偏低，鉴于目的层纵波速度影响因素的复杂多样性，包括受地层埋深起伏的影响，进一步改进公式(2)，引入横波速度，增加参与压力系数计算的敏感信息，提高计算精度，获得再优化Fillippone公式为公式(4)。根据公式(3)，拟合公式(4)的系数c2＝0.045，然后依据重点探井压力测试结果及纵横波速度信息，通过多元统计及回归的方法，拟合公式(4)中的系数为a2＝95.425、b2＝0.000222、d＝-0.001677，e＝-133.9418。基于公式(4)，预测地层压力系数。
图3示出了根据本发明的一个实施例的优化Fillippone公式、再优化Fillippone公式预测结果与实测结果的对比示意图，其中，灰线表示优化Fillippone公式预测结果，黑线表示再优化Fillippone公式预测结果，虚线表示实测结果，图中可明显看出再优化Fillippone公式预测精度更高，尤其是提高了复杂构造区钻井(丁页1、丁页3及焦页5)压力系数预测精度。
图4示出了根据本发明的一个实施例的丁山地区过丁页1井-丁页3井-丁页2井五峰组-龙马溪组一段地层纵波速度反演剖面的示意图。
图5示出了根据本发明的一个实施例的丁山地区过丁页1井-丁页3井-丁页2井五峰组-龙马溪组一段地层横波速度反演剖面的示意图。根据图4与图5，丁页1井到丁页2井地层的纵波速度与横波速度不断减小。
图6示出了根据本发明的一个实施例的丁山地区五峰组-龙马溪组一段页岩地层压力系数预测的示意图。由图可知，预测结果与实测吻合较好。由丁山地区东南部到盆内的西北部压力系数不断增大，其中丁页2井区压力系数最大，属于异常高压带，丁页2井压力系数预测1.52，实测1.49；丁山地区东南部大尺度高角度裂缝的发育破坏了页岩气保存条件，地层压力系数显著降低，丁页1及丁页3压力系数预测结果分别为0.9及1.1，实测分别为0.98及1.08。
综上所述，本发明通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区，对于提高页岩气勘探成功率，推动南方海相页岩气勘探开发进程具有重要的战略意义。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例，提供了一种海相泥页岩地层的压力系数预测系统，可以包括：存储器，存储有计算机可执行指令；处理器，处理器运行存储器中的计算机可执行指令，执行以下步骤：准备钻井数据、测井数据与地震资料；针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；基于钻井数据，分析钻井密度与纵波速度的关系，获得Gardner公式，基于关系和纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、纵波速度，拟合优化Fillippone公式中的系数；基于横波速度对优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据Gardner公式、所述钻井数据、纵波速度与横波速度，拟合再优化Fillippone公式中的系数；基于再优化Fillippone公式，预测地层压力系数。
在一个示例中，针对所述超压页岩地层，以偶极子声波测井或全波列声波测井，分析对所述超压页岩地层敏感的测井参数，选择敏感测井曲线为所述纵波速度与所述横波速度。
在一个示例中，Fillippone公式为：
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
在一个示例中，优化Fillippone公式为：
其中，Pc 1为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
本发明通过对页岩气保存条件进行定量评价研究，完善和深化页岩气甜点预测技术序列，优选页岩气高产富集区，对于提高页岩气勘探成功率，推动南方海相页岩气勘探开发进程具有重要的战略意义。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种海相泥页岩地层的压力系数预测方法，包括：
准备钻井数据、测井数据与地震资料；
针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；
基于所述钻井数据，分析钻井密度与所述纵波速度的关系，获得Gardner公式，基于所述关系和所述纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、所述纵波速度，拟合优化Fillippone公式中的系数；
基于所述横波速度对所述优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据所述Gardner公式、所述钻井数据、所述纵波速度与所述横波速度，拟合再优化Fillippone公式中的系数；
基于所述再优化Fillippone公式，预测地层压力系数。
2.根据权利要求1所述的海相泥页岩地层的压力系数预测方法，其中，针对所述超压页岩地层，以偶极子声波测井或全波列声波测井，分析对所述超压页岩地层敏感的测井参数，选择敏感测井曲线，获取所述纵波速度与所述横波速度。
3.根据权利要求1所述的海相泥页岩地层的压力系数预测方法，其中，所述Fillippone公式为：
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mover>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>o</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>min</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
4.根据权利要求3所述的海相泥页岩地层的压力系数预测方法，其中，所述优化Fillippone公式为：
<mrow>
<msubsup>
<mi>P</mi>
<mi>c</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
5.根据权利要求1所述的海相泥页岩地层的压力系数预测方法，其中，所述Gardner公式为：
<mrow>
<mover>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>o</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mi>&gamma;</mi>
<mo>*</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为上覆地层的平均密度，Vave为上覆地层平均速度，γ为经验系数，c1为地层平均速度的指数。
6.根据权利要求5所述的海相泥页岩地层的压力系数预测方法，其中，所述再优化Fillippone公式为：
<mrow>
<msubsup>
<mi>P</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>dV</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
</msubsup>
<mo>+</mo>
<mi>e</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为再优化地层压力系数，Vp为目的层纵波速度，Vs为横波速度，Vave为上覆地层平均速度，a2、b2、c2、d，e为再优化系数，其中a2、b2及d为纵波速度及横波速度的拟合经验系数，c2为Gardner公式中的地层平均速度的指数，e为多元统计及回归计算中的经验常数。
7.一种海相泥页岩地层的压力系数预测系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
准备钻井数据、测井数据与地震资料；
针对超压页岩地层，选择敏感测井曲线，获取纵波速度与横波速度；
基于所述钻井数据，分析钻井密度与所述纵波速度的关系，获得Gardner公式，基于所述关系和所述纵波速度对Fillippone公式进行优化，获取优化Fillippone公式，根据地层压力系数、所述纵波速度，拟合优化Fillippone公式中的系数；
基于所述横波速度对所述优化Fillippone公式进行再次优化，获取再优化Fillippone公式，根据所述Gardner公式、所述钻井数据、所述纵波速度与所述横波速度，拟合再优化Fillippone公式中的系数；
基于所述再优化Fillippone公式，预测地层压力系数。
8.根据权利要求7所述的海相泥页岩地层的压力系数预测系统，其中，针对所述超压页岩地层，以偶极子声波测井或全波列声波测井，分析对所述超压页岩地层敏感的测井参数，选择敏感测井曲线，获取所述纵波速度与所述横波速度。
9.根据权利要求7所述的海相泥页岩地层的压力系数预测系统，其中，所述Fillippone公式为：
<mrow>
<msub>
<mi>P</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mover>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>o</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>min</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Pc为地层压力系数，为上覆地层的平均密度，Vp为纵波速度，Vmax为孔隙为零时的地层速度，Vmin为刚性为零时的岩石速度。
10.根据权利要求9所述的海相泥页岩地层的压力系数预测系统，其中，所述优化Fillippone公式为：
<mrow>
<msubsup>
<mi>P</mi>
<mi>c</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为优化地层压力系数，Vp为纵波速度，Vave为上覆地层平均速度，a1、b1、c1为优化系数，其中a1及b1为Vmax及Vmin简化后的经验系数，c1为Gardner公式中的地层平均速度的指数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710860322.7A CN107817520B (zh) | 2017-09-20 | 2017-09-20 | 海相泥页岩地层的压力系数预测方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710860322.7A CN107817520B (zh) | 2017-09-20 | 2017-09-20 | 海相泥页岩地层的压力系数预测方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107817520A true CN107817520A (zh) | 2018-03-20 |
| CN107817520B CN107817520B (zh) | 2020-05-05 |
Family
ID=61607838
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710860322.7A Active CN107817520B (zh) | 2017-09-20 | 2017-09-20 | 海相泥页岩地层的压力系数预测方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107817520B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109188525A (zh) * | 2018-07-31 | 2019-01-11 | 中国地质大学（武汉） | 一种海相页岩层埋深数据的获取方法及系统 |
| CN110320574A (zh) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | 基于缓坡三角洲薄层砂体刻画的方法 |
| CN111060986A (zh) * | 2019-10-18 | 2020-04-24 | 中国石油化工股份有限公司 | 一种地层压力预测方法及岩性油藏评价方法 |
| CN112326803A (zh) * | 2020-09-17 | 2021-02-05 | 神华地质勘查有限责任公司 | 评估天然气储层可压性的方法及装置 |
| CN112649856A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 基于vsp资料的地层压力钻前预测方法及系统 |
| CN113514890A (zh) * | 2021-03-18 | 2021-10-19 | 中国石油大学(华东) | 一种使用地震数据预测地层压力系数的方法、装置及设备 |
| CN114893166A (zh) * | 2022-04-13 | 2022-08-12 | 中国石油大学(华东) | 一种地层压力系数计算方法 |
| CN115045646A (zh) * | 2022-06-07 | 2022-09-13 | 中国地质调查局油气资源调查中心 | 一种页岩气井位优选方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6826486B1 (en) * | 2000-02-11 | 2004-11-30 | Schlumberger Technology Corporation | Methods and apparatus for predicting pore and fracture pressures of a subsurface formation |
| CN104698492A (zh) * | 2013-12-09 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种计算异常地层压力的方法 |
| CN106199690A (zh) * | 2015-04-29 | 2016-12-07 | 中国石油化工股份有限公司 | 泥页岩裂缝的预测方法 |
| CN106324680A (zh) * | 2016-08-18 | 2017-01-11 | 中国石油天然气集团公司 | 地层破裂压力预测方法 |
| CN106814388A (zh) * | 2016-12-27 | 2017-06-09 | 中国石油大学(北京) | 砂泥岩储层地层压力的地震预测方法及装置 |
-
2017
- 2017-09-20 CN CN201710860322.7A patent/CN107817520B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6826486B1 (en) * | 2000-02-11 | 2004-11-30 | Schlumberger Technology Corporation | Methods and apparatus for predicting pore and fracture pressures of a subsurface formation |
| CN104698492A (zh) * | 2013-12-09 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种计算异常地层压力的方法 |
| CN106199690A (zh) * | 2015-04-29 | 2016-12-07 | 中国石油化工股份有限公司 | 泥页岩裂缝的预测方法 |
| CN106324680A (zh) * | 2016-08-18 | 2017-01-11 | 中国石油天然气集团公司 | 地层破裂压力预测方法 |
| CN106814388A (zh) * | 2016-12-27 | 2017-06-09 | 中国石油大学(北京) | 砂泥岩储层地层压力的地震预测方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 屈大鹏 等: "川东南地区基于海相泥页岩地层的压力系数预测——以丁山区块为例", 《物探与化探》 * |
| 王斌 等: "砂砾岩储层岩石物理特征及地层压力预测新模型", 《SPG/SEG北京2016国际地球物理会议电子文集》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110320574A (zh) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | 基于缓坡三角洲薄层砂体刻画的方法 |
| CN109188525A (zh) * | 2018-07-31 | 2019-01-11 | 中国地质大学（武汉） | 一种海相页岩层埋深数据的获取方法及系统 |
| CN109188525B (zh) * | 2018-07-31 | 2019-09-17 | 中国地质大学（武汉） | 一种海相页岩层埋深数据的获取方法及系统 |
| CN112649856A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 基于vsp资料的地层压力钻前预测方法及系统 |
| CN111060986A (zh) * | 2019-10-18 | 2020-04-24 | 中国石油化工股份有限公司 | 一种地层压力预测方法及岩性油藏评价方法 |
| CN112326803A (zh) * | 2020-09-17 | 2021-02-05 | 神华地质勘查有限责任公司 | 评估天然气储层可压性的方法及装置 |
| CN113514890A (zh) * | 2021-03-18 | 2021-10-19 | 中国石油大学(华东) | 一种使用地震数据预测地层压力系数的方法、装置及设备 |
| CN114893166A (zh) * | 2022-04-13 | 2022-08-12 | 中国石油大学(华东) | 一种地层压力系数计算方法 |
| CN114893166B (zh) * | 2022-04-13 | 2022-11-25 | 中国石油大学(华东) | 一种地层压力系数计算方法 |
| CN115045646A (zh) * | 2022-06-07 | 2022-09-13 | 中国地质调查局油气资源调查中心 | 一种页岩气井位优选方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107817520B (zh) | 2020-05-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817520A (zh) | 海相泥页岩地层的压力系数预测方法及系统 | |
| US11016214B2 (en) | Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium | |
| CN105445791B (zh) | 一种基于多种地震属性的地层孔隙压力预测方法 | |
| CN104914465B (zh) | 火山岩裂缝定量预测方法及装置 | |
| CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
| US20230083651A1 (en) | Method and system for analyzing filling for karst reservoir based on spectrum decomposition and machine learning | |
| EP3433643B1 (en) | Method and device for estimating sonic slowness in a subterranean formation | |
| CN111399044B (zh) | 一种储层渗透率预测方法、装置及存储介质 | |
| CN105652323B (zh) | 一种储层预测方法 | |
| CN108303510A (zh) | 页岩气储层性能的评价方法、装置及计算机存储介质 | |
| CN106597543A (zh) | 一种地层沉积相划分方法 | |
| CN107843927A (zh) | 基于井震联合速度的页岩地层压力预测方法及装置 | |
| CN107831540A (zh) | 储层物性参数直接提取新方法 | |
| CN112147677B (zh) | 油气储层参数标签数据生成方法及装置 | |
| CN104267430B (zh) | 确定地震流体敏感因子的方法及装置 | |
| CN104280773B (zh) | 利用随炮检距变化的时频谱交汇图预测薄层厚度的方法 | |
| RU2011148308A (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления | |
| CN110118994B (zh) | 一种基于地震反演和机器学习的陆相烃源岩定量预测方法 | |
| Naji et al. | Prediction of sonic shear wave using artificial neural network | |
| CN109339771A (zh) | 一种页岩油气层孔隙压力预测方法及系统 | |
| Bayuk | Why anisotropy is important for location of microearthquake events in shale? | |
| CN113138412A (zh) | 深层页岩气孔隙度地震预测方法及装置 | |
| CN112147676A (zh) | 一种煤层及夹矸厚度预测方法 | |
| Nath et al. | Prediction and analysis of geomechanical properties using deep learning: A Permian Basin case study | |
| Lee et al. | Basin‐scale prediction of S‐wave Sonic Logs using Machine Learning techniques from conventional logs |
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 |