CN110133722A - 一种应用声速定性识别泥质分布类型新方法 - Google Patents

一种应用声速定性识别泥质分布类型新方法 Download PDF

Info

Publication number
CN110133722A
CN110133722A CN201910493729.XA CN201910493729A CN110133722A CN 110133722 A CN110133722 A CN 110133722A CN 201910493729 A CN201910493729 A CN 201910493729A CN 110133722 A CN110133722 A CN 110133722A
Authority
CN
China
Prior art keywords
sandstone
argillaceous
velocity
longitudinal wave
distribution
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
CN201910493729.XA
Other languages
English (en)
Other versions
CN110133722B (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.)
Karamay Onco Energy Technology Co Ltd
Original Assignee
Karamay Onco Energy Technology 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 Karamay Onco Energy Technology Co Ltd filed Critical Karamay Onco Energy Technology Co Ltd
Priority to CN201910493729.XA priority Critical patent/CN110133722B/zh
Publication of CN110133722A publication Critical patent/CN110133722A/zh
Application granted granted Critical
Publication of CN110133722B publication Critical patent/CN110133722B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/303—Analysis for determining velocity profiles or travel times
    • 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/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • 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

本发明提供了一种应用声速定性识别泥质分布类型新方法，涉及石油和天然气勘探技术领域，该方法包括利用CCT模型计算接触泥质砂岩的纵波速度作为砂岩纵波速度分布上限，绘制接触泥质、结构泥质和层状泥质的泥质含量与纵波速度交会图；利用HMHS模型计算分散泥质砂岩纵波速度作为砂岩纵波速度分布下限，绘制分散泥质含量与纵波速度交会图；标记分散泥质砂岩、接触泥质砂岩、结构泥质砂岩和层状泥质砂岩声速的分布区域，得到泥质砂岩声速分布图版；将实际泥质砂岩的声速和泥质含量数据代入图版，确定泥质分布形式。该方法能够定性识别泥质砂岩中泥质的分布类型，代表性强、周期短，能够大规模应用于测井、地震评价等。

Description

一种应用声速定性识别泥质分布类型新方法
技术领域
本发明涉及石油和天然气勘探技术领域，具体涉及一种应用声速定性识别泥质分布类型新方法。
背景技术
在实际的地球物理评价中需要考虑泥质分布类型和含量对砂岩声速的影响。通常，可以采用钻井取芯的实验室镜下观察、粒度分析来确定泥质分布类型，但该方法存在样本代表性有限、实验周期长等技术缺陷，无法大规模应用于测井、地震评价。
中国专利CN109324171A公开了一种基于岩性统计的沉积相定性识别方法，具体包括：数据收集、数据处理、参数计算、建立沉积相定性识别图版以及根据各计算单元的参数计算结果，其中，数据收集包括收集工区探井岩性录井数据和已知的沉积相数据；数据处理以单井上所划分的一种沉积相类型所对应的连续不间断的地层段作为一个计算单元，剔除厚度小于20m的计算单元，剔除每个计算单元内的无效岩层，剩余的岩层分为砂岩、泥岩两大类；参数计算分别计算各计算单元的砂岩频率以及砂地比两个参数；根据参数计算结果则在砂岩频率-砂地比交汇图上进行投点，再根据计算单元的沉积相类型以及投点位置，绘制不同沉积相类型的识别区域。该发明大幅度提高了沉积相识别工作的工作效率与准确率；具有速度快、成本低、效率高、主观干预影响小的优点。但该方法主要针对连续性地层段，并需要对无效岩层进行剔除，同时仍需要提前对岩层类型进行主观判断，应用范围有局限性。
中国专利CN106842359B公开了一种利用波阻抗定性识别复杂砂砾岩体岩性的方法，包括：利用岩心资料，结合测井数据，确定不同岩性在阵列感应、密度、声波时差、中子孔隙度这些测井曲线数据上的分布范围；进行岩性的粗分；进行岩性细分，并利用波阻抗，进一步确定出灰质砂岩、灰质砾岩；通过以上步骤，精细判别出七类岩性。该方法具有较好的操作性，利于推广，为复杂砂砾岩体岩性的定性识别提供了切实可行的新方法。但该方法需要详细的岩心资料，包括阵列感应、密度、声波时差和中子孔隙度等数据，需要针对岩层进行各个方面的参数进行测量，较为复杂。
目前为止，泥质分布类型识别方法存在较为复杂、需依赖多项参数以及主观性判断等问题，因此需要寻找一种代表性强、周期短、主观性依赖小，能够大规模应用于测井、地震评价的识别方法，用以定性识别泥质砂岩中泥质的分布类型。
发明内容
本发明针对现有技术存在的问题，提供了一种应用声速定性识别泥质分布类型新方法。该方法代表性强、周期短，能够大规模应用于测井、地震评价的识别方法，用以定性识别泥质砂岩中泥质的分布类型。
为实现上述目的，本发明采用的技术方案如下：
一种应用声速定性识别泥质分布类型新方法，包括以下步骤：
(1)、利用模型计算接触泥质砂岩的纵波速度作为岩纵波速度分布上限，绘制接触泥质、结构泥质和层状泥质的泥质含量与纵波速度交会图；
(2)、利用模型计算分散泥质砂岩纵波速度作为砂岩纵波速度分布下限，绘制分散泥质含量与纵波速度交会图；
(3)、分别标记分散泥质砂岩、接触泥质砂岩、结构泥质砂岩和层状泥质砂岩声速的分布区域，得到泥质砂岩声速分布图版；
(4)、将实际泥质砂岩的声速和泥质含量数据代入图版，并根据实际泥质砂岩声速数据和边界相对位置确定泥质分布形式。
进一步地，步骤(1)中所述模型为CCT模型，其中，CCT模型的使用方法参照论文Xue-Hui H,Jun-Xin G,Feng-Bi L,et al.Generalization of the Expression ofCementation Radius in Contact Cement Theory and Its Application[J].ChineseJournal of Geophysics,2015,57(4):439-449.
进一步地，步骤(2)中所述模型为HMHS模型，其中，HMHS模型的使用方法参照论文Dvorkin J P,Nur A M.Elasticity of high-porous sandstones:Theory for two NorthSea data sets[J].Geophysics,1996,61(5):1363-1370.
进一步地，步骤(2)中所述分散泥质的纵波速度介于Vlow-aVlow之间；步骤(1)中所述接触泥质的纵波速度介于cVup-Vup之间；步骤(1)中所述结构泥质和层状泥质的纵波速度均介于dVup-eVup之间。其中，所述Vup为砂岩纵波速度分布上限，Vlow为砂岩纵波速度分布下限。所述接触泥质的泥质含量小于b；所述结构泥质和层状泥质的泥质含量大于b，其中，b是一个泥质含量的插值，在声速分布上限Vup上，b以前认为是接触泥质，b以后认为是结构或层状泥质。
进一步地，所述a、b、c、d和e为待定系数，通过经验确定，在一些具体的实施方式中，所述a的范围为1-1.6，所述b的范围为0-0.15，所述c的范围为0.8-1，所述d的范围为0.8-1，所述e的范围为1-1.3。优选地，所述a为1.5，所述b为0.1，所述c为0.9，所述d为0.9，所述e为1.1。
值得说明的是，本发明被图版中所围起来的空白区域为混合泥质，实际声速分布越靠近下限说明泥质分布以分散泥质为主；实际声速越靠近上限说明泥质分布以结构泥质或层状泥质为主；实际声速越靠近接触泥质区域说明泥质分布以接触泥质为主。
本发明所取得的技术效果是：
1.本发明的方法能够定性识别泥质砂岩中泥质的分布类型，适用于泥质分布类型较为单一的情况。
2.本发明的方法代表性强、周期短、主观性依赖小，能够大规模应用于测井、地震评价，为地层精细解释提供参考，降低解释风险。
附图说明
图1为泥质砂岩声速分布图版；
图2为泥质砂岩声速分布及泥质砂岩声波速度图版，其中，散点为韩德华(Han,De-hua)公开发表的泥质砂岩声波速度数据。
具体实施方式
一种应用声速定性识别泥质分布类型新方法，包括以下步骤：
(1)、根据CCT模型计算接触泥质砂岩的纵波速度作为砂岩纵波速度分布上限，绘制接触泥质、结构泥质和层状泥质的泥质含量与纵波速度交会图；
(2)、根据HMHS模型计算分散泥质砂岩声速作为砂岩纵波速度分布下限，绘制分散泥质含量与纵波速度交会图；
(3)、标记分散泥质砂岩、接触泥质砂岩、结构泥质砂岩和层状泥质砂岩声速的分布区域，得到泥质砂岩声速分布图版；其中，分散泥质的纵波速度介于Vlow-aVlow之间，接触泥质的纵波速度介于cVup-Vup之间，泥质含量小于b；结构泥质或层状泥质的纵波速度介于dVup-eVup之间，泥质含量大于b。其中，a、b、c、d和e为待定系数，a为1.5，b为0.1，c为0.9，d为0.9，e为1.1；Vlow代表砂岩纵波速度分布下限，Vup代表砂岩纵波速度分布上限。依据如上参数建立的泥质砂岩声速分布图版如图1所示；其中，各个参数的设定为：石英体积模量K＝37GPa；石英剪切模量G＝44GPa；石英密度ρ＝2.65g/cm3；泥质体积模量Kc＝25GPa；泥质剪切模量Gc＝9GPa；泥质密度ρc＝2.55g/cm3；临界孔隙度配位数n＝9；泥质含量其中，为砂岩的实际孔隙度。
(4)、将论文Han D,Nur A,Morgan D.Effects of porosity and clay contenton wave velocities in sandstones[J].Geophysics,1986,51(11):2093.中实际泥质砂岩的声速和泥质含量数据代入图版(即图1)中，得到图2，并根据实际泥质砂岩声速数据和边界相对位置确定泥质分布形式。
由图2可知，本发明中，实际泥质砂岩的声速和泥质含量数据落入的区域大部分为结构或层状泥质，泥质砂岩样品中泥质分布以结构泥质或层状为主，泥质分布类型较为单一，表现出了近似线性的关系，和韩德华总结出近似线性的泥质砂岩声速的经验公式相互验证，上述结果表明本发明的方法得到的结论与韩德华提出来的经验公式和实验数据相互验证和补充，本发明中的方法较为准确。
最后应当说明的是，以上内容仅用以说明本发明的技术方案，而非对本发明保护范围的限制，本领域的普通技术人员对本发明的技术方案进行的简单修改或者等同替换，均不脱离本发明技术方案的实质和范围。

Claims (10)

1.一种应用声速定性识别泥质分布类型新方法，其特征在于：包括以下步骤：
(1)、利用模型计算接触泥质砂岩的纵波速度作为砂岩纵波速度分布上限，绘制接触泥质、结构泥质和层状泥质的泥质含量与纵波速度交会图；
(2)、利用模型计算分散泥质砂岩纵波速度作为砂岩纵波速度分布下限，绘制分散泥质含量与纵波速度交会图；
(3)、分别标记分散泥质砂岩、接触泥质砂岩、结构泥质砂岩和层状泥质砂岩声速的分布区域，得到泥质砂岩声速分布图版；
(4)、将实际泥质砂岩的声速和泥质含量数据代入图版，并根据实际泥质砂岩声速数据和边界相对位置确定泥质分布形式。
2.根据权利要求1所述的方法，其特征在于：步骤(1)中所述模型为CCT模型。
3.根据权利要求1所述的方法，其特征在于：步骤(2)中所述模型为HMHS模型。
4.根据权利要求1所述的方法，其特征在于：步骤(2)中所述分散泥质的纵波速度介于Vlow-aVlow之间，所述Vlow为砂岩纵波速度分布下限。
5.根据权利要求1所述的方法，其特征在于：步骤(1)中所述接触泥质的纵波速度介于cVup-Vup之间，所述Vup为砂岩纵波速度分布上限。
6.根据权利要求1所述的方法，其特征在于：步骤(1)中所述结构泥质的纵波速度介于dVup-eVup之间，所述Vup为砂岩纵波速度分布上限。
7.根据权利要求1所述的方法，其特征在于：步骤(1)中所述层状泥质的纵波速度介于dVup-eVup之间，所述Vup为砂岩纵波速度分布上限。
8.根据权利要求1所述的方法，其特征在于：步骤(1)中所述接触泥质的泥质含量小于b。
9.根据权利要求1所述的方法，其特征在于：步骤(1)中所述结构泥质和层状泥质的泥质含量大于b。
10.根据权利要求4-9任意一项所述的方法，其特征在于：所述a、b、c、d和e为待定系数，所述a的范围为1-1.6，所述b的范围为0-0.15，所述c的范围为0.8-1，所述d的范围为0.8-1，所述e的范围为1-1.3。
CN201910493729.XA 2019-06-06 2019-06-06 一种应用声速定性识别泥质分布类型的方法 Active CN110133722B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910493729.XA CN110133722B (zh) 2019-06-06 2019-06-06 一种应用声速定性识别泥质分布类型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910493729.XA CN110133722B (zh) 2019-06-06 2019-06-06 一种应用声速定性识别泥质分布类型的方法

Publications (2)

Publication Number Publication Date
CN110133722A true CN110133722A (zh) 2019-08-16
CN110133722B CN110133722B (zh) 2019-12-24

Family

ID=67580471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910493729.XA Active CN110133722B (zh) 2019-06-06 2019-06-06 一种应用声速定性识别泥质分布类型的方法

Country Status (1)

Country Link
CN (1) CN110133722B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987096A (zh) * 2021-03-15 2021-06-18 中国石油大学(华东) 一种高含泥质砂岩声速的计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426390A (zh) * 2011-10-21 2012-04-25 中国石油大学(北京) 一种非均质泥砂岩储层储量确定方法
CN104570127A (zh) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 一种利用地震波速度同时反演孔隙度和泥质含量的方法
US20150120197A1 (en) * 2013-10-30 2015-04-30 Chevron U.S.A. Inc. Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN104950331A (zh) * 2015-06-30 2015-09-30 中国石油天然气集团公司 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法
CN108399270A (zh) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 一种确定页岩地层中各向异性泥质比例的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426390A (zh) * 2011-10-21 2012-04-25 中国石油大学(北京) 一种非均质泥砂岩储层储量确定方法
CN104570127A (zh) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 一种利用地震波速度同时反演孔隙度和泥质含量的方法
US20150120197A1 (en) * 2013-10-30 2015-04-30 Chevron U.S.A. Inc. Method of Transforming Reservoir Properties to a Seismic Attribute for Hydrocarbon and Lithology Identification
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN104950331A (zh) * 2015-06-30 2015-09-30 中国石油天然气集团公司 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法
CN108399270A (zh) * 2017-02-08 2018-08-14 中国石油化工股份有限公司 一种确定页岩地层中各向异性泥质比例的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王飞: "交叉偶极声波测井资料处理及在致密气层评价中的应用研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987096A (zh) * 2021-03-15 2021-06-18 中国石油大学(华东) 一种高含泥质砂岩声速的计算方法
CN112987096B (zh) * 2021-03-15 2021-11-26 中国石油大学(华东) 一种高含泥质砂岩声速的计算方法

Also Published As

Publication number Publication date
CN110133722B (zh) 2019-12-24

Similar Documents

Publication Publication Date Title
CN106054248B (zh) 一种基于大面积致密储层地震岩石物理反演方法
CN109653725A (zh) 基于沉积微相和岩石相的混积储层水淹程度测井解释方法
CN105651966A (zh) 一种页岩油气优质储层评价方法及参数确定方法
CN105317435A (zh) 一种水平井裂缝识别方法
CN105301647B (zh) 区分灰质泥岩和砂岩的方法
CN107515290A (zh) 岩石矿物组分含量定量计算方法
CN105467466B (zh) 一种基于多尺度信息的致密储层成岩相预测方法
CN103412335A (zh) 一种利用地震物相体预测储层的方法
Johansson et al. Spectral induced polarization of limestones: time domain field data, frequency domain laboratory data and physicochemical rock properties
CN109324171A (zh) 一种基于岩性统计的沉积相定量识别方法
CN102914797B (zh) 一种获得地层各向异性系数的方法及装置
AU2021105676A4 (en) Near-source electromagnetic seismoelectric combined GR inversion method
CN110133722B (zh) 一种应用声速定性识别泥质分布类型的方法
Wang et al. Seismic diagenetic facies prediction of tight sandstone in the offshore sparse well area: An example from the Xihu Depression of the East China Sea Basin
JP2022134085A (ja) 弾性波地震周波数共振探査方法
CN110821496B (zh) 有机质页岩相模式建立方法及有机质页岩评价方法
CN112784404B (zh) 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法
CN113720745B (zh) 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法
CN1015671B (zh) 一种用地震资料描述油藏的方法
Zhang et al. Prediction of Oil Production in a Tight Sandstone Reservoir: Triassic Chang 9 Member, Jiyuan Area, Ordos Basin, NW China
Jouniaux et al. Seismo-electrics, electro-seismics, and seismo-magnetics for earth sciences
Rolfs Integrated geomechanical, geophysical, and geochemical analysis of the Bakken Formation, Elm Coulee field, Williston Basin, Montana
CN115508890B (zh) 一种裂缝孔隙型储层叠前叠后反演方法
RU2201606C1 (ru) Способ типизации и корреляции нефтегазопродуктивных горных пород по скважинным спектрально-временным параметрам
Zhou et al. Study on the standard of the four relationships of reservoir and the lower limit of effective reservoir in Nantun Formation of Oilfield A

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