CN114895372B - 地下采空区探测综合解释方法及系统 - Google Patents
地下采空区探测综合解释方法及系统 Download PDFInfo
- Publication number
- CN114895372B CN114895372B CN202210342332.2A CN202210342332A CN114895372B CN 114895372 B CN114895372 B CN 114895372B CN 202210342332 A CN202210342332 A CN 202210342332A CN 114895372 B CN114895372 B CN 114895372B
- Authority
- CN
- China
- Prior art keywords
- data
- goaf
- detection
- determining
- initial
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了地下采空区探测综合解释方法及系统，首先在地下空间采空区的探测区域内，获取探测区域的地球物理数据、钻孔数据以及地质条件，地球物理数据包括电性数据、磁性数据以及重力数据；然后基于地球物理数据以及钻孔数据，确定出采空区在所述探测区域内的初始分布层位以及采空区的物性特征；此后对电性数据以及磁性数据进行反演，得到探测区域的三维电性结构，基于磁性数据、重力数据以及初始分布层位，确定采空区在所述三维电性结构内的初始范围；最后基于钻孔数据，对初始范围进行校正，并基于初始范围的校正结果，确定采空区的分布范围。该方法能够大大减少钻探工作量，提高采空区探测准确度，缩短勘探周期，为后期治理提供准确数据。
Description
技术领域
本发明是关于地球物理勘探技术领域，特别是关于一种地下采空区探测综合解释方法及系统。
背景技术
电磁法是资源和能源探查的重要手段之一，它以地壳中各种类型的物质的导电性、极化率等性质的差异为基本依据，通过观测电磁场的时空分布，实现地下目标体的有效探测。
地下地质体与围岩密度通常情况下存在一定的差异，引起重力异常，在地质体尺寸足够大，存在一定的剩余质量，附近干扰水平较低、埋藏条件有利的时候，地质体引起的重力异常能够被仪器探测，数据经过处理分析得到地下地质体空间分布以及物性特征等信息。
含有磁性的物质产生的磁场与地球磁场叠加，产生的新磁场发生畸变。通过探测分析被勘探物质引起的磁异常，研究被勘探物质的空间分布。磁力探测成本较低、效率高、应用范围广，应用范围较广。
钻探是直接获取地下物质的方法手段，能够准确的获得钻孔点的地下物质信息。获得的样品空间位置最为精确，能够通过测试得到对应的电性、磁性、重力信息和组成成分。
现有的技术中，虽然电磁法、磁法和重力采用的探测数据均可以从地面和空中获取，但是单独一个参数解释的结果往往与实际情况有较大差异。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
本发明的目的在于提供一种地下采空区探测综合解释方法及系统，其能够大大减少钻探工作量，提高采空区探测准确度，缩短勘探周期，为后期治理提供准确数据。
为实现上述目的，本发明提供了一种地下采空区探测综合解释方法，包括：
在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
在本发明的一实施方式中，所述基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征，包括：
基于所述钻孔数据，对所述地球物理数据进行聚类分析，概化出所述探测区域内各地层的物性特征以及所述采空区的物性特征；
基于所述各地层的物性数据以及所述钻孔数据，对所述采空区的层位进行初始化，确定所述初始分布层位。
在本发明的一实施方式中，所述基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征，之前包括：
对所述地球物理数据中的干扰数据进行剔除。
在本发明的一实施方式中，所述对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，包括：
采用横向约束反演方法，对所述电性数据以及所述磁性数据进行抽道、叠加、滤波以及水平调处理，得到处理结果；
剔除所述处理结果中的噪声点，得到所述三维电性结构。
在本发明的一实施方式中，所述基于所述钻孔数据，对所述初始范围进行校正，包括：
计算所述钻孔数据与所述三维电性结构之间的差异，并基于所述差异，确定所述钻孔数据的影响因子；
基于所述影响因子，对所述电性数据以及所述磁性数据进行约束反演，得到新的三维电性结构，并基于新的三维电性结构，对所述初始范围进行校正。
在本发明的一实施方式中，所述基于所述初始范围的校正结果，确定所述采空区的分布范围，包括：
判断新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据是否相符；
若相符，则基于所述校正结果，确定所述采空区的分布范围；
否则，确定划分所述采空区的物性边界条件，并迭代执行对所述初始范围的校正过程，直至新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据相符。
在本发明的一实施方式中，所述对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，之前包括：
基于所述电性数据，设置所述探测区域内各地层的电阻率分布。
本发明还提供一种地下采空区探测综合解释系统，包括：
数据获取模块，用于在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
初始化模块，用于基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
反演模块，用于对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
探测解释模块，用于基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述的地下采空区探测综合解释方法。
本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述的地下采空区探测综合解释方法。
本发明还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述的地下采空区探测综合解释方法。
与现有技术相比，根据本发明的地下采空区探测综合解释方法及系统，首先在地下空间采空区的探测区域内，获取探测区域的地球物理数据、钻孔数据以及地质条件，地球物理数据包括电性数据、磁性数据以及重力数据；然后基于地球物理数据以及钻孔数据，确定出采空区在所述探测区域内的初始分布层位以及采空区的物性特征；此后对电性数据以及磁性数据进行反演，得到探测区域的三维电性结构，基于磁性数据、重力数据以及初始分布层位，确定采空区在所述三维电性结构内的初始范围；最后基于钻孔数据，对初始范围进行校正，并基于初始范围的校正结果，确定采空区的分布范围。该方法能够大大减少钻探工作量，提高采空区探测准确度，缩短勘探周期，为后期治理提供准确数据。
附图说明
图1是本发明一实施方式提供的地下采空区探测综合解释方法的流程示意图之一；
图2是本发明一实施方式提供的典型的干扰数据随时间变化示意图；
图3是本发明一实施方式提供的地下采空区探测综合解释方法的流程示意图之二；
图4是本发明一实施方式提供的没有进行约束反演时得到的平面分布范围示意图；
图5是本发明一实施方式提供的带有约束反演时得到的平面分布范围示意图；
图6是本发明一实施方式提供的最终得到的采空区的分布范围示意图；
图7是本发明一实施方式提供的地下采空区探测综合解释系统的结构示意图；
图8是本发明一实施方式提供的电子设备的结构示意图。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
图1为本发明实施例中提供的一种地下采空区探测综合解释方法的流程示意图，如图1所示，该方法包括：
S1，在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
S2，基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
S3，对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
S4，基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
具体地，本发明实施例中提供的地下采空区探测综合解释方法，其执行主体为地下采空区探测综合解释系统，该系统可以配置于服务器内，该服务器可以是本地服务器，也可以是云端服务器，本地服务器具体可以是计算机等，本发明实施例中对此不作具体限定。
首先执行步骤S1，在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据。该地球物理数据可以包括探测区域的电性数据、磁性数据以及重力数据。钻孔数据可以包括钻孔获取的分层数据以及物性数据。可以理解的是，该地球物理数据以及钻孔数据均可以通过探测区域的区域资料中分类得到，地质条件可以通过经验确定。
然后执行步骤S2，通过地球物理数据以及钻孔数据，确定采空区在探测区域内的初始分布层位以及采空区的物性特征。可以理解的是，在步骤S2执行之前，可以先对地球物理数据中的干扰数据进行剔除，如此可以保证步骤S2执行的效率和准确性。如图2所示，为典型的干扰数据随时间变化示意图，图2横坐标为时间Time(单位为s)，纵坐标为dB/dt，单位为V/(A*M^4)。此处，需要对图2中右端的6个数据点进行剔除，同时还要注意地面建筑物干扰形成的信息。
在步骤S2执行之前，还可以通过电性数据，设置探测区域内各地层的电阻率分布，其目的是使反演结果与实际情况更接近，反演最终得出的电阻率误差更小。
进而，在执行步骤S2时，可以通过地球物理数据以及钻孔数据，确定采空区在探测区域内的初始分布层位以及采空区的物性特征。
此处，可以先通过钻孔数据，对地球物理数据进行聚类分析，归纳出采空区与围岩的差异，概化出探测区域内各地层的物性特征以及采空区的物性特征，以作为后续反演和解释的初始条件。此处，聚类分析结果即为探测区域内各地层的物性特征以及采空区的物性特征。然后，通过各地层的物性数据以及钻孔数据，对采空区的层位进行初始化，确定采空区的初始分布层位，即初步设定采空区可能存在的层位。
可以理解的是，层位初始化需要赋予初始层位界限、层间隔大小、层内物质的物性特征和组成成分。初始层位界限是指初始分布层位的最高高度以及最低高度。组成成分需要根据聚类分析结果进行修改，以达到反演所需效果。例如，可以分析各类岩石的电阻率，确定电阻率取值1-10为粘土，电阻率取值7-30为细沙，可以根据区分出的岩石类型，修改组成成分的赋值大小。
然后执行步骤S3，可以通过电磁法，对电性数据以及磁性数据进行反演，得到探测区域的三维电性结构，该三维电性结构可以通过电阻率数据进行表征。通过该三维电性结构可以初步确定采空区的平面分布范围和深度。然后，利用磁性数据以及重力数据，并结合初始分布层位，可以确定采空区在三维电性结构内的初始范围。该初始范围可以是采空区的平面分布范围，该初始范围可以是一个大概的范围，通过该大概的范围可以确定采空区的精确的范围。
最后执行步骤S4，可以根据钻孔数据，对初始范围进行校正。对初始范围进行校正的过程，即利用钻孔获取的分层数据，分析整理物性数据，剔除干扰数据后，约束电磁法再次反演得到初始范围的校正结果的过程。通过初始范围的校正结果是否符合约束条件，实现对采空区的分布范围的确定。
可以理解的是，最后得到的采空区的分布范围是三维分布范围，该分布范围是采空区的精确的范围。约束条件可以包括地质条件、采空区的物性特征以及钻孔数据等。
本发明实施例中提供的地下采空区探测综合解释方法，首先在地下空间采空区的探测区域内，获取探测区域的地球物理数据、钻孔数据以及地质条件，地球物理数据包括电性数据、磁性数据以及重力数据；然后基于地球物理数据以及钻孔数据，确定出采空区在所述探测区域内的初始分布层位以及采空区的物性特征；此后对电性数据以及磁性数据进行反演，得到探测区域的三维电性结构，基于磁性数据、重力数据以及初始分布层位，确定采空区在所述三维电性结构内的初始范围；最后基于钻孔数据，对初始范围进行校正，并基于初始范围的校正结果，确定采空区的分布范围。该方法能够大大减少钻探工作量，提高采空区探测准确度，缩短勘探周期，为后期治理提供准确数据。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释方法，所述对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，包括：
采用横向约束反演方法，对所述电性数据以及所述磁性数据进行抽道、叠加、滤波以及水平调处理，得到处理结果；
剔除所述处理结果中的噪声点，得到所述三维电性结构。
具体地，本发明实施例中，在对电性数据以及磁性数据进行反演，得到反演结果时，反演动作包括电性数据以及磁性数据抽道、叠加、滤波、水平调等处理步骤，得到处理结果，并剔除处理结果中能够检测到的噪声点，即可以得到探测区域的三维电性结构。本发明实施例中，可以采用横向约束反演方法实现反演动作。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释方法，所述基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围，包括：
具体地，本发明实施例中，可以利用磁性数据、重力数据对初始分布层位的地层连续性、平面分布等信息进行判别，得出采空区在三维电性结构内的初始范围，该初始范围即可能范围。
进而，在得到探测区域的三维电性结构之后，可以导入聚类分析结果，对采空区在三维电性结构内的初始范围进行约束。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释方法，所述基于所述钻孔数据，对所述初始范围进行校正，包括：
计算所述钻孔数据与所述三维电性结构之间的差异，并基于所述差异，确定所述钻孔数据的影响因子；
基于所述影响因子，对所述电性数据以及所述磁性数据进行约束反演，得到新的三维电性结构，并基于新的三维电性结构，对所述初始范围进行校正。
具体地，本发明实施例中，在通过钻孔数据对初始范围进行校正时，可以先计算钻孔数据与三维电性结构之间的差异。该差异可以是钻孔位置与三维电性结构内的参考点的位置之间的距离，通过该差异可以确定钻孔数据的可信度。可以理解的是，钻孔数据的可信度与参考点距离成反比。例如，如果参考点距离钻孔小于50m，则三维电性结构的分层与钻孔数据一致，即钻孔数据的可信度为100％，而参考点距离钻孔50m-100m，则钻孔数据的可信度为75％，参考点距离钻孔100m-200m，则钻孔数据的可信度为50％。
然后，可以通过钻孔数据的可信度，确定钻孔数据的影响因子，该影响因子可以是校正参数，用以对电性数据以及磁性数据的反演行为进行约束，并得到新的三维电性结构。进而，通过新的三维电性结构，对初始范围进行校正。特别地，可以配置设定距离，当上述差异超过设定距离时，钻孔数据的可信度将低于预设值，此时需要人为参与设计钻孔的影响范围，并对影响范围内的电阻率进行修正，进而对采空区的初始范围进行校正。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释方法，所述基于所述初始范围的校正结果，确定所述采空区的分布范围，包括：
判断新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据是否相符；
若相符，则基于所述校正结果，确定所述采空区的分布范围；
否则，确定划分所述采空区的物性边界条件，并迭代执行对所述初始范围的校正过程，直至新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据相符。
具体地，本发明实施例中，在根据采空区的初始范围的校正结果，确定采空区的分布范围时，需要判断新的三维电性结构与地质条件、采空区的物性特征以及钻孔数据是否相符，判断的过程可以是利用地质条件初步判别电阻率异常区厚度、范围是否与地层产状一致，是否符合采空区的工作设计等。
若相符，则基于校正结果，确定采空区的分布范围。此处，可以直接将校正结果作为采空区的分布范围，该分布范围可以包括平面分布和深度分布。
如果新的三维电性结构与地质条件、采空区的物性特征以及钻孔数据不相符，则需要确定划分采空区的物性边界条件，以限制边界外扩，并迭代执行对初始范围的校正过程，直至新的三维电性结构与地质条件、采空区的物性特征以及钻孔数据相符。
可以理解的是，物性边界条件是一种边界约束条件，是指划分边界时需要判别的地质干扰，如岩体、断裂带等因素，同时地表干扰影响的采集数据质量也会造成部分边界模糊的情况，物性边界条件可以通过人为确定。
在上述实施例的基础上，如图3所示，本发明实施例中提供的地下采空区探测综合解释方法，包括：
收集资料；
物性资料聚类分析；
电磁数据反演，即对电性数据以及磁性数据进行反演；
结合钻孔、物性约束，进行约束反演；
判断约束反演结果是否符合地质情况，如果符合则可以通过约束反演得到的三维电性结构对采空区的初始范围进行校正，得到采空区分布图；
如果约束反演结果不符合地质情况则重复进行约束反演，直至约束反演结果满足地质情况。
图4为没有进行约束反演时得到的平面分布范围示意图，图5为带有约束反演时得到的平面分布范围示意图，图6为最终得到的采空区的分布范围示意图。
如图7所示，在上述实施例的基础上，本发明实施例中提供了一种地下采空区探测综合解释系统，包括：
数据获取模块71，用于在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
初始化模块72，用于基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
反演模块73，用于对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
探测解释模块74，用于基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，所述初始化模块，用于：
基于所述钻孔数据，对所述地球物理数据进行聚类分析，概化出所述探测区域内各地层的物性特征以及所述采空区的物性特征；
基于所述各地层的物性数据以及所述钻孔数据，对所述采空区的层位进行初始化，确定所述初始分布层位。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，还包括数据处理模块，用于：
对所述地球物理数据中的干扰数据进行剔除。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，所述反演模块，用于：
采用横向约束反演方法，对所述电性数据以及所述磁性数据进行抽道、叠加、滤波以及水平调处理，得到处理结果；
剔除所述处理结果中的噪声点，得到所述三维电性结构。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，所述探测解释模块，用于：
计算所述钻孔数据与所述三维电性结构之间的差异，并基于所述差异，确定所述钻孔数据的影响因子；
基于所述影响因子，对所述电性数据以及所述磁性数据进行约束反演，得到新的三维电性结构，并基于新的三维电性结构，对所述初始范围进行校正。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，所述探测解释模块，还用于：
判断新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据是否相符；
若相符，则基于所述校正结果，确定所述采空区的分布范围；
否则，确定划分所述采空区的物性边界条件，并迭代执行对所述初始范围的校正过程，直至新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据相符。
在上述实施例的基础上，本发明实施例中提供的地下采空区探测综合解释系统，还包括设置模块，用于：
基于所述电性数据，设置所述探测区域内各地层的电阻率分布。
具体地，本发明实施例中提供的地下采空区探测综合解释系统中各模块的作用与上述方法类实施例中各步骤的操作流程是一一对应的，实现的效果也是一致的，具体参见上述实施例，本发明实施例中对此不再赘述。
图8示例了一种电子设备的实体结构示意图，如图8所示，该电子设备可以包括：处理器(processor)810、通信接口(Communications Interface)820、存储器(memory)830和通信总线840，其中，处理器810，通信接口820，存储器830通过通信总线840完成相互间的通信。处理器810可以调用存储器830中的逻辑指令，以执行上述各实施例中提供的地下采空区探测综合解释方法，该方法包括：在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
此外，上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各实施例中提供的地下采空区探测综合解释方法，该方法包括：在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例中提供的地下采空区探测综合解释方法，该方法包括：在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；对所述电性数据以及所述磁性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式，并且很显然，根据上述教导，可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用，从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。
Claims (9)
1.一种地下采空区探测综合解释方法，其特征在于，包括：
在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
对所述电性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围；所述基于所述钻孔数据，对所述初始范围进行校正，包括：计算所述钻孔数据与所述三维电性结构之间的差异，并基于所述差异，确定所述钻孔数据的影响因子；基于所述影响因子，对所述电性数据进行约束反演，得到新的三维电性结构，并基于新的三维电性结构，对所述初始范围进行校正。
2.如权利要求1所述的地下采空区探测综合解释方法，其特征在于，所述基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征，包括：
基于所述钻孔数据，对所述地球物理数据进行聚类分析，概化出所述探测区域内各地层的物性特征以及所述采空区的物性特征；
基于所述各地层的物性数据以及所述钻孔数据，对所述采空区的层位进行初始化，确定所述初始分布层位。
3.如权利要求1所述的地下采空区探测综合解释方法，其特征在于，所述基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征，之前包括：
对所述地球物理数据中的干扰数据进行剔除。
4.如权利要求1-3中任一项所述的地下采空区探测综合解释方法，其特征在于，所述对所述电性数据进行反演，得到所述探测区域的三维电性结构，包括：
采用横向约束反演方法，对所述电性数据进行抽道、叠加、滤波以及水平调处理，得到处理结果；
剔除所述处理结果中的噪声点，得到所述三维电性结构。
5.如权利要求1所述的地下采空区探测综合解释方法，其特征在于，所述基于所述初始范围的校正结果，确定所述采空区的分布范围，包括：
判断新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据是否相符；
若相符，则基于所述校正结果，确定所述采空区的分布范围；
否则，确定划分所述采空区的物性边界条件，并迭代执行对所述初始范围的校正过程，直至新的三维电性结构与所述地质条件、所述采空区的物性特征以及所述钻孔数据相符。
6.如权利要求1-3中任一项所述的地下采空区探测综合解释方法，其特征在于，所述对所述电性数据进行反演，得到所述探测区域的三维电性结构，之前包括：
基于所述电性数据，设置所述探测区域内各地层的电阻率分布。
7.一种地下采空区探测综合解释系统，其特征在于，包括：
数据获取模块，用于在地下空间采空区的探测区域内，获取所述探测区域的地球物理数据、钻孔数据以及地质条件，所述地球物理数据包括电性数据、磁性数据以及重力数据；
初始化模块，用于基于所述地球物理数据以及所述钻孔数据，确定所述采空区在所述探测区域内的初始分布层位以及所述采空区的物性特征；
反演模块，用于对所述电性数据进行反演，得到所述探测区域的三维电性结构，基于所述磁性数据、所述重力数据以及所述初始分布层位，确定所述采空区在所述三维电性结构内的初始范围；
探测解释模块，用于基于所述钻孔数据，对所述初始范围进行校正，并基于所述初始范围的校正结果，确定所述采空区的分布范围；所述基于所述钻孔数据，对所述初始范围进行校正，包括：计算所述钻孔数据与所述三维电性结构之间的差异，并基于所述差异，确定所述钻孔数据的影响因子；基于所述影响因子，对所述电性数据进行约束反演，得到新的三维电性结构，并基于新的三维电性结构，对所述初始范围进行校正。
8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述的地下采空区探测综合解释方法。
9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的地下采空区探测综合解释方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210342332.2A CN114895372B (zh) | 2022-04-02 | 2022-04-02 | 地下采空区探测综合解释方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210342332.2A CN114895372B (zh) | 2022-04-02 | 2022-04-02 | 地下采空区探测综合解释方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114895372A CN114895372A (zh) | 2022-08-12 |
| CN114895372B true CN114895372B (zh) | 2023-09-19 |
Family
ID=82714483
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210342332.2A Active CN114895372B (zh) | 2022-04-02 | 2022-04-02 | 地下采空区探测综合解释方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114895372B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115993103B (zh) * | 2023-03-22 | 2023-06-02 | 中国矿业大学(北京) | 一种采空区体积确定方法和装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111007578A (zh) * | 2019-12-06 | 2020-04-14 | 山东滨莱高速公路有限公司 | 一种路基岩溶及采空区的综合探测方法及系统 |
| AU2020102336A4 (en) * | 2020-09-18 | 2020-10-29 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method, system and device for identifying water-filled goaf area in coal mine |
-
2022
- 2022-04-02 CN CN202210342332.2A patent/CN114895372B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111007578A (zh) * | 2019-12-06 | 2020-04-14 | 山东滨莱高速公路有限公司 | 一种路基岩溶及采空区的综合探测方法及系统 |
| AU2020102336A4 (en) * | 2020-09-18 | 2020-10-29 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method, system and device for identifying water-filled goaf area in coal mine |
Non-Patent Citations (7)
| Title |
|---|
| 基于多源地学信息约束的繁昌盆地重磁联合反演及三维地质建模研究;丁文祥;中国优秀硕士学位论文全文数据库信息科技辑(第1期);摘要、第2-6章 * |
| 多参数综合物探方法在煤矿采空区勘探中的应用研究;罗霄;许智海;李正胜;李文;姜鹏;;煤炭工程;第52卷(第02期);全文 * |
| 综合物探反演解释方法在南京横溪-小丹阳示范区中的应用;黄宁, 刘国辉, 刘一;地质调查与研究;第39卷(第2期);第146-150页 * |
| 耦合地震约束信息的重磁电联合反演方法及其应用;杨敏, 于鹏, 朱光有, 张罗磊, 闫磊, 赵崇进, 马德波, 陈志勇;天然气地球科学;第33卷(第1期);全文 * |
| 重、磁、电联合反演在银额盆地定量解释中的应用;刘建利, 李西周, 张泉;物探与化探;第37卷(第5期);全文 * |
| 重磁电联合反演方法在天津基岩构造研究中的应用;马为, 李世斌, 徐新学, 郑军, 赵洪鹏, 李华强;地质调查与研究;第37卷(第3期);全文 * |
| 阳泉矿区煤田采空区综合地球物理探测模式应用;王森;王春辉;查恩来;翁爱华;;世界地质;第39卷(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114895372A (zh) | 2022-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8639443B2 (en) | Microseismic event monitoring technical field | |
| US11567232B2 (en) | Method for providing a calibrated rock-physics model of a subsoil | |
| CN114895372B (zh) | 地下采空区探测综合解释方法及系统 | |
| WO2016105765A1 (en) | Determining p-wave azimuthal anisotropy from walkaround vsp with offset dependent slowness corrections | |
| WO2017112033A1 (en) | A methodology for building realistic numerical forward stratigraphic models in data sparse environment | |
| CN111221048A (zh) | 基于跨孔电阻率ct多尺度反演的孤石边界识别与成像方法 | |
| WO2018067120A1 (en) | Geostatistical analysis of microseismic data in fracture modeling | |
| CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
| CN109630089B (zh) | 水平井地质构造识别方法及装置 | |
| CN111077577B (zh) | 井地联合储层描述方法及装置 | |
| CN106226217B (zh) | 一种确定地下裂缝张开度的方法、装置及其应用 | |
| CN112485841A (zh) | 深层地层岩性识别方法及装置 | |
| EP2929375B1 (en) | Method of using a data weighted electromagnetic source for time-lapse monitoring of a reservoir production or hydraulic fracturing | |
| CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
| CN115144929A (zh) | 油气藏储量确定方法、装置、存储介质及电子设备 | |
| AU2014209708B2 (en) | Method for determining fracture plane orientation using passive seismic signals | |
| CN112711076B (zh) | 提取石油钻井中泥浆侵入地层深度的方法和装置 | |
| Schrader et al. | Determining hydraulic parameters of a karst aquifer using unique historical data from large-scale dewatering by deep level mining–a case study from South Africa | |
| CN112485844B (zh) | 火山岩岩性岩相分布预测方法和预测装置 | |
| US20180306939A1 (en) | Correcting Biases In Microseismic-Event Data | |
| CN113311426A (zh) | 一种用于确定滑坡发生概率的方法及系统 | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
| CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
| US20240319392A1 (en) | Full waveform location technology and method to detect and locate microseismic events and characterise their moment tensor | |
| Doğru et al. | Comparison of EGM08 Bouguer Gravity with Ground Survey: A Case Study in Western Anatolian Region, Turkey |
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 |