CN102508310A - 一种煤田火区上部地层孔隙度分布的探测方法 - Google Patents
一种煤田火区上部地层孔隙度分布的探测方法 Download PDFInfo
- Publication number
- CN102508310A CN102508310A CN2011103351684A CN201110335168A CN102508310A CN 102508310 A CN102508310 A CN 102508310A CN 2011103351684 A CN2011103351684 A CN 2011103351684A CN 201110335168 A CN201110335168 A CN 201110335168A CN 102508310 A CN102508310 A CN 102508310A
- Authority
- CN
- China
- Prior art keywords
- porosity
- fire district
- coalfield
- survey line
- factor
- 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
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 37
- 238000001514 detection method Methods 0.000 title claims abstract description 16
- 239000003245 coal Substances 0.000 title abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims description 8
- 230000001939 inductive effect Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 239000004575 stone Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000011065 in-situ storage Methods 0.000 abstract 1
- 230000001902 propagating effect Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 23
- 238000002485 combustion reaction Methods 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种煤田火区上部地层孔隙度分布的探测方法,在煤田火区上部的地表面布置多条测线,在各测线间隔设置多个测点,逐一对每条测线上的测点施工钻孔,然后逐一在同一测线上相邻的两个钻孔中分别设置发射天线和接收天线,并依次对每条测线上的所有钻孔进行钻孔雷达多偏移距跨孔透射探测,采集在不同偏移距情况下电磁波从发射天线到接收天线的传播时间;通过计算机对所得电磁波的传播时间数据进行处理,并将每一条测线上相邻两钻孔之间二维剖面上的孔隙度数据分别进行拼接合并,作每条测线剖面上的地层孔隙度分布的等值线图,从而得出煤田火区上部地层孔隙度的分布,实现了在原位状态下快速、有效探测煤田火区上部空间内松散破碎地层的孔隙度分布情况。
Description
技术领域
本发明涉及一种探测方法,尤其是一种适用于煤田火区上部地层孔隙度分布的探测方法。
背景技术
煤田火区上部地层的孔隙度分布是研究煤田火区燃烧***的供风及出风通道以及灭火介质在火区上部松散地层中的运移规律的基础。但是目前还没有一种能够在原位状态下探测松散地层孔隙度分布的方法。地层孔隙度的测量在油气地质中应用较多,主要有测井法(如声波测井、密度测井和核磁共振测井)、地震反演法。但声波测井、密度测井和核磁共振测井法探测的范围非常有限,只能在有钻井的位置从纵向上对地层孔隙发育情况进行评价,而不能表现孔隙度在横向上的变化情况,因此不适用于空间上变化很大的煤田火区上部地层孔隙度的探测。地震反演法近年来在国内外油气地质勘探中得到了较多的应用,该方法通过岩石物理分析,建立声波阻抗或速度与储层孔隙度关系,从而求取孔隙度。地震反演法主要针对埋深较深的地层,但对于埋深较浅且较松散的火区上部松散区地层的分辨率较差,且该方法需要在测区进行放炮作业,会造成火区塌陷区的进一步塌陷,甚至导致人员伤亡。地震反演法结合测井资料联合反演地层孔隙度的方法也有应用,该方法是通过提取地震属性,建立地震属性与各测井孔隙度的关系,利用非线性神经网络和多元线性相关统计等,定量计算无井区的孔隙度分布从而得到全区的孔隙度变化的。但是由于所提取的地震属性与孔隙度很难直接建立有确定意义的关系式,所以对地层孔隙度的计算误差较大。通过以上分析可以看出,测井法和地震反演法均不能用于探测空间上变化较大且松散破碎地层的孔隙度分布。
发明内容
本发明的目的是提供一种探测煤田火区上部地层孔隙度分布的方法,以实现在原位状态下探测煤田火区上部地层孔隙度的分布。
本发明的煤田火区上部地层孔隙度分布的探测方法,包括如下步骤:
a、在煤田火区上部的地表面布置多条测线,在各测线间上隔设置多个测点,逐一对每条测线上的测点施工钻孔,然后逐一在同一测线上相邻的两个钻孔中分别设置发射天线和接收天线,通过发射天线和接收天线依次对每条测线上的所有钻孔进行钻孔雷达多偏移距跨孔透射探测,采集在不同偏移距情况下电磁波从发射天线到接收天线的传播时间;
b、将采集到的电磁波传播时间通过计算机处理,得到每条测线上相邻两个钻孔之间的二维剖面上各像元的慢度 ,根据公式: 计算出每条测线上相邻两个钻孔之间的二维剖面上各像元的体积介电常数;式中: 为电磁波在真空中的传播速度;
d、利用公式: 计算出在每相邻的两个钻孔之间的二维剖面上各像元的孔隙度,通过计算机将每一条测线上相邻两个钻孔之间二维剖面上的孔隙度数据分别进行拼接合并,并作每条测线剖面上的地层孔隙度分布的等值线图,从而得出煤田火区上部地层孔隙度的分布,孔隙度越大则地层越松散,孔隙度越小则地层越密实;式中: 为空气的介电常数,是一个与电场方向有关的因子。
所述在煤田火区上部的地表面布置多条测线的间距为50m;所述在各测线间隔设置多个测点的间距为40~50m。
有益效果:本发明在原位状态下采用钻孔雷达对火区上部地层进行多偏移距跨孔透射测量,采集在不同偏移距情况下电磁波从发射天线到接收天线的传播时间;通过计算机将所得数据进行处理,得到相邻两钻孔之间的二维剖面上的各像元的慢度及体积介电常数;采集火区上部未烧变地层的岩石样品,并测量其相对真空的介电常数;利用公式:计算出在相邻的两钻孔之间的二维剖面上各像元的孔隙度,通过计算机将每一条测线上相邻钻孔之间二维剖面上的孔隙度数据分别进行拼接合并,并作每条测线剖面上的地层孔隙度分布的等值线图,从而探测出煤田火区上部松散地层孔隙度的分布,孔隙度越大代表地层越松散,孔隙度越小则表示地层越密实。本发明方法简便、数据准确,与声波测井、密度测井和核磁共振测井等测井方法相比,克服了测井法不能表现孔隙度在横向上的变化情况的缺点;与地震反演法相比,克服了地震反演法只适用于探测埋深较深的地层孔隙度,而对于埋深较浅的地层的分辨率较差的缺点;另外,该方法无需建立地震属性与各井孔隙度的关系,因此避免了地震结合测井资料的联合反演法对地层孔隙度的计算误差较大的问题;同时,该方法无需在探测区进行放炮等危险作业,避免了人员伤亡等事故的发生。因此,本发明解决了已有测试方法的缺陷,实现了在原位状态下快速、有效探测煤田火区上部的空间上变化较大且松散破碎地层的孔隙度分布,具有广泛的适用性。
具体实施方式
实施例1:我国西部某煤田火区,火区面积100898m2,煤层均总厚度23m,倾角26°,因煤层燃烧(火区中心温度达1200℃),导致了地表煤层露头和上部地层整体沉降,形成了松散破碎带,为了解煤田火区燃烧***的供风、出风通道的分布以及灭火介质在火区上部松散地层中的运移规律,利用钻孔雷达对该火区上部地层孔隙度的分布进行了多偏移距跨孔透射探测,探测步骤如下:
a、在煤田火区上部的地表面布置多条测线,多条测线的间距为50m;在各测线上间隔设置多个测点,多个测点的间距为40~50m;逐一对每条测线上的测点施工钻孔,然后逐一在同一测线上相邻的两个钻孔中分别设置发射天线和接收天线,通过发射天线和接收天线依次对每条测线上的所有钻孔进行钻孔雷达多偏移距跨孔透射探测,采集在不同偏移距情况下电磁波从发射天线到接收天线的传播时间;
探测所用的钻孔雷达天线采用60MHz频率的偶极天线。探测时,首先将发射天线固定在孔口不动,接收天线从孔口位置向下移动,每3m采样一次,直至孔底;之后将发射天线向下移动3m,接收天线从下向上移动每3m采样一次;这样反复进行,直至发射天线移动到孔底,完成钻孔雷达多偏移距跨孔透射探测。
b、将采集到的电磁波传播时间通过计算机采用共轭梯度法进行迭代求解,得到每条测线上相邻两个钻孔之间的二维剖面上各像元的慢度,根据公式: 计算出每条测线上相邻两个钻孔之间的二维剖面上各像元的体积介电常数;式中: 为电磁波在真空中的传播速度;
d、利用公式: 计算出在每相邻的两个钻孔之间的二维剖面上各像元的孔隙度,通过计算机将每一条测线上相邻钻孔之间二维剖面上的孔隙度数据分别进行拼接合并,并作每条测线剖面上的地层孔隙度分布的等值线图,从而得出煤田火区上部地层孔隙度的分布,孔隙度越大则地层越松散,孔隙度越小则地层越密实;式中: 为空气的介电常数,是一个与电场方向有关的因子。
通过上述探测搞清煤田火区上部地层的孔隙度分布情况,为研究火区燃烧***的供风、出风通道以及在灭火施工过程中灭火介质在火区上部地层中的运移规律提供基础。
Claims (3)
1.一种煤田火区上部地层孔隙度分布的探测方法,其特征在于包括如下步骤:
a、在煤田火区上部的地表面布置多条测线,在各测线上间隔设置多个测点,逐一对每条测线上的测点施工钻孔,然后逐一在同一测线上相邻的两个钻孔中分别设置发射天线和接收天线,通过发射天线和接收天线依次对每条测线上的所有钻孔进行钻孔雷达多偏移距跨孔透射探测,采集在不同偏移距情况下电磁波从发射天线到接收天线的传播时间;
b、将采集到的电磁波传播时间通过计算机处理,得到每条测线上相邻两个钻孔之间的二维剖面上各像元的慢度 ,根据公式: 计算出每条测线上相邻两个钻孔之间的二维剖面上各像元的体积介电常数;式中: 为电磁波在真空中的传播速度;
c、采集煤田火区上部地层的未烧变岩石样品,用介电常数测量仪测量岩石样品相对真空的介电常数;
2.根据权利要求1中所述的煤田火区上部地层孔隙度分布的探测方法,其特征在于:所述在煤田火区上部的地表面布置多条测线的间距为50m。
3.根据权利要求1中所述的煤田火区上部地层孔隙度分布的探测方法,其特征在于:所述在各测线间隔设置多个测点的间距为40~50m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110335168 CN102508310B (zh) | 2011-10-31 | 2011-10-31 | 一种煤田火区上部地层孔隙度分布的探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110335168 CN102508310B (zh) | 2011-10-31 | 2011-10-31 | 一种煤田火区上部地层孔隙度分布的探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508310A true CN102508310A (zh) | 2012-06-20 |
CN102508310B CN102508310B (zh) | 2013-08-28 |
Family
ID=46220415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110335168 Expired - Fee Related CN102508310B (zh) | 2011-10-31 | 2011-10-31 | 一种煤田火区上部地层孔隙度分布的探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508310B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005706A (zh) * | 2015-07-31 | 2015-10-28 | 陕西煤业化工技术研究院有限责任公司 | 一种单一煤层火烧区边界的估算方法 |
CN105352996A (zh) * | 2015-10-20 | 2016-02-24 | 中国矿业大学 | 一种测试地下煤田火区覆岩温度变化的模型试验方法 |
CN105588870A (zh) * | 2016-01-04 | 2016-05-18 | 中国矿业大学 | 一种基于硫同位素组分比的煤田火区探测装置及方法 |
CN105628574A (zh) * | 2014-10-29 | 2016-06-01 | 中国石油天然气股份有限公司 | 岩样孔隙度的测量方法 |
CN109187643A (zh) * | 2018-06-25 | 2019-01-11 | 国网福建省电力有限公司宁德供电公司 | 一种基于广角测量法的输电杆塔塔基气隙检测装置 |
CN111596372A (zh) * | 2019-12-22 | 2020-08-28 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于电磁波测量***的含水孔隙度及矿化度反演方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种确定储层孔隙各向异性的瞬变电磁方法 |
US20110054795A1 (en) * | 2009-08-27 | 2011-03-03 | Conocophillips Company | Petrophysical Evaluation of Subterranean Formations |
US20110068788A1 (en) * | 2009-09-23 | 2011-03-24 | Chanh Cao Minh | Estimating porosity and fluid volume |
CN102052070A (zh) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性碎屑岩地层组分测井处理方法 |
-
2011
- 2011-10-31 CN CN 201110335168 patent/CN102508310B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种确定储层孔隙各向异性的瞬变电磁方法 |
CN101382599B (zh) * | 2007-09-03 | 2011-02-09 | 中国石油天然气集团公司 | 一种确定储层孔隙各向异性的瞬变电磁方法 |
US20110054795A1 (en) * | 2009-08-27 | 2011-03-03 | Conocophillips Company | Petrophysical Evaluation of Subterranean Formations |
US20110068788A1 (en) * | 2009-09-23 | 2011-03-24 | Chanh Cao Minh | Estimating porosity and fluid volume |
CN102052070A (zh) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性碎屑岩地层组分测井处理方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105628574A (zh) * | 2014-10-29 | 2016-06-01 | 中国石油天然气股份有限公司 | 岩样孔隙度的测量方法 |
CN105628574B (zh) * | 2014-10-29 | 2019-04-09 | 中国石油天然气股份有限公司 | 岩样孔隙度的测量方法 |
CN105005706A (zh) * | 2015-07-31 | 2015-10-28 | 陕西煤业化工技术研究院有限责任公司 | 一种单一煤层火烧区边界的估算方法 |
CN105005706B (zh) * | 2015-07-31 | 2018-06-29 | 陕西煤业化工技术研究院有限责任公司 | 一种单一煤层火烧区边界的估算方法 |
CN105352996A (zh) * | 2015-10-20 | 2016-02-24 | 中国矿业大学 | 一种测试地下煤田火区覆岩温度变化的模型试验方法 |
CN105352996B (zh) * | 2015-10-20 | 2018-08-10 | 中国矿业大学 | 一种测试地下煤田火区覆岩温度变化的模型试验方法 |
CN105588870A (zh) * | 2016-01-04 | 2016-05-18 | 中国矿业大学 | 一种基于硫同位素组分比的煤田火区探测装置及方法 |
CN105588870B (zh) * | 2016-01-04 | 2018-07-20 | 中国矿业大学 | 一种基于硫同位素组分比的煤田火区探测装置及方法 |
CN109187643A (zh) * | 2018-06-25 | 2019-01-11 | 国网福建省电力有限公司宁德供电公司 | 一种基于广角测量法的输电杆塔塔基气隙检测装置 |
CN111596372A (zh) * | 2019-12-22 | 2020-08-28 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于电磁波测量***的含水孔隙度及矿化度反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102508310B (zh) | 2013-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508310B (zh) | 一种煤田火区上部地层孔隙度分布的探测方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN102426384B (zh) | 一种探测地下采空区和岩溶分布的方法 | |
CN108957521B (zh) | 一种用于隧道长距离三维超前地质预报方法 | |
CN102288990B (zh) | 超浅层三维地震反射波勘探方法 | |
CN102866417A (zh) | 一种地下溶洞地震跨孔ct探测及层析成像装置及方法 | |
Elkarmoty et al. | In-situ GPR test for three-dimensional mapping of the dielectric constant in a rock mass | |
CN202837558U (zh) | 一种地下溶洞地震跨孔ct探测及层析成像装置 | |
CN108957548B (zh) | 一种多波多分量联合观测地震页岩气富集区预测方法 | |
CN103424769A (zh) | 一种采空区联合多波地震勘探方法 | |
CN107065014A (zh) | 一种基于微动探测结合地质钻探的地铁区间孤石群探测方法 | |
CN106772554A (zh) | 一种用于复杂地形条件下的多道瞬态面波勘探方法 | |
CN108490485A (zh) | 双线隧道微震监测***及其定位精度评估方法 | |
CN105093314B (zh) | 一种测定微地震震源的方法 | |
CN112882125B (zh) | 一种隧道-滑坡体综合探测和稳定性分析方法 | |
CN111812707A (zh) | 一种基于微动和面波勘探的地铁沿线场地内孤石探测方法 | |
CN102590874A (zh) | 一种山地采煤沉陷水田地表裂缝探测方法 | |
CN106371137A (zh) | 用于浅层异常体地震探测装置及其三维观测方法 | |
CN108614307B (zh) | 一种利用综合物探技术的仰拱质量检测方法及*** | |
Araffa et al. | Integration of geophysical techniques to detect geotechnical hazards: a case study in Mokattam, Cairo, Egypt | |
Long et al. | Locating geothermal resources using seismic exploration in Xian county, China | |
CN113050085A (zh) | 一种超前地质预报的方法 | |
CN116774276A (zh) | 一种基于量子测量的地震波采集方法与*** | |
Isakova et al. | GPR for mapping fractures for the extraction of facing granite from a quarry: A case study from Republic of Karelia | |
Nolan et al. | Near-surface salt dissolution void identification using passive MASW |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130828 |
|
CF01 | Termination of patent right due to non-payment of annual fee |