CN106896401A - 一种研究岩石裂隙发育特征的测井曲线反演方法 - Google Patents

一种研究岩石裂隙发育特征的测井曲线反演方法 Download PDF

Info

Publication number
CN106896401A
CN106896401A CN201710115151.5A CN201710115151A CN106896401A CN 106896401 A CN106896401 A CN 106896401A CN 201710115151 A CN201710115151 A CN 201710115151A CN 106896401 A CN106896401 A CN 106896401A
Authority
CN
China
Prior art keywords
shale content
log
data
rock
inversion method
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.)
Pending
Application number
CN201710115151.5A
Other languages
English (en)
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.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN201710115151.5A priority Critical patent/CN106896401A/zh
Publication of CN106896401A publication Critical patent/CN106896401A/zh
Pending legal-status Critical Current

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

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

本发明公开了一种研究岩石裂隙发育特征的测井曲线反演方法，具体指岩石泥质含量的反演方法，包括：利用物探手段获取天然放射性测井曲线资料；根据天然放射性测井曲线导出岩石天然放射性实际量化数据；依据实际数据，利用公式反演泥质含量数据；根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线；分析泥质含量反演曲线，推断岩石裂隙发育程度。本发明将石油勘探领域中裂隙发育常用的测井反演方法应用到煤炭领域中，充分利用测井数据，更好反映矿区某深度范围内泥质含量的分布与大小，在一定程度上更直观表征裂隙发育情况，这为煤炭领域研究裂隙发育提供了一种新的思路和方法。

Description

一种研究岩石裂隙发育特征的测井曲线反演方法
技术领域
本发明涉及煤炭勘探技术领域，特别涉及一种岩石泥质含量的反演方法。
背景技术
获取岩石孔隙度的主要方法是取岩心测定和物探方法，目前国内外学者针对孔隙度与地震数据关系做了大量研究，比如利用伽马-伽马、视电阻率测井曲线来反演孔隙度，进而反映裂隙发育情况，但大多应用于石油勘探领域，很少将测井曲线反演应用到煤炭勘探技术领域。有少数学者在文章中提及用伽马-伽马、视电阻率测井曲线来反演孔隙度，但没有学者利用天然放射性曲线反演泥质含量，再利用泥质含量与孔隙发育的关系，来表征裂隙发育程度，本发明就是利用这一思路，充分利用测井数据资料，为煤炭领域研究裂隙发育提供了一种新的思路和方法。
发明内容
1.本发明的目的
在煤炭勘探技术领域中为研究裂隙发育提供一种新的思路和方法。
2.本发明的技术方案
为实现上述目的，本发明一种研究岩石裂隙发育特征的测井曲线反演方法，该方法包括：步骤A，利用物探手段获取天然放射性测井曲线资料；步骤B，根据天然放射性测井曲线导出岩石天然放射性实际量化数据；步骤C，依据实际数据，利用公式反演泥质含量数据；步骤D，根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线；步骤E，分析泥质含量反演曲线，推断岩石裂隙发育程度。
本发明的岩石泥质含量反演方法，首次将石油勘探领域中裂隙发育常用的测井反演方法应用到煤炭领域中，充分利用测井数据，更好反映矿区某深度范围内泥质含量的分布与大小，在一定程度上更直观表征裂隙发育情况，这为煤炭领域研究裂隙发育提供了一种新的思路和方法。
附图说明
此处用来说明的附图是为了对本发明的进一步解释和说明，为本申请的一部分，但并不能限定本发明。
图1为本发明岩石泥质含量的反演方法流程图。
图2A及图2B分别为潘谢外围矿区深、浅部根据反演公式及利用suffer软件做出的岩石平均泥质含量等值线图。
图3A、3B及3C为潘谢外围矿区二叠系煤系砂岩主要含水层重构泥质含量曲线反演图。
具体实施方案
下面结合附图及本发明所应用的实例，对本发明进行进一步说明。
图1为本发明岩石泥质含量的反演方法流程图。如图1所示，该方法包括：
步骤A，利用物探手段获取天然放射性测井曲线资料；
步骤B，根据天然放射性测井曲线导出岩石天然放射性实际量化数据；
步骤C，依据实际数据，利用公式反演泥质含量数据；
步骤D，根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线；
步骤E，分析泥质含量反演曲线，推断岩石裂隙发育程度。
以下结合具体实例，对上述每一步骤进行解释说明。
步骤A，利用物探手段获取天然放射性测井曲线资料。
步骤B，根据天然放射性测井曲线导出岩石天然放射性实际量化数据。在钻孔柱状图中截取某深度范围内天然放射性曲线，利用“getpl”插件，导出实际量化数据文本。
步骤C，依据实际数据，利用公式反演泥质含量数据。将岩石天然放射性实际量化数据文本导入excel表格，利用泥质含量Vsh与自然伽马相对值IGR之间的经验公式：
IGR＝(GR-GRmin)/(GRmax-GRmin) (1)
其中，Vsh为泥质含量，％；
IGR为自然伽马相对值；
GRmin为纯灰岩的自然伽马测井值，API；
GRmax为邻近泥岩层段的自然伽马测井值，API；
G为Hilchie指数，随地层的地质年代而改变(由经验值确定，对第三系以后的地层取3.7，对老地层取2.0)，注意IGR的取值范围为0～1。
在本次实例中，由于本实例研究对象主要为二叠系及以前的地层，因此本实例中反演G值取值为2.0，取自然伽马曲线最大值为GRmax，以免IGR计算为负值，GRmin为所提取的自然伽马离散数据的最低值。由图2A及图2B所示，利用反演公式将天然放射性曲线反演为泥质含量数据及利用suffer软件将泥质含量数据做出的岩石平均泥质含量等值线图。
步骤D，根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线。根据泥质含量数值，校正X坐标，与原Y坐标组合为修正后的坐标；将修正后坐标输入AutoCAD软件中，软件进行自动绘图，完成对泥质含量曲线的重构。如图3A、3B及3C所示，为潘谢外围矿区二叠系煤系砂岩主要含水层重构泥质含量曲线反演图。
步骤E，分析泥质含量反演曲线，推断岩石裂隙发育程度。根据泥质含量反演曲线的波动程度和分布范围，分析不同区域下岩石裂隙发育程度。
根据图2A、图2B对比所示，潘谢外围深浅部裂隙发育特征总体上呈现一定的差异性，在陈桥-潘集背斜南翼一侧，外围区域平均泥质含量整体表现出由西向东逐渐增大的趋势，在背斜北翼一侧，区域平均泥质含量表现出不同的变化趋势，由南向北开始逐渐减小。总体来说，研究区深部岩石泥质含量普遍较小于浅部岩石，体现出深部泥质岩类含量较少，岩性以砂岩为主，裂隙发育程度较高，侧面反映出深部富水程度较高于浅部区域。
由图3A、图3B、图3C可知，13-1煤与11-2煤顶底板含水层泥质含量随埋深的增加，其变化趋势基本一致，总体上呈现先减小后增大的趋势，8～4～1煤顶底板含水层泥质含量的变化则表现一定的差异性，随着深度的增加，其泥质含量呈现增大-减小-增大的趋势。综合上述曲线分析，可以看出随着深度的逐渐增加，13-1煤、11-2煤、8～4～1煤顶底板含水层泥质含量整体上是一个增大的过程，侧面反映岩石裂隙发育程度逐渐降低。
本发明的岩石泥质含量反演方法，首次将石油勘探领域中裂隙发育常用的测井反演方法应用到煤炭领域中，充分利用测井数据，更好反映矿区某深度范围内泥质含量的分布与大小，在一定程度上更直观表征裂隙发育情况，这为煤炭领域研究裂隙发育提供了一种新的思路和方法。
以上所述具体实例，对本发明的目的，过程和有益效果进行详细说明，并不用于限定本发明的限定范围，凡在本发明的精神原则之内，所做的任何修改、等同替换等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种研究岩石裂隙发育特征的测井曲线反演方法，其特征在于，该方法包括：
步骤A，利用物探手段获取天然放射性测井曲线资料；
步骤B，根据天然放射性测井曲线导出岩石天然放射性实际量化数据；
步骤C，依据实际数据，利用公式反演泥质含量数据；
步骤D，根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线；
步骤E，分析泥质含量反演曲线，推断岩石裂隙发育程度。
2.如权利要求1所述的岩石裂隙发育特征的测井曲线反演方法，其特征在于，步骤B中，根据天然放射性测井曲线导出岩石天然放射性实际量化数据。在钻孔柱状图中截取某深度范围内天然放射性曲线，输入“getpl”指令，导出实际量化数据文本。
3.如权利要求1所述的岩石裂隙发育特征的测井曲线反演方法，其特征在于，步骤C中，依据实际数据，利用公式反演泥质含量数据。将岩石天然放射性实际量化数据文本导入excel表格，利用泥质含量Vsh与自然伽马相对值IGR之间的经验公式：
IGR＝(GR-GRmin)/(GRmax-GRmin) (1)
V s h = 100 × ( 2 G × I G R - 1 ) / ( 2 G - 1 ) - - - ( 2 )
其中，Vsh为泥质含量，％；
IGR为自然伽马相对值；
GRmin为纯灰岩的自然伽马测井值，API；
GRmax为邻近泥岩层段的自然伽马测井值，API；
G为Hilchie指数，随地层的地质年代而改变(由经验值确定，对第三系以后的地层取3.7，对老地层取2.0)，注意IGR的取值范围为0～1。
4.如权利要求1所述的岩石裂隙发育特征的测井曲线反演方法，其特征在于，步骤D中，根据泥质含量数据，利用AutoCAD软件重构泥质含量曲线。根据泥质含量数值，校正X坐标，与原Y坐标组合为修正后的坐标；将修正后坐标输入AutoCAD软件中，软件进行自动绘图，完成对泥质含量曲线的重构。
CN201710115151.5A 2017-03-01 2017-03-01 一种研究岩石裂隙发育特征的测井曲线反演方法 Pending CN106896401A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710115151.5A CN106896401A (zh) 2017-03-01 2017-03-01 一种研究岩石裂隙发育特征的测井曲线反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710115151.5A CN106896401A (zh) 2017-03-01 2017-03-01 一种研究岩石裂隙发育特征的测井曲线反演方法

Publications (1)

Publication Number Publication Date
CN106896401A true CN106896401A (zh) 2017-06-27

Family

ID=59184384

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710115151.5A Pending CN106896401A (zh) 2017-03-01 2017-03-01 一种研究岩石裂隙发育特征的测井曲线反演方法

Country Status (1)

Country Link
CN (1) CN106896401A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107748814A (zh) * 2017-10-13 2018-03-02 中国煤炭地质总局七三勘探队 在平面图上钻孔方法
CN108446532A (zh) * 2018-03-20 2018-08-24 安徽理工大学 一种基于层次分析与gis的底板岩层阻、隔水性评价方法
CN109598049A (zh) * 2018-11-28 2019-04-09 中化地质矿山总局地质研究院 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105114068A (zh) * 2015-09-07 2015-12-02 中国地质大学（北京） 一种利用测井资料预测煤层气区块高产水区的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105114068A (zh) * 2015-09-07 2015-12-02 中国地质大学（北京） 一种利用测井资料预测煤层气区块高产水区的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高丽军: "柳州地区太原组灰岩含水富水性研究", 《2013年煤层气学术研讨会论文集》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107748814A (zh) * 2017-10-13 2018-03-02 中国煤炭地质总局七三勘探队 在平面图上钻孔方法
CN108446532A (zh) * 2018-03-20 2018-08-24 安徽理工大学 一种基于层次分析与gis的底板岩层阻、隔水性评价方法
CN109598049A (zh) * 2018-11-28 2019-04-09 中化地质矿山总局地质研究院 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法
CN109598049B (zh) * 2018-11-28 2023-04-07 中化地质矿山总局地质研究院 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法

Similar Documents

Publication Publication Date Title
Eaton et al. Breakdown of the Gutenberg‐Richter relation for microearthquakes induced by hydraulic fracturing: Influence of stratabound fractures
Ding et al. Fractures in continental shale reservoirs: a case study of the Upper Triassic strata in the SE Ordos Basin, Central China
Seithel et al. Probability of fault reactivation in the Bavarian Molasse Basin
Toro et al. Thermal and kinematic evolution of the Eastern Cordillera fold and thrust belt, Colombia
Beamish et al. The Cornubian geothermal province: heat production and flow in SW England: estimates from boreholes and airborne gamma-ray measurements
Tavani et al. Transpressive inversion of a Mesozoic extensional forced fold system with an intermediate décollement level in the Basque‐Cantabrian Basin (Spain)
Oncken et al. Strain accumulation controls failure of a plate boundary zone: Linking deformation of the Central Andes and lithosphere mechanics
CN106896401A (zh) 一种研究岩石裂隙发育特征的测井曲线反演方法
Masaferro et al. Kinematic evolution and fracture prediction of the Valle Morado structure inferred from 3-D seismic data, Salta province, northwest Argentina
Frost III et al. Platform-margin trajectory as a control on syndepositional fracture patterns, Canning Basin, Western Australia
Wigger et al. Logging solutions for completion optimization in unconventional resource plays
Zhu et al. Seismic characterization of hypogenic karst systems associated with deep hydrothermal fluids in the middle‐lower ordovician Yingshan Formation of the Shunnan Area, Tarim Basin, NW China
Han et al. Well log evaluation of shale gas reservoirs and preservation conditions of Lower Cambrian shale succession in Cengong Block of southeast Sichuan basin, south China
Giacomuzzi et al. Tracking transient changes in the plumbing system at Campi Flegrei Caldera
Mercier et al. Evolution of a block cave from time-lapse passive source body-wave traveltime tomography
Wang et al. Exploration process and genesis mechanism of deep geothermal resources in the North Jiangsu Basin, East China: From nothing to something
Sun et al. Characteristics and prediction of weathered volcanic rock reservoirs: A case study of Carboniferous rocks in Zhongguai paleouplift of Junggar Basin, China
Wang et al. Tectonic controls on the Pleistocene–Holocene Wudalianchi volcanic field (northeastern China)
CN112363210A (zh) 基于透射槽波波速和衰减系数联合反演的煤厚定量预测方法
Kaven et al. Micro-seismicity, fault structure and hydraulic compartmentalization within the Coso Geothermal Field, California
Thompson et al. Using core and outcrop analogues to predict flow pathways in the subsurface: examples from the Triassic sandstones of north Cheshire, UK
Wang et al. Compactive deformation of incoming calcareous pelagic sediments, northern Hikurangi subduction margin, New Zealand: Implications for subduction processes
Wang et al. Heat flow correction for the high-permeability formation: A case study for xiong’an new area
Ruiz-Constán et al. Gravity anomalies and orthogonal box fold development on heterogeneous basement in the Neogene Ronda Depression (Western Betic Cordillera)
Jiang et al. The characteristics of the neotectonic movement and their effects on the formation of gas reservoirs in the marginal depressions of Tarim basin, NW China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB03 Change of inventor or designer information

Inventor after: Liu Qimeng

Inventor after: Li Pengfei

Inventor after: Chen Xiuyan

Inventor after: Ye Mei

Inventor after: Zhang Dandan

Inventor before: Li Pengfei

Inventor before: Chen Xiuyan

Inventor before: Liu Qimeng

Inventor before: Ye Mei

Inventor before: Zhang Dandan

CB03 Change of inventor or designer information
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170627

RJ01 Rejection of invention patent application after publication