CN114487125B - 一种用于确定煤体各向异性的三维监测综合方法 - Google Patents
一种用于确定煤体各向异性的三维监测综合方法 Download PDFInfo
- Publication number
- CN114487125B CN114487125B CN202210064568.4A CN202210064568A CN114487125B CN 114487125 B CN114487125 B CN 114487125B CN 202210064568 A CN202210064568 A CN 202210064568A CN 114487125 B CN114487125 B CN 114487125B
- Authority
- CN
- China
- Prior art keywords
- layer
- coal
- drilling
- ultrasonic
- penetrating
- 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
Links
- 239000003245 coal Substances 0.000 title claims abstract description 136
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012544 monitoring process Methods 0.000 title claims abstract description 22
- 238000005553 drilling Methods 0.000 claims abstract description 162
- 230000035699 permeability Effects 0.000 claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 20
- 239000011148 porous material Substances 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 230000035945 sensitivity Effects 0.000 claims abstract description 14
- 238000009533 lab test Methods 0.000 claims abstract description 12
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 11
- 238000010998 test method Methods 0.000 claims abstract description 10
- 230000002159 abnormal effect Effects 0.000 claims abstract description 6
- 208000008918 voyeurism Diseases 0.000 claims abstract description 6
- 238000012216 screening Methods 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 claims description 56
- 238000012360 testing method Methods 0.000 claims description 34
- 239000007788 liquid Substances 0.000 claims description 33
- 230000000149 penetrating effect Effects 0.000 claims description 24
- 239000002775 capsule Substances 0.000 claims description 23
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 238000007789 sealing Methods 0.000 claims description 14
- 239000011435 rock Substances 0.000 claims description 11
- 238000002347 injection Methods 0.000 claims description 9
- 239000007924 injection Substances 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 239000007822 coupling agent Substances 0.000 claims description 8
- 238000005119 centrifugation Methods 0.000 claims description 6
- 238000013499 data model Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000007872 degassing Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract 1
- 230000006978 adaptation Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/11—Analysing solids by measuring attenuation of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/02—Devices for withdrawing samples
- G01N1/04—Devices for withdrawing samples in the solid state, e.g. by cutting
- G01N1/08—Devices for withdrawing samples in the solid state, e.g. by cutting involving an extracting tool, e.g. core bit
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Theoretical Computer Science (AREA)
- Dispersion Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Biomedical Technology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种用于确定煤体各向异性的三维监测综合方法，采用现场探测和实验室测试相结合的方式，在实验室测试方面借助低场核磁共振技术、应力孔渗测试方法和三向纵波测试方法对不同深度煤芯的孔隙结构、应力、渗透率敏感性进行测试并获取对应数据；在现场探测方面通过钻孔窥视仪获取对不同深度钻孔壁面上的层理或弱面结构位置及周向分布数据，并通过特定的穿层钻孔位置布设及超声波收发系统的探测方式及位置布设，从而能获得探测钻孔间不同深度煤体基质信息，然后上述进行现场探测数据和实验室测试获得数据结果通过重叠数据筛选、异常数据剔除及综合数据反演后，计算得出煤层不同深度位置的非均质性系数，最终获得煤层各向异性的评估准则。
Description
技术领域
本发明涉及一种用于确定煤体各向异性的三维监测综合方法，尤其适用于特厚煤层全局结构分布及空间三维各向异性的精细表征。
背景技术
随着矿井开采深度的不断增加，煤层厚度和地应力也随之增大。煤基质的属性精准监测成为深部较厚煤层安全开采的首要前提。煤层属性主要包括煤体变质程度、煤体孔隙率、煤体矿物组分、煤体渗透性等，这些属性参数均与各向异性存在密切关系。比如，煤体变质程度低，基质成熟度相对较低，孔隙率相对较大，形成煤体在孔隙结构、受力方向、渗透率等方面的各向异性。煤体基质的各向异性对煤体瓦斯含量分布及渗流运移同样也会产生一定程度的影响。孔隙结构连通性较差时，容易在煤基质内部形成多个独立的“瓦斯包”，可能会对煤层瓦斯含量及压力测试造成一定程度的影响，同时，各向异性较大的煤体基质在地应力的作用下，采用强化致裂增透措施所引起的裂隙发育、起裂具有很大的方向随机性，很难做到裂隙发育的方向可控。因此，煤体各向异性的确定对于探明煤体资源、瓦斯含量分布及宏观裂隙结构具有重要意义。
目前主要通过钻探和物探的方式对煤体各向异性进行测试，其中钻探方式是指通过钻取大量岩心，并根据对岩心孔隙结构等基本性质进行测试确定煤层的各向异性，但是由于该方式样品取样位置随机性较强，且样品取芯位置不能代表煤层整体，因此获得数据对于煤层整体来说存在一定误差。同时，对岩心的测试手段主要包括压汞法、SEM、CT扫描法等，这些测试手段对样品尺寸具有一定的要求，测试范围较为有限，不能够反映所有的地层信息。而物探主要是借助巷道超前探测仪等检测手段对煤层结构进行探测，其探测范围较广，但是由于其施工地点一般布置在煤巷内，使得物探技术的信号接收容易受到煤层矿物质、气体、水等因素的影响发生较大衰减，接收信号强度及测试精度存在较大的误差。同时监测地点布置及探头的下放等都具有距离要求，震源与探头之间及相邻探头之间若距离较远，则接收信号强度较小，可能会因为门槛值设置遗漏掉某种信号的重要信息；而若距离较近，则会使测试工作量和测试成本大幅增加。
因此，如何提供一种方法，能在提高对煤体各向异性的监测精度及增加对煤体的探测范围情况下，还能有效节省监测成本，为煤层、页岩气储层、油气储层的非均质性探测提供数据支撑，是本行业的研究方向之一。
发明内容
针对上述现有技术存在的问题，本发明提供一种用于确定煤体各向异性的三维监测综合方法，能在提高对煤体各向异性的监测精度及增加对煤体的探测范围情况下，还能有效节省监测成本，为煤层、页岩气储层、油气储层的非均质性探测提供数据支撑。
为了实现上述目的，本发明采用的技术方案是：一种用于确定煤体各向异性的三维监测综合方法，具体步骤为：
A、从岩巷顶板向待测厚煤层钻设多组穿层钻孔组，且各个穿层钻孔组的钻孔角度不同，每个穿层钻孔组均由多个穿层钻孔组成，且同组内各个穿层钻孔的钻孔方向平行；各穿层钻孔组均穿过煤层厚度直至煤层顶部与岩体结合处，然后利用岩心取芯机对各穿层钻孔组处于煤层内不同深度的位置进行取芯，接着将获取的煤芯作好标记并送至样品养护室内进行定期养护；
B、利用钻孔窥视仪对各穿层钻孔组取芯位置的钻孔壁依次进行可视化探测，且在取芯位置每移动0.1m，输出一次拍摄的钻孔壁图像，最后通过多次图像集成，获取整个取芯位置的钻孔壁图像，进而识别出钻孔壁上存在的层理、割理结构信息；
C、选择一组穿层钻孔组，将超声波发射装置放入该组内其中一个穿层钻孔直至钻孔最深处，所述超声波发射装置上设有多个超声波发射探头，超声波发射探头的数量与该组内其余穿层钻孔的数量相同、且各个超声波发射探头在该穿层钻孔内分别一一对应朝向该组内其余穿层钻孔，该组内其余穿层钻孔分别放入一个超声波接收探头，各个超声波接收探头分别达到各个穿层钻孔的最深处、且各个超声波接收探头均朝向放入超声波发射装置的穿层钻孔；然后对该组内所有穿层钻孔均采用胶囊封孔器封孔，使各个穿层钻孔内在靠近孔口的位置形成胶囊封孔段；所述胶囊封孔器上装有穿过胶囊封孔器的进液管和回液管，进液管一端和回液管一端均设有阀门、且回液管另一端伸入至穿层钻孔的最深处；打开进液管和回液管的阀门，采用注入泵将耦合剂以一定压力通过进液管分别注入各个穿层钻孔内，直至回液管一端流出耦合剂时，关闭回液管的阀门，并持续保压一段时间后停止耦合剂注入，并将进液管的阀门关闭，使各个穿层钻孔内均分成耦合剂充填段、胶囊封孔段和剩余钻孔段；
D、开启超声波发射装置，使其各个超声波发射探头同时向该组内其余穿层钻孔发射频率相同的超声波，该组内其余穿层钻孔内各个超声波接收探头分别接收各自对应的超声波发射探头发出的超声波信号，均传递给超声波数据分析仪，超声波数据分析仪对超声波信号衰减信息进行反演，从而得出超声波发射装置所处穿层钻孔与其余各个穿层钻孔间在该深度位置的煤体基质信息；接着将超声波发射装置及各个超声波接收探头在各自穿层钻孔内均向孔口方向移动相同距离，并重复上述超声波发射及接收过程，从而得出超声波发射装置所处穿层钻孔与其余各个穿层钻孔间在当前深度位置的煤体基质信息；如此重复，直至超声波发射装置及各个超声波接收探头在各自穿层钻孔内到达胶囊封孔器时，停止超声波监测过程，进而获得步骤A所取该穿层钻孔组内各个不同深度位置煤芯周围的煤体基质信息；对其余穿层钻孔组依次重复步骤C和D，获得各穿层钻孔组内各个不同深度位置煤芯周围的煤体基质信息；
E、将步骤A标记的各穿层钻孔组煤芯分别采用低场核磁共振技术、应力孔渗测试方法及三向纵波测试方法对煤芯的孔隙结构、应力敏感性及渗透率各向异性进行实验室测试，获取各穿层钻孔组不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据；
F、将步骤B和D通过现场测试获得的数据信息与步骤E通过实验室测试获得数据结果进行综合分析，将所有数据通过重叠数据筛选、异常数据剔除处理，将筛选后的数据随机分成训练组和测试组，通过BP优化算法对训练组进行训练后创建BP神经网络，接着利用创建的BP神经网络对测试组进行测试，最终经过测试后获得非均质特性数据模型，后续获得不同深度位置的参数输入至该数据模型即能计算得到对应位置的非均质特性，从而确定煤体各向异性。
进一步，所述步骤A中每个穿层钻孔组均由9个穿层钻孔组成，且9个穿层钻孔呈3×3矩形阵列布置，相邻钻孔之间的间距为1m；每个穿层钻孔直径为150mm；取芯时在各个穿层钻孔内每隔0.5m取一次煤芯，直至岩心取芯机达到煤层顶部与岩体结合处时，停止取芯；所述钻取煤芯的尺寸为直径25mm、高度为50～60mm的圆柱体和直径80mm～100mm的圆柱体。
进一步，所述步骤C中胶囊封孔段的长度为穿层钻孔长度的三分之一；所述耦合剂的注入压力为2～3MPa。
进一步，所述超声波发射装置上设有8个超声波发射探头，8个超声波发射探头能同时释放同等频率的超声波信号；各个超声波发射探头与其对应的超声波接收探头之间的对向角度偏移量不超过5°，超声波发射装置放入各穿层钻孔组内中心位置的穿层钻孔内；所述步骤D中超声波发射装置及各个超声波接收探头在各自穿层钻孔内每次移动距离为0.5m。
进一步，所述步骤E的具体过程为：
对钻取煤芯分别进行真空脱水脱气及完全饱水处理，利用低场核磁共振技术测试得到饱水状态下的T2弛豫分布曲线；对饱水煤芯进行离心处理，同样利用低场核磁共振技术测试获得离心状态下的T2弛豫分布曲线；通过饱水-离心联测方法处理T2弛豫幅值差后，获得自由水空间占比，通过T2弛豫幅值差与孔隙度之间的线性关系，得到饱水、离心两种状态下的总孔隙率和残余孔隙率，二者之差等于有效孔隙率，最后根据SDR模型计算得到不同深度煤芯的核磁渗透率；
将煤芯置于应力三向非均质性的夹持器腔体中，通过施加不同荷载，获得煤芯在三个不同方向上的应力、渗透率敏感性数据；
对煤芯三个不同正交方向上不同测点位置的纵波进行测试，通过计算声波各向异性系数、波速比值、波列图分布参数获得煤体结构各向异性特征；其中，声波各向异性系数：
波速比值：
最终得到各穿层钻孔组不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据。
与现有技术相比，本发明采用现场探测和实验室测试相结合的方式，对煤体各向异性特征进行全面表征；在实验室测试方面借助低场核磁共振技术、应力孔渗测试方法和三向纵波测试方法对不同深度煤芯的孔隙结构、应力、渗透率敏感性进行测试并获取对应数据，揭示不同深度煤层孔隙结构空间分布差异性、三向渗透率和应力受载非均一性与煤体各向异性之间的对应关系；在现场探测方面通过钻孔窥视仪获取对不同深度钻孔壁面上的层理或弱面结构位置及周向分布数据，并通过特定的穿层钻孔位置布设及超声波收发系统的探测方式及位置布设，从而能获得探测钻孔间不同深度煤体基质信息，该信息包括不同深度煤体基质断层、大型割理层理、溶洞数据，然后上述进行现场探测数据和实验室测试获得数据结果进行综合分析，将所有数据通过重叠数据筛选、异常数据剔除及综合数据反演后，计算得出煤层不同深度位置的非均质性系数，最终获得煤层各向异性的评估准则。综上所述，本发明通过对同一钻孔先进行不同深度煤芯钻取，接着对该钻孔进行二次利用，通过钻孔窥视技术和超声波探测技术进行煤层探测，无需再进行钻孔作业，不仅能提高对煤体各向异性的监测精度及增加对煤体的探测范围，而且能有效节省监测成本，为煤层、页岩气储层、油气储层的非均质性探测提供数据支撑。
附图说明
图1是本发明的整体布设示意图。
图2是本发明中穿层钻孔组内9个穿层钻孔的空间布置图。
图3是本发明中进行封孔及耦合剂注入后各个穿层钻孔的内部示意图。
图4是本发明中超声波发射装置的整体结构示意图。
图5是本发明中超声波发射及接收的原理结构图。
图中：1-岩巷；2-待测厚煤层；3-穿层钻孔组；3-1-中心位置的穿层钻孔；3-2-其余穿层钻孔；4-超声波发射装置；5-信号传输线；6-超声波数据分析仪；6-1-超声波发射探头；6-2-超声波接收探头；6-3-超声波信号显示器；7-胶囊封孔器；A-胶囊封孔段；B-耦合剂充填段；C-剩余钻孔段。
具体实施方式
下面将对本发明作进一步说明。
如图1所示，本发明具体步骤为：
A、从岩巷1顶板向待测厚煤层2钻设多组穿层钻孔组3，且各穿层钻孔组3的钻孔角度不同，每个穿层钻孔组3均由9个穿层钻孔组成，且同组内9个穿层钻孔的钻孔方向平行；如图2所示，9个穿层钻孔呈3×3矩形阵列布置，相邻钻孔之间的间距为1m；每个穿层钻孔直径为150mm；各穿层钻孔组3均穿过煤层厚度直至煤层顶部与岩体结合处，然后利用岩心取芯机对各穿层钻孔组处于煤层内不同深度的位置进行取芯，取芯时在各个穿层钻孔内每隔0.5m取一次煤芯，直至岩心取芯机达到煤层顶部与岩体结合处时，停止取芯；所述钻取煤芯的尺寸为直径25mm、高度为50～60mm的圆柱体和直径80mm～100mm的圆柱体，接着将获取的煤芯作好标记并送至样品养护室内进行定期养护；
B、利用KDVJ-400型矿用钻孔窥视仪对各穿层钻孔组取芯位置的钻孔壁依次进行可视化探测，且在取芯位置每移动0.1m，输出一次拍摄的钻孔壁图像，最后通过多次图像集成，获取整个取芯位置的钻孔壁图像，进而识别出钻孔壁上存在的层理、割理结构信息；
C、选择一组穿层钻孔组3，将超声波发射装置4放入该组内中心位置的穿层钻孔3-1直至钻孔最深处，并通过信号传输线5与钻孔外部控制装置连接，如图4所示，所述超声波发射装置4上设有8个超声波发射探头6-1，8个超声波发射探头6-1能同时释放同等频率的超声波信号；且各个超声波发射探头6-1在该穿层钻孔内分别一一对应朝向该组内8个其余穿层钻孔3-2，8个其余穿层钻孔3-2分别放入一个超声波接收探头6-2，各个超声波接收探头6-2分别达到各个穿层钻孔的最深处、且各个超声波接收探头6-2均朝向该组中心位置的穿层钻孔3-1，各个超声波发射探头6-1与其对应的超声波接收探头6-2之间的对向角度偏移量不超过5°；然后对该组内所有穿层钻孔均采用胶囊封孔器7封孔，使各个穿层钻孔内在靠近孔口的位置形成胶囊封孔段A，胶囊封孔段A的长度为穿层钻孔长度的三分之一；所述胶囊封孔器7上装有穿过胶囊封孔器7的进液管和回液管，进液管一端和回液管一端均设有阀门、且回液管另一端伸入至穿层钻孔的最深处；打开进液管和回液管的阀门，采用注入泵将耦合剂以2～3MPa压力通过进液管分别注入各个穿层钻孔内，直至回液管一端流出耦合剂时，关闭回液管的阀门，并持续保压一段时间后停止耦合剂注入，并将进液管的阀门关闭，如图3所示，使各个穿层钻孔内均分成耦合剂充填段B、胶囊封孔段A和剩余钻孔段C；
D、如图5所示，开启超声波发射装置4，使其各个超声波发射探头6-1同时向该组内其余穿层钻孔3-2发射频率相同的超声波，该组内其余穿层钻孔3-2内各个超声波接收探头6-2分别接收各自对应的超声波发射探头6-1发出的超声波信号，均传递给超声波数据分析仪6，超声波数据分析仪6对超声波信号衰减信息进行反演，并通过超声波信号显示器6-3显示接收的超声波信号，从而得出超声波发射装置4所处穿层钻孔与其余各个穿层钻孔间在该深度位置的煤体基质信息；接着将超声波发射装置4及各个超声波接收探头6-2在各自穿层钻孔内均向孔口方向移动0.5m，并重复上述超声波发射及接收过程，从而得出超声波发射装置4所处穿层钻孔与其余各个穿层钻孔间在当前深度位置的煤体基质信息；如此重复，直至超声波发射装置4及各个超声波接收探头6-2在各自穿层钻孔内到达胶囊封孔器7时，停止超声波监测过程，进而获得步骤A所取该穿层钻孔组3内各个不同深度位置煤芯周围的煤体基质信息；对其余穿层钻孔组依次重复步骤C和D，获得各穿层钻孔组3内各个不同深度位置煤芯周围的煤体基质信息；
E、将步骤A标记的各穿层钻孔组3煤芯分别采用低场核磁共振技术、应力孔渗测试方法及三向纵波测试方法对煤芯的孔隙结构、应力敏感性及渗透率各向异性进行实验室测试，获取各穿层钻孔组3不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据，具体为：
对钻取煤芯分别进行真空脱水脱气及完全饱水处理，利用低场核磁共振技术测试得到饱水状态下的T2弛豫分布曲线；对饱水煤芯进行离心处理，同样利用低场核磁共振技术测试获得离心状态下的T2弛豫分布曲线；通过饱水-离心联测方法处理T2弛豫幅值差后，获得自由水空间占比，通过T2弛豫幅值差与孔隙度之间的线性关系，得到饱水、离心两种状态下的总孔隙率和残余孔隙率，二者之差等于有效孔隙率，最后根据SDR模型计算得到不同深度煤芯的核磁渗透率；
将煤芯置于应力三向非均质性的夹持器腔体中，通过施加不同荷载，获得煤芯在三个不同方向上的应力、渗透率敏感性数据；
对煤芯三个不同正交方向上不同测点位置的纵波进行测试，通过计算声波各向异性系数、波速比值、波列图分布参数获得煤体结构各向异性特征；其中，声波各向异性系数：
波速比值：
最终得到各穿层钻孔组不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据；
F、将步骤B和D通过现场测试获得的数据信息与步骤E通过实验室测试获得数据结果进行综合分析，由于同一穿层钻孔组内不同深度煤芯周围的岩体之间因为断层、多尺度孔裂隙、流体等信息的存在会表现出不同的基础力学特征，在钻孔现场测试过程当中，会由于地层结构相同出现部分重叠数据和异常数据，为获得钻孔测试的有效数据，需要将重叠数据和异常数据进行剔除，因此将所有数据进行重叠数据筛选、异常数据剔除，将筛选后的数据随机分成训练组和测试组，通过BP优化算法对训练组进行训练后创建BP神经网络，接着利用创建的BP神经网络对测试组进行测试，最终经过测试后获得非均质特性数据模型，为后续获得不同深度位置的参数输入至该数据模型即能计算得到对应位置的非均质特性，从而确定煤体各向异性。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种用于确定煤体各向异性的三维监测综合方法，其特征在于，具体步骤为：
A、从岩巷顶板向待测厚煤层钻设多组穿层钻孔组，且各个穿层钻孔组的钻孔角度不同，每个穿层钻孔组均由多个穿层钻孔组成，且同组内各个穿层钻孔的钻孔方向平行；各穿层钻孔组均穿过煤层厚度直至煤层顶部与岩体结合处，然后利用岩心取芯机对各穿层钻孔组处于煤层内不同深度的位置进行取芯，接着将获取的煤芯作好标记并送至样品养护室内进行定期养护；
B、利用钻孔窥视仪对各穿层钻孔组取芯位置的钻孔壁依次进行可视化探测，且在取芯位置每移动0.1m，输出一次拍摄的钻孔壁图像，最后通过多次图像集成，获取整个取芯位置的钻孔壁图像，进而识别出钻孔壁上存在的层理、割理结构信息；
C、选择一组穿层钻孔组，将超声波发射装置放入该组内其中一个穿层钻孔直至钻孔最深处，所述超声波发射装置上设有多个超声波发射探头，超声波发射探头的数量与该组内其余穿层钻孔的数量相同、且各个超声波发射探头在该穿层钻孔内分别一一对应朝向该组内其余穿层钻孔，该组内其余穿层钻孔分别放入一个超声波接收探头，各个超声波接收探头分别达到各个穿层钻孔的最深处、且各个超声波接收探头均朝向放入超声波发射装置的穿层钻孔；然后对该组内所有穿层钻孔均采用胶囊封孔器封孔，使各个穿层钻孔内在靠近孔口的位置形成胶囊封孔段；所述胶囊封孔器上装有穿过胶囊封孔器的进液管和回液管，进液管一端和回液管一端均设有阀门、且回液管另一端伸入至穿层钻孔的最深处；打开进液管和回液管的阀门，采用注入泵将耦合剂以一定压力通过进液管分别注入各个穿层钻孔内，直至回液管一端流出耦合剂时，关闭回液管的阀门，并持续保压一段时间后停止耦合剂注入，并将进液管的阀门关闭，使各个穿层钻孔内均分成耦合剂充填段、胶囊封孔段和剩余钻孔段；
D、开启超声波发射装置，使其各个超声波发射探头同时向该组内其余穿层钻孔发射频率相同的超声波，该组内其余穿层钻孔内各个超声波接收探头分别接收各自对应的超声波发射探头发出的超声波信号，均传递给超声波数据分析仪，超声波数据分析仪对超声波信号衰减信息进行反演，从而得出超声波发射装置所处穿层钻孔与其余各个穿层钻孔间在该深度位置的煤体基质信息；接着将超声波发射装置及各个超声波接收探头在各自穿层钻孔内均向孔口方向移动相同距离，并重复上述超声波发射及接收过程，从而得出超声波发射装置所处穿层钻孔与其余各个穿层钻孔间在当前深度位置的煤体基质信息；如此重复，直至超声波发射装置及各个超声波接收探头在各自穿层钻孔内到达胶囊封孔器时，停止超声波监测过程，进而获得步骤A所取该穿层钻孔组内各个不同深度位置煤芯周围的煤体基质信息；对其余穿层钻孔组依次重复步骤C和D，获得各穿层钻孔组内各个不同深度位置煤芯周围的煤体基质信息；
E、将步骤A标记的各穿层钻孔组煤芯分别采用低场核磁共振技术、应力孔渗测试方法及三向纵波测试方法对煤芯的孔隙结构、应力敏感性及渗透率各向异性进行实验室测试，获取各穿层钻孔组不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据；
F、将步骤B和D通过现场测试获得的数据信息与步骤E通过实验室测试获得数据结果进行综合分析，将所有数据通过重叠数据筛选、异常数据剔除处理，将筛选后的数据随机分成训练组和测试组，通过BP优化算法对训练组进行训练后创建BP神经网络，接着利用创建的BP神经网络对测试组进行测试，最终经过测试后获得非均质特性数据模型，后续获得不同深度位置的参数输入至该数据模型即能计算得到对应位置的非均质特性，从而确定煤体各向异性。
2.根据权利要求1所述的用于确定煤体各向异性的三维监测综合方法，其特征在于，所述步骤A中每个穿层钻孔组均由9个穿层钻孔组成，且9个穿层钻孔呈3×3矩形阵列布置，相邻钻孔之间的间距为1m；每个穿层钻孔直径为150mm；取芯时在各个穿层钻孔内每隔0.5m取一次煤芯，直至岩心取芯机达到煤层顶部与岩体结合处时，停止取芯；所述钻取煤芯的尺寸为直径25mm、高度为50～60mm的圆柱体和直径80mm～100mm的圆柱体。
3.根据权利要求1所述用于确定煤体各向异性的三维监测综合方法，其特征在于，所述步骤C中胶囊封孔段的长度为穿层钻孔长度的三分之一；所述耦合剂的注入压力为2～3MPa。
4.根据权利要求1所述用于确定煤体各向异性的三维监测综合方法，其特征在于，所述超声波发射装置上设有8个超声波发射探头，8个超声波发射探头能同时释放同等频率的超声波信号；各个超声波发射探头与其对应的超声波接收探头之间的对向角度偏移量不超过5°，超声波发射装置放入各穿层钻孔组内中心位置的穿层钻孔内；所述步骤D中超声波发射装置及各个超声波接收探头在各自穿层钻孔内每次移动距离为0.5m。
5.根据权利要求1所述用于确定煤体各向异性的三维监测综合方法，其特征在于，所述步骤E的具体过程为：
对钻取煤芯分别进行真空脱水脱气及完全饱水处理，利用低场核磁共振技术测试得到饱水状态下的T2弛豫分布曲线；对饱水煤芯进行离心处理，同样利用低场核磁共振技术测试获得离心状态下的T2弛豫分布曲线；通过饱水-离心联测方法处理T2弛豫幅值差后，获得自由水空间占比，通过T2弛豫幅值差与孔隙度之间的线性关系，得到饱水、离心两种状态下的总孔隙率和残余孔隙率，二者之差等于有效孔隙率，最后根据SDR模型计算得到不同深度煤芯的核磁渗透率；
将煤芯置于应力三向非均质性的夹持器腔体中，通过施加不同荷载，获得煤芯在三个不同方向上的应力、渗透率敏感性数据；
对煤芯三个不同正交方向上不同测点位置的纵波进行测试，通过计算声波各向异性系数、波速比值、波列图分布参数获得煤体结构各向异性特征；其中，声波各向异性系数：
波速比值：
最终得到各穿层钻孔组不同深度位置煤芯的孔隙结构、应力敏感性及渗透率各向异性数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210064568.4A CN114487125B (zh) | 2022-01-20 | 2022-01-20 | 一种用于确定煤体各向异性的三维监测综合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210064568.4A CN114487125B (zh) | 2022-01-20 | 2022-01-20 | 一种用于确定煤体各向异性的三维监测综合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114487125A CN114487125A (zh) | 2022-05-13 |
| CN114487125B true CN114487125B (zh) | 2023-08-11 |
Family
ID=81471788
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210064568.4A Active CN114487125B (zh) | 2022-01-20 | 2022-01-20 | 一种用于确定煤体各向异性的三维监测综合方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114487125B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115110921B (zh) * | 2022-07-18 | 2023-07-14 | 中煤科工集团重庆研究院有限公司 | 一种用于渗透各向异性本煤层瓦斯保压驱替抽采方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019071755A1 (zh) * | 2017-10-12 | 2019-04-18 | 中国矿业大学 | 煤层突出危险性随钻测试方法及装置 |
| CN109716169A (zh) * | 2016-07-01 | 2019-05-03 | 澳大利亚Nmr服务有限公司 | 井下扩散系数测量 |
| CN110174463A (zh) * | 2018-10-09 | 2019-08-27 | 天地科技股份有限公司 | 一种工作面三向采动应力场的无损定量测试方法 |
| CN112598986A (zh) * | 2020-12-21 | 2021-04-02 | 中国石油大学(北京) | 一种碳酸盐岩储层微观模型及其制备方法和应用 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8760657B2 (en) * | 2001-04-11 | 2014-06-24 | Gas Sensing Technology Corp | In-situ detection and analysis of methane in coal bed methane formations with spectrometers |
-
2022
- 2022-01-20 CN CN202210064568.4A patent/CN114487125B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109716169A (zh) * | 2016-07-01 | 2019-05-03 | 澳大利亚Nmr服务有限公司 | 井下扩散系数测量 |
| WO2019071755A1 (zh) * | 2017-10-12 | 2019-04-18 | 中国矿业大学 | 煤层突出危险性随钻测试方法及装置 |
| CN110174463A (zh) * | 2018-10-09 | 2019-08-27 | 天地科技股份有限公司 | 一种工作面三向采动应力场的无损定量测试方法 |
| CN112598986A (zh) * | 2020-12-21 | 2021-04-02 | 中国石油大学(北京) | 一种碳酸盐岩储层微观模型及其制备方法和应用 |
Non-Patent Citations (1)
| Title |
|---|
| 单轴加载煤孔裂隙各向异性核磁共振特征;卢方超;张玉贵;江林华;;煤田地质与勘探(01);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114487125A (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10884084B2 (en) | Systems and methods for tri-axial NMR testing | |
| Chen et al. | Investigation of hydraulic fracture propagation using a post-peak control system coupled with acoustic emission | |
| CN103513280B (zh) | 一种微地震监测模拟系统 | |
| CN101581701A (zh) | 一种地表注浆效果检测的方法 | |
| CN104532886B (zh) | 一种灌注桩桩底沉渣和桩端地基检测装置和方法 | |
| CN204401664U (zh) | 一种灌注桩桩底沉渣和桩端地基检测装置 | |
| CN108776175B (zh) | 一种冻结壁平均抗压强度超声波检测方法 | |
| CN105863626A (zh) | 一种钻井液与泥页岩地层理化作用的评价方法 | |
| CN109239779B (zh) | 一种隧道围岩松动圈的测试方法及围岩损伤的分级方法 | |
| CN104374827A (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
| CN114487125B (zh) | 一种用于确定煤体各向异性的三维监测综合方法 | |
| EP3811069A1 (en) | Systems and methods for tri-axial nmr testing | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN106401576A (zh) | 一种复杂页岩地层地应力的测试方法 | |
| CN105738215A (zh) | 一种通过声发射和差应变联合测试地应力的新方法 | |
| CN106869843B (zh) | 一种油水井物理法增产增注效果评价模拟实验方法 | |
| Yuzhang et al. | Experimental investigation of hydraulic fracture propagation in acoustic monitoring inside a large-scale polyaxial test | |
| CN110905012B (zh) | 一种建筑物桩基检测方法 | |
| CN113899811B (zh) | 一种煤矿巷道岩体累积性损伤的声波法测试系统 | |
| CN106869904B (zh) | 一种利用钻机运行参数原位实时确定岩体损伤状态的方法 | |
| Bearce et al. | Characterization of simulated soilcrete column curing using acoustic tomography | |
| Zhao | Application of the Kaiser effect of acoustic emission to measure vertical stress in an underground mine | |
| CN111596377B (zh) | 一种高地应力软岩隧道松动圈的联合测试方法 | |
| Falcon-Suarez et al. | Experimental rig to improve the geophysical and geomechanical understanding of CO2 reservoirs | |
| Hill et al. | Techniques for determining subsurface stress direction and assessing hydraulic fracture azimuth |
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 | ||
| CB02 | Change of applicant information |
Address after: 467000 No. 21 miner Road, Henan, Pingdingshan Applicant after: PINGDINGSHAN TIANAN COAL MINING Co.,Ltd. Applicant after: CHINA University OF MINING AND TECHNOLOGY Applicant after: China Pingmei Shenma Holding Group Co.,Ltd. Address before: 467000 No. 21 miner Road, Henan, Pingdingshan Applicant before: PINGDINGSHAN TIANAN COAL MINING Co.,Ltd. Applicant before: CHINA University OF MINING AND TECHNOLOGY Applicant before: CHINA PINGMEI SHENMA ENERGY CHEMICAL GROUP Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |