CN117991348A - Space self-adaptive unstructured element resistivity tomography method and system - Google Patents

Space self-adaptive unstructured element resistivity tomography method and system Download PDF

Info

Publication number
CN117991348A
CN117991348A CN202410126376.0A CN202410126376A CN117991348A CN 117991348 A CN117991348 A CN 117991348A CN 202410126376 A CN202410126376 A CN 202410126376A CN 117991348 A CN117991348 A CN 117991348A
Authority
CN
China
Prior art keywords
resistivity
acquisition
depth
structural
inversion
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
CN202410126376.0A
Other languages
Chinese (zh)
Other versions
CN117991348B (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.)
Guangdong Heli Engineering Investigation Institute
Original Assignee
Guangdong Heli Engineering Investigation Institute
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 Guangdong Heli Engineering Investigation Institute filed Critical Guangdong Heli Engineering Investigation Institute
Priority to CN202410126376.0A priority Critical patent/CN117991348B/en
Publication of CN117991348A publication Critical patent/CN117991348A/en
Application granted granted Critical
Publication of CN117991348B publication Critical patent/CN117991348B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking

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

A space self-adaptive non-structural element resistivity tomography method and a system thereof relate to the technical field of geological exploration. The method comprises the following steps: acquiring resistivity and acquisition depth of a geological condition of a complex structure in a plurality of acquisition layers, and generating a first non-structural element corresponding to each acquisition layer according to each resistivity and each acquisition depth; performing inversion operation on each first unstructured element to obtain a first resistivity of each first unstructured element, and adjusting the first unstructured elements according to the first resistivity to obtain a second unstructured element; performing inversion operation on each second unstructured element to obtain a second resistivity of each second unstructured element; and analogically, finally generating the apparent resistivity layer image of the geological condition of the complex structure until the iteration is finished. The effect of improving the accuracy of resistivity tomography is achieved.

Description

空间自适应非结构元电阻率层析成像方法及系统Spatial adaptive non-structured element resistivity tomography method and system

技术领域Technical Field

本申请涉及地质勘探技术领域,具体涉及一种空间自适应非结构元电阻率层析成像方法及系统。The present application relates to the field of geological exploration technology, and in particular to a spatially adaptive non-structured element resistivity tomography method and system.

背景技术Background technique

随着地下资源勘查和环境监测领域对更精确、更高分辨率的地下模型的需求持续升高,复杂结构地质探测已经成为当前一个重要且具有挑战性的问题。在复杂地质环境中,介质属性表现为各项异性,我们通过各种工具和方法来理解环境、收集数据和处理信息,以便将这种环境变得更稳定、可认知和可描述。电阻率层析成像(ERT)是电法勘探的一种技术,它对具有不同电阻率差异的地电体响应敏感。由于这一特性,ERT技术一直以来都被广泛应用于水文地质、岩土工程、矿产资源、环境工程和考古的探测中。然而,考虑到非结构化环境的复杂性和可变性,我们需要进一步改进和优化当前的ERT技术,以更好地适应对复杂结构地质条件探测需求As the demand for more accurate and higher-resolution underground models continues to increase in the fields of underground resource exploration and environmental monitoring, complex structural geological detection has become an important and challenging issue today. In complex geological environments, the properties of the medium are anisotropic. We use various tools and methods to understand the environment, collect data, and process information in order to make this environment more stable, recognizable, and describable. Electrical resistivity tomography (ERT) is a technique for electrical exploration that is sensitive to geoelectric bodies with different resistivity differences. Due to this characteristic, ERT technology has been widely used in the detection of hydrogeology, geotechnical engineering, mineral resources, environmental engineering, and archaeology. However, considering the complexity and variability of unstructured environments, we need to further improve and optimize the current ERT technology to better meet the needs of detecting complex structural geological conditions.

目前,现有电阻率层析成像方法通常需要假设地下电阻率分布状况,并根据地下电阻率分布状况进行反演计算,最终将反演结果转换成像,现有技术多采用规则化非结构元单元成像,但是在实际应用中,地下电阻率分布可能会在不同的深度存在突变或者不连续的情况,规则化非结构元对于突变界面或曲面形态较难刻画,因此,现有技术得到的电阻率分布结果与地质体结构形态存在较大差异,导致电阻率层图像不准确。At present, the existing resistivity tomography methods usually need to assume the underground resistivity distribution, and perform inversion calculations based on the underground resistivity distribution, and finally convert the inversion results into images. The existing technologies mostly use regularized non-structured element unit imaging, but in actual applications, the underground resistivity distribution may be abrupt or discontinuous at different depths. Regularized non-structured elements are difficult to characterize abrupt interfaces or curved surface morphology. Therefore, the resistivity distribution results obtained by the existing technology are greatly different from the structural morphology of the geological body, resulting in inaccurate resistivity layer images.

发明内容Summary of the invention

本申请提供了一种空间自适应非结构元电阻率层析成像方法及系统,具有提高电阻率层析成像准确性的效果。The present application provides a spatially adaptive non-structured element resistivity tomography method and system, which has the effect of improving the accuracy of resistivity tomography.

第一方面,本申请提供了一种空间自适应非结构元电阻率层析成像方法,包括:In a first aspect, the present application provides a spatially adaptive non-structured element resistivity tomography method, comprising:

获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各所述电阻率和各所述采集深度,生成各所述采集层对应的第一非结构元;Obtaining resistivity and acquisition depths of multiple acquisition layers under complex structural geological conditions, and generating first non-structural elements corresponding to each acquisition layer according to each resistivity and each acquisition depth;

对各所述第一非结构元进行反演操作,得到各所述第一非结构元的第一电阻率,并根据所述第一电阻率调整所述第一非结构元,得到第二非结构元;Performing an inversion operation on each of the first non-structure elements to obtain a first resistivity of each of the first non-structure elements, and adjusting the first non-structure element according to the first resistivity to obtain a second non-structure element;

对各所述第二非结构元进行反演操作,得到各所述第二非结构元的第二电阻率;Performing an inversion operation on each of the second non-structural elements to obtain a second resistivity of each of the second non-structural elements;

重复所述反演操作,直到所述反演操作达到预设迭代终止条件,根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成所述复杂结构地质条件的视电阻率层图像。The inversion operation is repeated until the inversion operation reaches a preset iteration termination condition, and an apparent resistivity layer image of the complex structural geological condition is generated according to the second resistivity of each of the second non-structural elements and the acquisition depth.

通过采用上述技术方案,在非均质地层的多个采集层获取电阻率数据,结合采集深度信息建立初始非结构元,可以实现对复杂地质环境中电阻率体分布的模拟。然后利用迭代反演计算并自动优化非结构元的方式,可以不断提升电阻率计算的精度,使其结果更加准确地反映地层中的电阻率分布状态。电阻率计算结果可视化生成电阻率层析图像,与实际地质信息相结合,形成直观而富有针对性的解析表达,大大提高电阻率技术在复杂环境中的应用能力,提高了电阻率层析成像的准确性。By adopting the above technical solution, resistivity data is obtained in multiple acquisition layers of heterogeneous strata, and the initial non-structural element is established in combination with the acquisition depth information, so that the distribution of resistivity volume in complex geological environments can be simulated. Then, by using iterative inversion calculation and automatic optimization of non-structural elements, the accuracy of resistivity calculation can be continuously improved, so that the results can more accurately reflect the resistivity distribution state in the stratum. The resistivity calculation results are visualized to generate resistivity tomography images, which are combined with actual geological information to form an intuitive and targeted analytical expression, greatly improving the application ability of resistivity technology in complex environments and improving the accuracy of resistivity tomography.

可选的,根据各所述第一电阻率和各所述采集深度,匹配各所述采集层的深度分辨率和水平分辨率;根据各所述第一电阻率在各所述采集层中的分布位置,确定网格形状;根据所述深度分辨率、所述水平分辨率和所述网格形状,生成各所述采集层对应的所述第一非结构元。Optionally, the depth resolution and horizontal resolution of each acquisition layer are matched according to each first resistivity and each acquisition depth; the grid shape is determined according to the distribution position of each first resistivity in each acquisition layer; and the first non-structural element corresponding to each acquisition layer is generated according to the depth resolution, the horizontal resolution and the grid shape.

通过采用上述技术方案,匹配深度分辨率和水平分辨率,可以使网格节点设置符合实际测量条件,有利于提高计算精度。考虑电阻率分布位置,可以覆盖各电阻率体,避免重要信息被忽略。并进行网格形状优化,实现对复杂地层结构的合理离散化。该方案全面考虑了测量特征、地层分布和计算需要等多方面因素,使第一非结构元质量得到显著提升,为后续精准反演计算奠定了基础。这种精细化和智能化的网格生成技术,推进了电阻率测量与解释融合的创新,增强了方法的适应性。By adopting the above technical solution and matching the depth resolution and horizontal resolution, the grid node settings can be made to meet the actual measurement conditions, which is conducive to improving the calculation accuracy. Considering the resistivity distribution position, each resistivity body can be covered to avoid important information being ignored. And the grid shape is optimized to achieve reasonable discretization of complex stratigraphic structures. This solution comprehensively considers multiple factors such as measurement characteristics, stratigraphic distribution and calculation needs, which significantly improves the quality of the first non-structural element and lays the foundation for subsequent precise inversion calculations. This refined and intelligent grid generation technology has promoted the innovation of the fusion of resistivity measurement and interpretation and enhanced the adaptability of the method.

可选的,获取各所述采集层中水平方向的初始采集点位置和初始深度位置;根据各所述初始采集点位置和各所述分布位置中所述电阻率的水平位置,确定各所述电阻率的第一网格坐标;根据各所述初始深度位置和各所述分布位置中所述电阻率的深度位置,确定各所述电阻率的第二网格坐标;根据各所述第一网格坐标和各所述第二网格坐标,将各所述电阻率匹配至对应的所述第一非结构元。Optionally, the initial acquisition point position and the initial depth position in the horizontal direction in each acquisition layer are obtained; the first grid coordinates of each resistivity are determined according to the initial acquisition point position and the horizontal position of the resistivity in each distribution position; the second grid coordinates of each resistivity are determined according to the initial depth position and the depth position of the resistivity in each distribution position; and each resistivity is matched to the corresponding first non-structural element according to each first grid coordinate and each second grid coordinate.

通过采用上述技术方案,在实现电阻率数据到第一非结构元的映射时,不仅考虑了初始采集点的坐标信息,还综合利用了电阻率在采集层中的具体分布位置,通过匹配空间坐标与深度坐标,实现电阻率值与非结构元单元的精确对应。这种匹配方式充分考虑了电阻率的空间分布状态,避免了简单插值或整体映射可能带来的数据误差。电阻率值可以明确投影到反映其实际分布位置的非结构元单元上,保证了数据的准确性。该方案从数据层面提升了电阻率信息与计算非结构元的契合程度,使复杂环境中的电阻率分布特征可以无失真地映射到非结构元模型中,为后续精确反演计算提供了可靠的数据基础。By adopting the above technical solution, when realizing the mapping of resistivity data to the first non-structural element, not only the coordinate information of the initial acquisition point is considered, but also the specific distribution position of the resistivity in the acquisition layer is comprehensively utilized, and the resistivity value and the non-structural element unit are accurately corresponded by matching the spatial coordinates and the depth coordinates. This matching method fully considers the spatial distribution state of the resistivity and avoids the data errors that may be caused by simple interpolation or overall mapping. The resistivity value can be clearly projected onto the non-structural element unit that reflects its actual distribution position, ensuring the accuracy of the data. This solution improves the fit between the resistivity information and the calculated non-structural element from the data level, so that the resistivity distribution characteristics in complex environments can be mapped to the non-structural element model without distortion, providing a reliable data basis for subsequent precise inversion calculations.

可选的,根据各所述第一非结构元中的各所述电阻率和所述第一非结构元的空间范围,确定各所述第一非结构元的目标电阻率;根据所述目标电阻率和所述采集深度,确定地表电位测量值;若所述地表电位测量值未超出标准电位值范围,则将所述目标电阻率作为所述第一电阻率。Optionally, the target resistivity of each of the first non-structure elements is determined according to the resistivity in each of the first non-structure elements and the spatial range of the first non-structure element; the surface potential measurement value is determined according to the target resistivity and the acquisition depth; if the surface potential measurement value does not exceed the standard potential value range, the target resistivity is used as the first resistivity.

通过采用上述技术方案,在确定第一电阻率时,充分考虑了第一次的反演结果的准确性与合理性。一方面,设置目标电阻率作为约束条件,防止反演偏离实际情况;另一方面,通过计算地表电位并检验,找出在满足实测要求前提下的最佳电阻率分布情况。这种加入约束与检验的计算方式,可明显提升第一阶段反演的质量,找到更加接近实际的第一电阻率分布,为第二阶段反演提供可靠的输入。它避免了简单使用第一次的反演结果所带来的误差积累。By adopting the above technical solution, the accuracy and rationality of the first inversion result are fully considered when determining the first resistivity. On the one hand, the target resistivity is set as a constraint to prevent the inversion from deviating from the actual situation; on the other hand, the surface potential is calculated and tested to find the optimal resistivity distribution under the premise of meeting the actual measurement requirements. This calculation method that adds constraints and tests can significantly improve the quality of the first stage inversion, find a first resistivity distribution that is closer to the actual situation, and provide reliable input for the second stage inversion. It avoids the error accumulation caused by simply using the first inversion result.

可选的,根据所述第一电阻率和所述电阻率,确定目标差值;若所述目标差值大于预设差值,则将所述第一非结构元作为所述待调整非结构元;根据所述目标差值,确定修正系数,并根据所述修正系数和所述待调整非结构元,确定所述第二非结构元。Optionally, a target difference is determined based on the first resistivity and the resistivity; if the target difference is greater than a preset difference, the first non-structural element is used as the non-structural element to be adjusted; a correction coefficient is determined based on the target difference, and the second non-structural element is determined based on the correction coefficient and the non-structural element to be adjusted.

通过采用上述技术方案,建立了电阻率反演结果与非结构元调整之间的主动反馈机制。通过设置目标差值来评判第一阶段反演结果的准确性,引入修正系数进行定量调整,实现了第二非结构元的自动优化。这种反演结果主动反馈非结构元的调整方式,可以显著提升电阻率模拟计算的精度和自动化程度。它不再单纯依赖一次反演,而是加入连续迭代优化的理念,使复杂地层条件下电阻率体的刻画更加准确。该方案充分融合了电阻率反演与非结构元建模技术,推动电阻率测量解释向智能化、精细化方向发展,大幅提高电阻率技术适应复杂地质的能力。By adopting the above technical scheme, an active feedback mechanism between the resistivity inversion results and the adjustment of non-structured elements is established. By setting the target difference to judge the accuracy of the first-stage inversion results, the correction coefficient is introduced for quantitative adjustment, and the automatic optimization of the second non-structured element is realized. This inversion result actively feeds back the adjustment of non-structured elements, which can significantly improve the accuracy and automation of resistivity simulation calculations. It no longer relies solely on a single inversion, but incorporates the concept of continuous iterative optimization to make the characterization of resistivity bodies under complex formation conditions more accurate. This scheme fully integrates resistivity inversion and non-structured element modeling technology, promotes the development of resistivity measurement interpretation in the direction of intelligence and refinement, and greatly improves the ability of resistivity technology to adapt to complex geology.

可选的,根据所述目标差值,确定空间范围修正值;根据所述空间范围修正值和所述第一非结构元,生成所述第二非结构元。Optionally, a spatial range correction value is determined according to the target difference; and the second non-structure element is generated according to the spatial range correction value and the first non-structure element.

通过采用上述技术方案,建立了电阻率反演结果与非结构元空间范围优化之间的主动校正模型。通过分析第一阶段反演误差,确定空间范围的定量修正值,并根据修正值调整第二非结构元的覆盖范围,实现非结构元模型的自动优化。这种基于反演反馈进行的空间范围修正,可以使复杂地层中各电阻率体都被第二结构元有效覆盖,提升后续模拟计算的准确度。其避免了人工判断的主观影响,使结构元范围调整更加智能化。该方案属于电阻率反演与结构元建模深度融合的创新设计,推动电阻率层析成像技术向精细化与智能化方向发展,可显著提升电阻率在复杂环境应用的效果。By adopting the above technical solution, an active correction model between the resistivity inversion results and the optimization of the spatial range of non-structural elements was established. By analyzing the inversion error of the first stage, the quantitative correction value of the spatial range is determined, and the coverage range of the second non-structural element is adjusted according to the correction value to achieve automatic optimization of the non-structural element model. This spatial range correction based on inversion feedback can effectively cover each resistivity body in the complex formation with the second structural element, thereby improving the accuracy of subsequent simulation calculations. It avoids the subjective influence of manual judgment and makes the adjustment of the structural element range more intelligent. This solution is an innovative design that deeply integrates resistivity inversion and structural element modeling, which promotes the development of resistivity tomography technology in the direction of refinement and intelligence, and can significantly improve the effect of resistivity application in complex environments.

可选的,获取各所述采集层的地质信息;根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成各所述采集层的电阻率分布子图像;根据各所述电阻率分布子图像和各所述地质信息,生成所述复杂结构地质条件的视电阻率层图像。Optionally, geological information of each of the acquisition layers is obtained; a resistivity distribution sub-image of each of the acquisition layers is generated based on the second resistivity of each of the second non-structural elements and the acquisition depth; and an apparent resistivity layer image of the complex structural geological conditions is generated based on each of the resistivity distribution sub-images and the geological information.

通过采用上述技术方案,在电阻率计算的基础上,进一步获取各采集层的地质信息,并将计算得到的三维电阻率结果表达为与地质信息相结合的可视化电阻率层析图像。这种可视化表达形式,使复杂环境下的电阻率计算结果更加直观、形象。电阻率体的空间分布与各层地质信息结合在一起,图像既具备电阻率的精细特征,也反映出实际物理属性,提高了解释的效率和准确性。By adopting the above technical solution, on the basis of resistivity calculation, the geological information of each acquisition layer is further obtained, and the calculated three-dimensional resistivity results are expressed as a visual resistivity tomography image combined with geological information. This visual expression makes the resistivity calculation results in complex environments more intuitive and vivid. The spatial distribution of the resistivity body is combined with the geological information of each layer. The image not only has the fine characteristics of the resistivity, but also reflects the actual physical properties, which improves the efficiency and accuracy of the interpretation.

在本申请的第二方面提供了一种空间自适应非结构元电阻率层析成像系统。In a second aspect of the present application, a spatially adaptive non-structured element resistivity tomography system is provided.

数据获取模块,用于获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各所述电阻率和各所述采集深度,生成各所述采集层对应的第一非结构元;A data acquisition module, used to acquire the resistivity and acquisition depth of multiple acquisition layers under complex structural geological conditions, and generate a first non-structural element corresponding to each acquisition layer according to each resistivity and each acquisition depth;

第一反演模块,用于对各所述第一非结构元进行反演操作,得到各所述第一非结构元的第一电阻率,并根据所述第一电阻率调整所述第一非结构元,得到第二非结构元;A first inversion module is used to perform an inversion operation on each of the first non-structure elements to obtain a first resistivity of each of the first non-structure elements, and adjust the first non-structure element according to the first resistivity to obtain a second non-structure element;

第二反演模块,用于对各所述第二非结构元进行反演操作,得到各所述第二非结构元的第二电阻率;A second inversion module is used to perform an inversion operation on each of the second non-structure elements to obtain a second resistivity of each of the second non-structure elements;

图像生成模块,用于重复所述反演操作,直到所述反演操作达到预设迭代终止条件,根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成所述复杂结构地质条件的视电阻率层图像。An image generation module is used to repeat the inversion operation until the inversion operation reaches a preset iteration termination condition, and generate an apparent resistivity layer image of the complex structural geological condition according to the second resistivity of each second non-structural element and the acquisition depth.

一种空间自适应非结构元电阻率层析成像系统,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的程序,该程序能够被处理器加载执行时实现一种空间自适应非结构元电阻率层析成像方法。A spatially adaptive non-structured element resistivity tomography system comprises a memory, a processor and a program stored in the memory and executable on the processor. The program can realize a spatially adaptive non-structured element resistivity tomography method when loaded and executed by the processor.

在本申请的第三方面提供了一种计算机可读存储介质。In a third aspect of the present application, a computer-readable storage medium is provided.

一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时使所述处理器实现一种空间自适应非结构元电阻率层析成像方法。A computer-readable storage medium stores a computer program, which, when executed by a processor, enables the processor to implement a spatially adaptive non-structured element resistivity tomography method.

综上所述,本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:In summary, one or more technical solutions provided in the embodiments of the present application have at least the following technical effects or advantages:

1、本申请通过在非均质地层的多个采集层获取电阻率数据,结合采集深度信息建立初始非结构元,可以实现对复杂地质环境中电阻率体分布的模拟。然后利用迭代反演计算并自动优化非结构元的方式,可以不断提升电阻率计算的精度,使其结果更加准确地反映地层中的电阻率分布状态。电阻率计算结果可视化生成电阻率层析图像,与实际地质信息相结合,形成直观而富有针对性的解析表达,大大提高电阻率技术在复杂环境中的应用能力,提高了电阻率层析成像的准确性。1. This application obtains resistivity data from multiple acquisition layers in heterogeneous strata, and establishes initial non-structural elements in combination with acquisition depth information, so as to simulate the distribution of resistivity volumes in complex geological environments. Then, by using iterative inversion calculations and automatic optimization of non-structural elements, the accuracy of resistivity calculations can be continuously improved, so that the results can more accurately reflect the resistivity distribution state in the strata. The resistivity calculation results are visualized to generate resistivity tomography images, which are combined with actual geological information to form an intuitive and targeted analytical expression, greatly improving the application ability of resistivity technology in complex environments and improving the accuracy of resistivity tomography.

2、本申请通过在实现电阻率数据到第一非结构元的映射时,不仅考虑了初始采集点的坐标信息,还综合利用了电阻率在采集层中的具体分布位置,通过匹配空间坐标与深度坐标,实现电阻率值与非结构元单元的精确对应。这种匹配方式充分考虑了电阻率的空间分布状态,避免了简单插值或整体映射可能带来的数据误差。电阻率值可以明确投影到反映其实际分布位置的非结构元单元上,保证了数据的准确性。该方案从数据层面提升了电阻率信息与计算非结构元的契合程度,使复杂环境中的电阻率分布特征可以无失真地映射到非结构元模型中,为后续精确反演计算提供了可靠的数据基础。2. This application not only considers the coordinate information of the initial acquisition point when realizing the mapping of resistivity data to the first non-structure element, but also comprehensively utilizes the specific distribution position of the resistivity in the acquisition layer, and realizes the accurate correspondence between the resistivity value and the non-structure element unit by matching the spatial coordinates and the depth coordinates. This matching method fully considers the spatial distribution state of the resistivity and avoids the data errors that may be caused by simple interpolation or overall mapping. The resistivity value can be clearly projected onto the non-structure element unit that reflects its actual distribution position, ensuring the accuracy of the data. This scheme improves the fit between the resistivity information and the calculated non-structure element from the data level, so that the resistivity distribution characteristics in complex environments can be mapped to the non-structure element model without distortion, providing a reliable data foundation for subsequent precise inversion calculations.

3、本申请通过建立电阻率反演结果与非结构元空间范围优化之间的主动校正模型。通过分析第一阶段反演误差,确定空间范围的定量修正值,并根据修正值调整第二非结构元的覆盖范围,实现非结构元模型的自动优化。这种基于反演反馈进行的空间范围修正,可以使复杂地层中各电阻率体都被第二结构元有效覆盖,提升后续模拟计算的准确度。其避免了人工判断的主观影响,使非结构元范围调整更加智能化。该方案属于电阻率反演与结构元建模深度融合的创新设计,推动电阻率层析成像技术向精细化与智能化方向发展,可显著提升电阻率在复杂环境应用的效果。3. This application establishes an active correction model between the resistivity inversion results and the optimization of the spatial range of non-structural elements. By analyzing the inversion error of the first stage, the quantitative correction value of the spatial range is determined, and the coverage range of the second non-structural element is adjusted according to the correction value to achieve automatic optimization of the non-structural element model. This spatial range correction based on inversion feedback can effectively cover each resistivity body in the complex formation with the second structural element, thereby improving the accuracy of subsequent simulation calculations. It avoids the subjective influence of manual judgment and makes the adjustment of the non-structural element range more intelligent. This scheme is an innovative design that deeply integrates resistivity inversion and structural element modeling, which promotes the development of resistivity tomography technology in the direction of refinement and intelligence, and can significantly improve the effect of resistivity application in complex environments.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1是本申请实施例提供的一种空间自适应非结构元电阻率层析成像方法的流程示意图;FIG1 is a schematic flow chart of a spatially adaptive non-structured element resistivity tomography method provided in an embodiment of the present application;

图2是本申请实施例公开的一种空间自适应非结构元电阻率层析成像系统的结构示意图;FIG2 is a schematic diagram of the structure of a spatially adaptive non-structured element resistivity tomography system disclosed in an embodiment of the present application;

图3是本申请实施例的公开的一种电子设备的结构示意图。FIG. 3 is a schematic diagram of the structure of an electronic device disclosed in an embodiment of the present application.

附图标记说明:300、电子设备;301、处理器;302、通信总线;303、用户接口;304、网络接口;305、存储器。Description of reference numerals: 300, electronic device; 301, processor; 302, communication bus; 303, user interface; 304, network interface; 305, memory.

具体实施方式Detailed ways

为了使本领域的技术人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below in conjunction with the drawings in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments.

在本申请实施例的描述中,“例如”或者“举例来说”等词用于表示作例子、例证或说明。本申请实施例中被描述为“例如”或者“举例来说”的任何实施例或设计方案不应被解释为比其他实施例或设计方案更优选或更具优势。确切而言,使用“例如”或者“举例来说”等词旨在以具体方式呈现相关概念。In the description of the embodiments of the present application, words such as "for example" or "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "for example" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "for example" or "for example" is intended to present related concepts in a specific way.

在本申请实施例的描述中,术语“多个”的含义是指两个或两个以上。例如,多个系统是指两个或两个以上的系统,多个屏幕终端是指两个或两个以上的屏幕终端。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”,除非是以其他方式另外特别强调。In the description of the embodiments of the present application, the meaning of the term "multiple" refers to two or more. For example, multiple systems refer to two or more systems, and multiple screen terminals refer to two or more screen terminals. In addition, the terms "first" and "second" are used for descriptive purposes only and cannot be understood as indicating or implying relative importance or implicitly indicating the indicated technical features. Thus, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. The terms "include", "comprise", "have" and their variations all mean "including but not limited to", unless otherwise specifically emphasized.

为了便于理解本申请实施例提供的方法及系统,在介绍本申请实施例之前,先对本申请实施例的背景进行介绍。In order to facilitate understanding of the method and system provided by the embodiments of the present application, before introducing the embodiments of the present application, the background of the embodiments of the present application is first introduced.

目前,现有电阻率层析成像方法通常需要假设地下电阻率分布状况,并根据地下电阻率分布状况进行反演计算,最终将反演结果转换成像,现有技术多采用规则化非结构元单元成像,但是在实际应用中,地下电阻率分布可能会在不同的深度存在突变或者不连续的情况,规则化非结构元对于突变界面或曲面形态较难刻画,因此,现有技术得到的电阻率分布结果与地质体结构形态存在较大差异,导致电阻率层图像不准确。At present, the existing resistivity tomography methods usually need to assume the underground resistivity distribution, and perform inversion calculations based on the underground resistivity distribution, and finally convert the inversion results into images. The existing technologies mostly use regularized non-structured element unit imaging, but in actual applications, the underground resistivity distribution may be abrupt or discontinuous at different depths. Regularized non-structured elements are difficult to characterize abrupt interfaces or curved surface morphology. Therefore, the resistivity distribution results obtained by the existing technology are greatly different from the structural morphology of the geological body, resulting in inaccurate resistivity layer images.

本申请实施例公开了一种空间自适应非结构元电阻率层析成像方法,通过获取非结构化地质环境中各采集层的第一电阻率和采集深度,进行反演操作,然后根据反演操作的反演结果,调整各第一非结构元,并进行反演操作,从而根据反演操作的反演结果生成非结构化地质环境的电阻率层图像。主要用于解决地下电阻率分布可能会在不同的深度存在突变或者不连续的情况,往往通过简单的单次反演计算,得到的电阻率分布结果存在差异,导致电阻率成像不准确的问题。The embodiment of the present application discloses a spatially adaptive non-structured element resistivity tomography method, which obtains the first resistivity and acquisition depth of each acquisition layer in an unstructured geological environment, performs an inversion operation, and then adjusts each first non-structured element according to the inversion result of the inversion operation, and performs an inversion operation, thereby generating a resistivity layer image of the unstructured geological environment according to the inversion result of the inversion operation. It is mainly used to solve the problem that the underground resistivity distribution may be abrupt or discontinuous at different depths. Often, through a simple single inversion calculation, the resistivity distribution results obtained are different, resulting in inaccurate resistivity imaging.

经过上述背景内容相关介绍,本领域技术人员可以了解现有技术中存在的问题,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行详细的描述,描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。After the above background content is introduced, those skilled in the art can understand the problems existing in the prior art. The technical solutions in the embodiments of the present application will be described in detail below in conjunction with the drawings in the embodiments of the present application. The described embodiments are only part of the embodiments of the present application, not all of the embodiments.

参照图1,一种空间自适应非结构元电阻率层析成像方法,该方法包括S10至S40,具体包括以下步骤:Referring to FIG. 1 , a spatially adaptive non-structured element resistivity tomography method is provided, the method comprising S10 to S40, specifically comprising the following steps:

S10:获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各电阻率和各采集深度,生成各采集层对应的第一非结构元。S10: Obtain resistivity and acquisition depths of multiple acquisition layers under complex structural geological conditions, and generate first non-structural elements corresponding to each acquisition layer according to each resistivity and each acquisition depth.

其中,复杂结构地质条件指地质体中存在不规则变化的地层结构和地质条件的环境。其主要特征包括:地层结构复杂,存在不规则的褶曲、断裂等地质构造。岩性变化复杂,如砂岩、泥岩等岩性互层或不连续分布。地下水条件复杂,存在孔隙水、裂隙水等多种类型地下水。矿产、能源等资源分布不均匀,存在丰富区和贫乏区交替。地质体物性参数差异大,如电阻率、密度、声波速度等参量变化剧烈。地层中存在各向异性,地球物理参数随方向变化。地质尺度跨度大,存在多尺度的地质体特征。时间效应显著,地质条件随时间演化变化。Among them, complex structural geological conditions refer to the environment in which there are irregular changes in the stratigraphic structure and geological conditions in the geological body. Its main characteristics include: complex stratigraphic structure, irregular folds, fractures and other geological structures. Complex lithology changes, such as interlayers or discontinuous distribution of sandstone, mudstone and other lithologies. Complex groundwater conditions, with multiple types of groundwater such as pore water and fissure water. Mineral, energy and other resources are unevenly distributed, with alternations of rich and poor areas. The physical parameters of geological bodies vary greatly, such as resistivity, density, acoustic wave velocity and other parameters that change dramatically. There is anisotropy in the strata, and geophysical parameters change with direction. The geological scale span is large, and there are multi-scale geological body characteristics. The time effect is significant, and the geological conditions evolve and change with time.

具体的,通过电法勘探设备,在多个采集层按一定间距布设电极,传入不同频率的电流,测量各层之间的电位差,并结合电流强度参数计算出第一电阻率。记录各层电极的深度坐标作为采集深度。随后,根据测得的各层第一电阻率值的分布范围和位置信息,以及对应采集深度,可以确定各采集层的电阻率变化情况和深度分布特征。根据这些第一电阻率和采集深度的数据,就可以生成各层对应的自适应第一非结构元。非结构元生成原则是,依据常规网格剖面方法计算的灵敏度矩阵对成像区域修正剖分,修正方法采用非结构有限元法,电阻率梯度变化剧烈的位置网格细分,变化平缓的位置网格较粗。结合采集深度,在垂直方向上也进行适当细分。这样生成的结果称为第一非结构元,能够根据各采集层的电阻率分布情况进行自适应调整,为后续精确反演成像打下基础。Specifically, through electrical exploration equipment, electrodes are arranged at a certain interval in multiple acquisition layers, currents of different frequencies are transmitted, the potential difference between the layers is measured, and the first resistivity is calculated in combination with the current intensity parameter. The depth coordinates of the electrodes in each layer are recorded as the acquisition depth. Subsequently, according to the distribution range and position information of the measured first resistivity values of each layer, as well as the corresponding acquisition depth, the resistivity change and depth distribution characteristics of each acquisition layer can be determined. Based on these first resistivity and acquisition depth data, the adaptive first non-structural element corresponding to each layer can be generated. The principle of non-structural element generation is to correct and divide the imaging area according to the sensitivity matrix calculated by the conventional grid profile method. The correction method adopts the non-structural finite element method. The grid is subdivided at the location where the resistivity gradient changes violently, and the grid is coarser at the location where the change is gentle. Combined with the acquisition depth, appropriate subdivision is also performed in the vertical direction. The result generated in this way is called the first non-structural element, which can be adaptively adjusted according to the resistivity distribution of each acquisition layer, laying the foundation for subsequent accurate inversion imaging.

在上述实施例的基础上,生成第一非结构元的具体步骤还包括S11至S13:Based on the above embodiment, the specific steps of generating the first non-structural element further include S11 to S13:

S11:根据各电阻率和各采集深度,匹配各采集层的深度分辨率和水平分辨率。S11: Matching the depth resolution and horizontal resolution of each acquisition layer according to each resistivity and each acquisition depth.

其中,深度分辨率和水平分辨率是构建电阻率非结构元时的两个关键参数,代表非结构元在垂直方向和水平方向上的细分程度。Among them, depth resolution and horizontal resolution are two key parameters when constructing resistivity non-structural elements, representing the degree of subdivision of non-structural elements in the vertical and horizontal directions.

深度分辨率指非结构元在垂直方向上相邻层间的距离或非结构元层数。深度分辨率越高,垂直方向非结构元越密集,可以识别更精细的垂直结构变化。Depth resolution refers to the distance between adjacent layers of non-structural elements in the vertical direction or the number of non-structural element layers. The higher the depth resolution, the denser the non-structural elements in the vertical direction, and the more subtle vertical structural changes can be identified.

水平分辨率指非结构元在水平方向上相邻非结构元节点的距离或单元边长。水平分辨率越高,水平方向非结构元越密集,可以识别更精细的水平结构变化。Horizontal resolution refers to the distance between adjacent non-structural element nodes in the horizontal direction or the length of the unit side. The higher the horizontal resolution, the denser the non-structural elements in the horizontal direction, and the more subtle horizontal structural changes can be identified.

示例性地,根据已获得的各层第一电阻率数据,分析其水平和垂直方向的变化梯度。当水平变化梯度较大时,需要提高该层的水平分辨率,也就是在水平方向设计更密集的非结构元,以捕捉电阻率的变化细节。当垂直变化梯度较大时,需要提高垂直分辨率,在垂直方向设计更多的细分非结构元。根据各层的电阻率分布特点,确定其对应的优化深度分辨率和水平分辨率。进行匹配后,在电阻率变化剧烈的位置,通过提高分辨率,可以实现更精细的非结构元设计。这可以减少因非结构元过于粗略导致的反演误差,提高成像的精度。而在变化较为平缓的位置,则可以采用较低的分辨率,从而优化非结构元总数,减少计算量。匹配深度与水平分辨率,生成针对各采集层地质情况定制的自适应非结构元,是实现精确电阻率层析成像的重要环节,可以显著提升复杂环境中非结构元的划分准确性。Exemplarily, according to the first resistivity data of each layer that has been obtained, the horizontal and vertical gradients are analyzed. When the horizontal gradient is large, it is necessary to improve the horizontal resolution of the layer, that is, to design denser non-structural elements in the horizontal direction to capture the details of the resistivity change. When the vertical gradient is large, it is necessary to improve the vertical resolution and design more subdivided non-structural elements in the vertical direction. According to the resistivity distribution characteristics of each layer, the corresponding optimized depth resolution and horizontal resolution are determined. After matching, at locations where the resistivity changes dramatically, a more refined non-structural element design can be achieved by improving the resolution. This can reduce the inversion error caused by the non-structural elements being too rough and improve the accuracy of imaging. In locations where the changes are relatively gentle, a lower resolution can be used to optimize the total number of non-structural elements and reduce the amount of calculation. Matching depth and horizontal resolution to generate adaptive non-structural elements customized for the geological conditions of each acquisition layer is an important link in achieving accurate resistivity tomography, which can significantly improve the accuracy of the division of non-structural elements in complex environments.

S12:根据各电阻率在各采集层中的分布位置,确定网格形状。S12: Determine the grid shape according to the distribution position of each resistivity in each acquisition layer.

示例性地,根据已经获得的各层电阻率的分布范围和位置坐标,可以分析出电阻率分布的轮廓形态。例如电阻率存在环形分布、层状分布或不规则分布。根据这些分布形态,可以定义出不同形状的非结构元进行匹配。对于环形分布,可以生成圆形或者椭圆形的局部密集非结构元;对于层状分布,可以生成杆状或者带状的非结构元。综合考虑所有采集层的电阻率分布形态,设计出轮廓符合其分布形态的网格形状。这样生成的符合电阻率分布形态的不规则网格,可以最大限度地覆盖电阻率的变化范围,在其变化剧烈的位置进行密集逼近。这比规则非结构元可以更好地反映电阻率分布的真实形式,提高反演的精确性,也可降低计算量。For example, according to the distribution range and position coordinates of the resistivity of each layer that have been obtained, the contour morphology of the resistivity distribution can be analyzed. For example, the resistivity has an annular distribution, a layered distribution or an irregular distribution. According to these distribution forms, non-structural elements of different shapes can be defined for matching. For annular distribution, circular or elliptical local dense non-structural elements can be generated; for layered distribution, rod-shaped or strip-shaped non-structural elements can be generated. Taking into account the resistivity distribution morphology of all acquisition layers, a grid shape whose contour conforms to its distribution morphology is designed. The irregular grid generated in this way that conforms to the resistivity distribution morphology can cover the range of resistivity changes to the maximum extent and perform dense approximation at the location where it changes drastically. This can better reflect the true form of the resistivity distribution than regular non-structural elements, improve the accuracy of inversion, and reduce the amount of calculation.

S13:根据深度分辨率、水平分辨率和网格形状,生成各采集层对应的第一非结构元。S13: Generate a first non-structural element corresponding to each acquisition layer according to the depth resolution, the horizontal resolution and the grid shape.

示例性地,为综合考虑非结构化环境中各采集层电阻率分布的个性化特征,设计出针对各层地质条件定制的自适应非结构元。只有这样生成的非结构元,才能够有效提高复杂环境中电阻率数据的采集效率和反演精度。具体过程为根据前面步骤分析得到的各采集层的优化深度分辨率、水平分辨率以及匹配电阻率分布的网格形状。在第一采集层中,根据其深度分辨率和水平分辨率要求,生成轮廓形状匹配其电阻率分布的二维非结构元;在第二采集层中,重复该过程,生成符合其分辨率要求和分布形态的二维非结构元;以此类推,逐层设计。在垂直方向上,根据各层之间的边界关系,堆叠各层二维非结构元,构建成三维的第一非结构元。这样生成的非结构元,能够根据每个采集层所特有的地质条件变化规律进行定制,使反演模型达到对复杂环境的高精度表达。这种自适应非结构元的应用,将可以大幅提升非结构化环境中电阻率层析成像的效果。For example, in order to comprehensively consider the personalized characteristics of the resistivity distribution of each acquisition layer in an unstructured environment, an adaptive unstructured element customized for the geological conditions of each layer is designed. Only the unstructured element generated in this way can effectively improve the acquisition efficiency and inversion accuracy of resistivity data in a complex environment. The specific process is to optimize the depth resolution, horizontal resolution and grid shape matching the resistivity distribution of each acquisition layer obtained by the analysis of the previous steps. In the first acquisition layer, according to its depth resolution and horizontal resolution requirements, a two-dimensional unstructured element with a contour shape matching its resistivity distribution is generated; in the second acquisition layer, the process is repeated to generate a two-dimensional unstructured element that meets its resolution requirements and distribution morphology; and so on, layer by layer. In the vertical direction, according to the boundary relationship between the layers, the two-dimensional unstructured elements of each layer are stacked to construct a three-dimensional first unstructured element. The unstructured element generated in this way can be customized according to the change law of the geological conditions unique to each acquisition layer, so that the inversion model can achieve high-precision expression of the complex environment. The application of this adaptive unstructured element will greatly improve the effect of resistivity tomography in an unstructured environment.

在本申请一种可选实施例中,还存在将第一电阻率添至第一非结构元的过程,具体步骤包括S14至S16:In an optional embodiment of the present application, there is also a process of adding the first resistivity to the first non-structural element, and the specific steps include S14 to S16:

S14:获取各采集层中水平方向的初始采集点位置和初始深度位置;根据各初始采集点位置和各分布位置中电阻率的水平位置,确定各电阻率的第一网格坐标。S14: Acquire the initial acquisition point position and the initial depth position in the horizontal direction in each acquisition layer; determine the first grid coordinates of each resistivity according to the initial acquisition point position and the horizontal position of the resistivity in each distribution position.

示例性地,在进行电阻率测量时,需要精确记录各层初始采集点的平面坐标和深度坐标,以确定每个电阻率数据点的空间位置。在生成非结构元后,根据电阻率点的水平坐标,判断它在非结构元的哪个单元内;根据其深度坐标,确定电阻率位于非结构元的第几层。给每个电阻率点赋予确定的第一网格坐标。这样每个电阻率数据都会正确映射到非结构元的对应单元中,确保电阻率分布信息可以完整灌入非结构元,为后续精确反演建模奠定基础。这种处理可以消除数据配准误差。For example, when performing resistivity measurement, it is necessary to accurately record the plane coordinates and depth coordinates of the initial acquisition points of each layer to determine the spatial position of each resistivity data point. After the non-structure element is generated, the horizontal coordinates of the resistivity point are used to determine which unit of the non-structure element it is in; and the layer of the non-structure element in which the resistivity is located is determined based on its depth coordinates. Each resistivity point is assigned a certain first grid coordinate. In this way, each resistivity data will be correctly mapped to the corresponding unit of the non-structure element, ensuring that the resistivity distribution information can be fully injected into the non-structure element, laying the foundation for subsequent accurate inversion modeling. This processing can eliminate data registration errors.

S15:根据各初始深度位置和各分布位置中电阻率的深度位置,确定各电阻率的第二网格坐标。S15: Determine the second grid coordinates of each resistivity according to each initial depth position and the depth position of the resistivity in each distribution position.

示例性地,根据测量获得电阻率在各采集层的初始深度值。然后,检测每个电阻率实际出现的深度位置,与初始深度可能存在一定偏差。将此深度差值与非结构元的深度分辨率相结合,判断电阻率的具体深度位置相当于非结构元的第几层,即第二网格坐标的深度值。电阻率的深度信息可以通过第二坐标准确配准到非结构元的对应层位中。这避免了仅依据初始深度带来的映射偏差,提高了电阻率深度分布到非结构元的匹配精度,对保证后续精确反演至关重要。该方案克服了复杂环境中电阻率深度变化的不确定性,实现准确的三维映射,大大提高电阻率数据到非结构元的配准效果,是电阻率层析成像获取高精度结果的关键环节。Exemplarily, the initial depth value of the resistivity in each acquisition layer is obtained based on the measurement. Then, the actual depth position of each resistivity is detected, which may deviate from the initial depth. This depth difference is combined with the depth resolution of the non-structural element to determine which layer of the non-structural element the specific depth position of the resistivity is equivalent to, that is, the depth value of the second grid coordinate. The depth information of the resistivity can be accurately aligned to the corresponding layer of the non-structural element through the second coordinate. This avoids the mapping deviation caused by relying solely on the initial depth, improves the matching accuracy of the resistivity depth distribution to the non-structural element, and is crucial to ensure subsequent accurate inversion. This scheme overcomes the uncertainty of resistivity depth changes in complex environments, realizes accurate three-dimensional mapping, greatly improves the registration effect of resistivity data to non-structural elements, and is a key link in obtaining high-precision results from resistivity tomography.

S16:根据各第一网格坐标和各第二网格坐标,将各电阻率匹配至对应的第一非结构元。S16: Match each resistivity to a corresponding first non-structural element according to each first grid coordinate and each second grid coordinate.

示例性地,根据已测量的各电阻率点的第一网格坐标,判断其在非结构元的对应位置;然后根据第二网格坐标,确定其在非结构元的垂直层位。综合这两个坐标信息,就可直接匹配出每个第一电阻率对应非结构元的三维单元。可以进行自动批处理,迅速实现大规模电阻率到非结构元的匹配。这样,第一电阻率数据就可以无偏差地置入非结构元相应的单元。这避免了复杂环境下人工映射的误差,大大提高了电阻率信息在非结构元中的表示准确性,为后续精细反演建模奠定了基础。For example, according to the first grid coordinates of each measured resistivity point, its corresponding position in the non-structure element is determined; then according to the second grid coordinates, its vertical layer in the non-structure element is determined. Combining these two coordinate information, the three-dimensional unit of each first resistivity corresponding to the non-structure element can be directly matched. Automatic batch processing can be performed to quickly achieve large-scale resistivity to non-structure element matching. In this way, the first resistivity data can be placed in the corresponding unit of the non-structure element without deviation. This avoids the error of manual mapping in complex environments, greatly improves the representation accuracy of resistivity information in non-structure elements, and lays the foundation for subsequent fine inversion modeling.

S20:对各第一非结构元进行反演操作,得到各第一非结构元的第一电阻率,并根据第一电阻率调整第一非结构元,得到第二非结构元。S20: performing an inversion operation on each first non-structure element to obtain a first resistivity of each first non-structure element, and adjusting the first non-structure element according to the first resistivity to obtain a second non-structure element.

其中,反演操作是指在电阻率层析成像方法中,基于构建的初始非结构元,进行的第一次数值反演计算,以求解电阻率分布的过程。The inversion operation refers to the process of performing the first numerical inversion calculation based on the constructed initial non-structured element in the resistivity tomography method to solve the resistivity distribution.

具体的,将设计好的第一非结构元输入计算机,根据非结构元中各单元的第一电阻率值,通过数值计算的方式解出每个单元的第一电阻率。此过程需要运用到有限元法、有限差分法等数值法,进行复杂的非结构元计算。得到第一电阻率后,需要根据其分布情况,判断第一非结构元是否需要进一步优化。如果电阻率变化剧烈的区域非结构元不够密集,则需要局部调整非结构元,获得第二非结构元。这一调整是实现第二次反演精度的关键。通过第一反演和根据其结果调整非结构元,可以实现电阻率分布的初步解算,并进一步优化计算非结构元,为第二次反演建立优化模型,提高最终成像效果。Specifically, the designed first non-structure element is input into the computer, and the first resistivity of each unit in the non-structure element is solved by numerical calculation according to the first resistivity value of each unit in the non-structure element. This process requires the use of numerical methods such as the finite element method and the finite difference method to perform complex non-structure element calculations. After obtaining the first resistivity, it is necessary to determine whether the first non-structure element needs to be further optimized based on its distribution. If the non-structure elements in the area where the resistivity changes dramatically are not dense enough, it is necessary to locally adjust the non-structure elements to obtain the second non-structure element. This adjustment is the key to achieving the accuracy of the second inversion. Through the first inversion and adjusting the non-structure elements according to its results, the preliminary solution of the resistivity distribution can be achieved, and the non-structure elements can be further optimized and calculated to establish an optimization model for the second inversion and improve the final imaging effect.

在上述实施例的基础上,确定第一电阻率的具体步骤还包括S21至S22:Based on the above embodiment, the specific step of determining the first resistivity further includes S21 to S22:

S21:根据各第一非结构元中的各电阻率和第一非结构元的空间范围,确定各第一非结构元的目标电阻率。S21: Determine the target resistivity of each first non-structure element according to each resistivity in each first non-structure element and the spatial range of the first non-structure element.

示例性地,根据第一非结构元中各单元格的初始电阻率值,结合非结构元的空间范围,进行整体分析。以非结构元为单位,计算各单元的电阻率平均值作为目标电阻率。也可以根据场地的已知地质情况,确定电阻率目标范围。设置这些目标电阻率后,在进行第二次反演计算时,以目标值约束反演运算的收敛结果。这可以避免第一反演带来的误差累积扩大,使第二次反演结果更准确地接近实际电阻率分布。通过设置目标电阻率,可在一定程度上纠正第一反演偏差,将第二次反演导向期望结果,显著提高电阻率层析成像的精确性和可靠性。Exemplarily, an overall analysis is performed based on the initial resistivity values of each cell in the first non-structural element and the spatial range of the non-structural element. Taking the non-structural element as a unit, the resistivity average value of each unit is calculated as the target resistivity. The resistivity target range can also be determined based on the known geological conditions of the site. After setting these target resistivities, when performing the second inversion calculation, the convergence result of the inversion operation is constrained by the target value. This can avoid the accumulation and expansion of errors caused by the first inversion, and make the second inversion result more accurately close to the actual resistivity distribution. By setting the target resistivity, the deviation of the first inversion can be corrected to a certain extent, and the second inversion can be directed to the expected result, which significantly improves the accuracy and reliability of resistivity tomography.

S22:根据目标电阻率和采集深度,确定地表电位测量值;若地表电位测量值未超出标准电位值范围,则将目标电阻率作为第一电阻率。S22: Determine the surface potential measurement value according to the target resistivity and the acquisition depth; if the surface potential measurement value does not exceed the standard potential value range, use the target resistivity as the first resistivity.

示例性地,将目标电阻率分布带入预先建立的非结构元模型,根据采集深度情形,通过前述运算模式计算出理论上的地表电位分布。然后,将其与实际在对应测线位置测得的地表电位值进行对比。如果计算结果符合标准电位值的误差范围,则认为目标电阻率计算正确,可以作为第二电阻率结果。如果对比发现超出误差范围,则需要重新调整目标电阻率,重复运算,直至测量值验证通过。通过地表电位验证,可以有效检验第二电阻率结果的准确性,避免反演误差累积扩大,大大提高结果的可信度。这是电阻率层析成像方法实现高精度的重要保障。Exemplarily, the target resistivity distribution is brought into the pre-established unstructured element model, and the theoretical surface potential distribution is calculated through the aforementioned operation mode according to the acquisition depth. Then, it is compared with the actual surface potential value measured at the corresponding survey line position. If the calculated result meets the error range of the standard potential value, it is considered that the target resistivity calculation is correct and can be used as the second resistivity result. If the comparison finds that it exceeds the error range, it is necessary to readjust the target resistivity and repeat the operation until the measurement value is verified. Through surface potential verification, the accuracy of the second resistivity result can be effectively tested, the accumulation and expansion of inversion errors can be avoided, and the credibility of the results can be greatly improved. This is an important guarantee for the resistivity tomography method to achieve high precision.

在上述实施例的基础上,确定第二非结构元的具体步骤还包括S23至S24:Based on the above embodiment, the specific step of determining the second non-structural element further includes S23 to S24:

S23:根据第一电阻率和电阻率,确定目标差值;若目标差值大于预设差值,则将第一非结构元作为待调整非结构元。S23: determining a target difference according to the first resistivity and the resistivity; if the target difference is greater than a preset difference, taking the first non-structural element as a non-structural element to be adjusted.

示例性地,将第二电阻率数据与第一电阻率数据逐点进行比较,计算两者间的绝对误差或相对误差,即目标差值。如果差值普遍大于预设的允许误差范围,说明存在系统偏差,对非结构元模型进行调整,以提高计算精度。当目标差值超过阈值时,第一非结构元就会被标记为待调整非结构元。接下来,可以检查非结构元中电阻率变化剧烈的区域是否分辨率不足,从而针对性地提高这些局部区域的水平或垂直分辨率,加入细分的非结构元单元,进行局部调整,以消除先前的偏差。通过这种自动化的误差检验和非结构元优化机制,可以有效消除反演误差的累积,保证电阻率层析成像的准确性和可靠性。Exemplarily, the second resistivity data is compared point by point with the first resistivity data, and the absolute error or relative error between the two, that is, the target difference, is calculated. If the difference is generally greater than the preset allowable error range, it means that there is a systematic deviation, and the non-structural element model is adjusted to improve the calculation accuracy. When the target difference exceeds the threshold, the first non-structural element will be marked as a non-structural element to be adjusted. Next, it can be checked whether the area where the resistivity changes drastically in the non-structural element has insufficient resolution, so as to improve the horizontal or vertical resolution of these local areas in a targeted manner, add subdivided non-structural element units, and make local adjustments to eliminate previous deviations. Through this automated error detection and non-structural element optimization mechanism, the accumulation of inversion errors can be effectively eliminated, ensuring the accuracy and reliability of resistivity tomography.

S24:根据目标差值,确定修正系数,并根据修正系数和待调整非结构元,确定第二非结构元。S24: determining a correction coefficient according to the target difference, and determining a second non-structural element according to the correction coefficient and the non-structural element to be adjusted.

示例性地,根据目标差值的大小确定修正系数,即本次非结构元调整的幅度。在待调整的第一非结构元的基础上,根据修正系数来确定该如何调整非结构元参数,比如在电阻率变化剧烈的局部区域增加水平和垂直方向的节点密度,进行细化分区,或者扩大某方向的非结构元范围。经过这种针对性调整后,即可获得第二非结构元。新生成的第二非结构元可以更好地逼近地层的电阻率分布状况。在该优化非结构元的基础上进行第二次反演,可以大大提高计算结果的准确性。这种自动修正机制使得电阻率层析成像可不断迭代优化,从而实现高精度反演与成像。Exemplarily, the correction coefficient, that is, the amplitude of the non-structural element adjustment this time, is determined according to the size of the target difference. On the basis of the first non-structural element to be adjusted, the correction coefficient is used to determine how to adjust the non-structural element parameters, such as increasing the node density in the horizontal and vertical directions in the local area where the resistivity changes dramatically, performing refined partitioning, or expanding the range of non-structural elements in a certain direction. After this targeted adjustment, the second non-structural element can be obtained. The newly generated second non-structural element can better approximate the resistivity distribution of the formation. Performing a second inversion based on the optimized non-structural element can greatly improve the accuracy of the calculation results. This automatic correction mechanism enables resistivity tomography to be continuously iteratively optimized, thereby achieving high-precision inversion and imaging.

在上述实施例的基础上,确定修正系数的具体步骤还包括S241至S243:Based on the above embodiment, the specific steps of determining the correction coefficient further include S241 to S243:

S241:根据目标差值,确定空间范围修正值。S241: Determine the spatial range correction value according to the target difference.

示例性地,如果目标差值较大,说明存在非结构元范围不足以覆盖电阻率异常变化区域的情况。这时可以根据差值大小确定一个空间范围修正值,例如沿某一方向延展10米。然后基于第一非结构元的范围扩充这个修正值,重新确定第二非结构元的空间范围。如果目标差值较小,则可以适当缩小非结构元范围,减少计算量。经过这种根据误差调整非结构元范围的过程,可以缩小反演的空间边界误差,有助于提高电阻率分布结果的准确性。For example, if the target difference is large, it means that the range of the non-structure element is not enough to cover the area of abnormal resistivity change. At this time, a spatial range correction value can be determined according to the difference size, for example, extending 10 meters in a certain direction. Then, based on the range of the first non-structure element, this correction value is expanded to redefine the spatial range of the second non-structure element. If the target difference is small, the range of the non-structure element can be appropriately narrowed to reduce the amount of calculation. Through this process of adjusting the range of the non-structure element according to the error, the spatial boundary error of the inversion can be reduced, which helps to improve the accuracy of the resistivity distribution results.

S242:根据空间范围修正值和第一非结构元,生成第二非结构元。S242: Generate a second non-structure element according to the spatial range correction value and the first non-structure element.

示例性地,在确定了空间范围修正值后,以第一非结构元为基础,按照修正值的方向和大小进行范围扩充或缩小。例如,如果确定沿x轴正向扩充10米,则将第一非结构元在该方向上进行拓展,增设10米范围的非结构元单元。经过这种针对性调整后,就可以生成范围得到优化的第二非结构元。这样,新生成的第二非结构元可以更准确匹配电阻率体的空间分布范围。在该范围优化非结构元上进行电阻率反演,可以大大减少反演边界的截断误差,提高成像效果。通过自动调整非结构元范围,可实现模型的自适应优化,使复杂电阻率异常体的成像更加准确。Exemplarily, after determining the spatial range correction value, the range is expanded or reduced based on the first non-structural element according to the direction and size of the correction value. For example, if it is determined that the expansion is to be 10 meters along the positive direction of the x-axis, the first non-structural element is expanded in this direction, and a non-structural element unit with a range of 10 meters is added. After this targeted adjustment, a second non-structural element with an optimized range can be generated. In this way, the newly generated second non-structural element can more accurately match the spatial distribution range of the resistivity body. Performing resistivity inversion on the non-structural element optimized in this range can greatly reduce the truncation error of the inversion boundary and improve the imaging effect. By automatically adjusting the range of the non-structural element, the adaptive optimization of the model can be achieved, making the imaging of complex resistivity anomalies more accurate.

S30:对各第二非结构元进行反演操作,得到各第二非结构元的第二电阻率。S30: performing an inversion operation on each second non-structure element to obtain a second resistivity of each second non-structure element.

具体的,将经过调整优化的第二非结构元作为计算模型。采用有限元法或有限差分法等数值计算方法,进行复杂的三维电阻率反演运算。根据设定的目标电阻率约束条件,经迭代计算,即可获得每个非结构元单元的第三电阻率值。由于非结构元质量得到提升,计算结果可以真实反映地层中电阻率体的三维分布状况,达到很高的精度和分辨率。通过后处理,即可生成清晰的电阻率层析剖面成像结果。Specifically, the adjusted and optimized second non-structural element is used as the calculation model. Numerical calculation methods such as the finite element method or the finite difference method are used to perform complex three-dimensional resistivity inversion operations. According to the set target resistivity constraints, the third resistivity value of each non-structural element unit can be obtained through iterative calculation. Since the quality of non-structural elements is improved, the calculation results can truly reflect the three-dimensional distribution of resistivity bodies in the formation, achieving high accuracy and resolution. Through post-processing, clear resistivity tomography profile imaging results can be generated.

S40:重复反演操作,直到反演操作达到预设迭代终止条件,根据各第二非结构元的第二电阻率和采集深度,生成复杂结构地质条件的视电阻率层图像。S40: Repeat the inversion operation until the inversion operation reaches a preset iteration termination condition, and generate an apparent resistivity layer image of complex structural geological conditions according to the second resistivity and acquisition depth of each second non-structural element.

其中,复杂结构地质条件的视电阻率层图像是指在复杂、不规则的地质环境中,根据计算得到的三维电阻率分布,将其划分为不同电阻率区间的层,并通过可视化生成电阻率剖面图,以直观表示各个电阻率体的空间分布形式。该非结构化地质环境的电阻率层图像反映复杂地质条件下电阻率分布的大致状况,可快速识别电阻率异常体,是电阻率层析成像的关键结果表达形式。Among them, the apparent resistivity layer image of complex structural geological conditions refers to the division of the complex and irregular geological environment into layers of different resistivity intervals according to the calculated three-dimensional resistivity distribution, and the generation of resistivity profiles through visualization to intuitively represent the spatial distribution of each resistivity body. The resistivity layer image of the unstructured geological environment reflects the general situation of resistivity distribution under complex geological conditions, can quickly identify resistivity anomalies, and is a key result expression form of resistivity tomography.

具体的,进行该步骤的目的是将计算得到的三维电阻率分布结果,转换成清晰的电阻率剖面图,直观表示复杂地质条件下不同电阻率体的空间分布状况。Specifically, the purpose of this step is to convert the calculated three-dimensional resistivity distribution results into a clear resistivity profile diagram, which intuitively represents the spatial distribution of different resistivity bodies under complex geological conditions.

重复反演操作,并在每次反演操作后检测该第二电阻率,当该第二电阻率与标准电阻率的误差率在预设阈值内时,判断达到预设迭代终止条件。结束反演操作后,自动读取各个第二非结构元单元的第三电阻率数据,并根据采集深度,将具有相近电阻率的单元划分到对应的阻率区间内。根据非结构元单元的空间坐标,将不同阻率区间的单元渲染到二维剖面平面上,形成色阶表达的电阻率剖面图。这样就可以直观表达复杂三维电阻率异质体在二维剖面上的分布效果。不同颜色对应的是不同电阻率的值和范围,呈现出非结构化地质条件下电阻率体的形态。这种成像结果更加直观,便于后续地质解析。Repeat the inversion operation, and detect the second resistivity after each inversion operation. When the error rate between the second resistivity and the standard resistivity is within the preset threshold, it is judged that the preset iteration termination condition is reached. After the inversion operation is completed, the third resistivity data of each second non-structural element unit is automatically read, and the units with similar resistivity are divided into corresponding resistivity intervals according to the acquisition depth. According to the spatial coordinates of the non-structural element unit, the units with different resistivity intervals are rendered on the two-dimensional section plane to form a resistivity profile expressed in color scale. In this way, the distribution effect of complex three-dimensional resistivity heterogeneity on the two-dimensional section can be intuitively expressed. Different colors correspond to different resistivity values and ranges, showing the morphology of the resistivity body under unstructured geological conditions. This imaging result is more intuitive and convenient for subsequent geological analysis.

在上述实施例的基础上,生成非结构化地质环境的电阻率层图像的具体步骤还包括S41至S42:Based on the above embodiment, the specific steps of generating the resistivity layer image of the unstructured geological environment also include S41 to S42:

S41:获取各采集层的地质信息;根据各第二非结构元的第二电阻率和采集深度,生成各采集层的电阻率分布子图像。S41: Acquire geological information of each acquisition layer; generate a resistivity distribution sub-image of each acquisition layer according to the second resistivity and acquisition depth of each second non-structural element.

其中,地质信息是指对某个地理区域或地质体的构造、性质、成分等方面的信息,主要包括以下几类:地层信息,如地层的名称、序号,岩性,地层关系等。构造信息,如褶曲、断裂、晶间裂隙等地质构造等信息。Among them, geological information refers to information on the structure, properties, composition, etc. of a certain geographical area or geological body, mainly including the following categories: Stratigraphic information, such as the name, sequence number, lithology, stratigraphic relationship, etc. Structural information, such as folds, faults, intercrystalline cracks and other geological structures.

示例性地,从钻孔剖面等获取研究区域各采集层的岩性、物性等地质信息。计算机将第二非结构元中的第三电阻率数据根据采集深度进行分层,并结合各层岩性特征,选择与之对应电阻率值的颜色标识。生成电阻率针对各个采集层的分布子图像。这样在电阻率层析图像上就可以清楚看到不同采集层的电阻率分布状况,各层图像背后的地质属性也可以对应上,大大提高了图像的可解析性。例如哪些颜色表示煤层,哪些为砂岩等。通过这种结合已知地质信息与电阻率的表达手段,可以使电阻率层析成像对复杂环境的解析更具针对性,为后续地质资源的开发利用提供更实用的信息。Exemplarily, geological information such as lithology and physical properties of each acquisition layer in the study area is obtained from the borehole profile. The computer stratifies the third resistivity data in the second non-structural element according to the acquisition depth, and selects the color identification corresponding to the resistivity value in combination with the lithological characteristics of each layer. Generate a resistivity distribution sub-image for each acquisition layer. In this way, the resistivity distribution of different acquisition layers can be clearly seen on the resistivity tomography image, and the geological attributes behind each layer of the image can also be matched, which greatly improves the image's parseability. For example, which colors represent coal seams and which are sandstones. Through this expression method that combines known geological information with resistivity, the resistivity tomography can be made more targeted in analyzing complex environments, providing more practical information for the subsequent development and utilization of geological resources.

S42:根据各电阻率分布子图像和各地质信息,生成复杂结构地质条件的视电阻率层图像。S42: Generate an apparent resistivity layer image of complex structural geological conditions based on each resistivity distribution sub-image and each geological information.

示例性地,在各个采集层的电阻率分布子图像中可以标注出对应地层的岩性或物性信息。然后,计算机集成这些子图像到统一的三维坐标下,并根据它们的空间位置关系重新构建一个综合的电阻率剖面图。这样在一个电阻率层析图像上就同时包含了定量电阻率结果和定性地质信息。这种结合使图像既具有电阻率反演的高分辨率精细特征,也具有地质信息的直接可读性。直观地展示了不同地层与电阻率体的对应关系,大大提高了成像的准确性。For example, the lithology or physical property information of the corresponding strata can be marked in the resistivity distribution sub-images of each acquisition layer. Then, the computer integrates these sub-images into a unified three-dimensional coordinate system and reconstructs a comprehensive resistivity profile according to their spatial position relationship. In this way, both quantitative resistivity results and qualitative geological information are included in a resistivity tomography image. This combination enables the image to have both the high-resolution and fine features of resistivity inversion and the direct readability of geological information. The corresponding relationship between different strata and resistivity bodies is intuitively displayed, greatly improving the accuracy of imaging.

参照图2,为本申请实施例提供的一种空间自适应非结构元电阻率层析成像系统,该系统包括:数据获取模块、第一反演模块、第二反演模块,图像生成模块,其中:2 , a spatially adaptive non-structured element resistivity tomography system is provided in an embodiment of the present application. The system includes: a data acquisition module, a first inversion module, a second inversion module, and an image generation module, wherein:

数据获取模块,用于获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各电阻率和各采集深度,生成各采集层对应的第一非结构元;A data acquisition module, used to acquire the resistivity and acquisition depth of multiple acquisition layers under complex structural geological conditions, and generate a first non-structural element corresponding to each acquisition layer according to each resistivity and each acquisition depth;

第一反演模块,用于对各第一非结构元进行反演操作,得到各第一非结构元的第一电阻率,并根据第一电阻率调整第一非结构元,得到第二非结构元;A first inversion module is used to perform an inversion operation on each first non-structure element to obtain a first resistivity of each first non-structure element, and adjust the first non-structure element according to the first resistivity to obtain a second non-structure element;

第二反演模块,用于对各第二非结构元进行反演操作,得到各第二非结构元的第二电阻率;A second inversion module is used to perform an inversion operation on each second non-structure element to obtain a second resistivity of each second non-structure element;

图像生成模块,用于重复反演操作,直到反演操作达到预设迭代终止条件,根据各第二非结构元的第二电阻率和采集深度,生成复杂结构地质条件的视电阻率层图像。The image generation module is used to repeat the inversion operation until the inversion operation reaches a preset iteration termination condition, and generate an apparent resistivity layer image of complex structural geological conditions according to the second resistivity and acquisition depth of each second non-structural element.

在上述实施例的基础上,数据获取模块还用于根据各第一电阻率和各采集深度,匹配各采集层的深度分辨率和水平分辨率;根据各第一电阻率在各采集层中的分布位置,确定网格形状;根据深度分辨率、水平分辨率和网格形状,生成各采集层对应的第一非结构元。On the basis of the above embodiment, the data acquisition module is also used to match the depth resolution and horizontal resolution of each acquisition layer according to each first resistivity and each acquisition depth; determine the grid shape according to the distribution position of each first resistivity in each acquisition layer; generate the first non-structural element corresponding to each acquisition layer according to the depth resolution, horizontal resolution and grid shape.

在上述实施例的基础上,数据获取模块还包括获取各采集层中水平方向的初始采集点位置和初始深度位置;根据各初始采集点位置和各分布位置中电阻率的水平位置,确定各电阻率的第一网格坐标;根据各初始深度位置和各分布位置中电阻率的深度位置,确定各电阻率的第二网格坐标;根据各第一网格坐标和各第二网格坐标,将各电阻率匹配至对应的第一非结构元。On the basis of the above embodiment, the data acquisition module also includes acquiring the initial acquisition point position and the initial depth position in the horizontal direction in each acquisition layer; determining the first grid coordinates of each resistivity according to each initial acquisition point position and the horizontal position of the resistivity in each distribution position; determining the second grid coordinates of each resistivity according to each initial depth position and the depth position of the resistivity in each distribution position; and matching each resistivity to the corresponding first non-structural element according to each first grid coordinate and each second grid coordinate.

在上述实施例的基础上,第一反演模块还用于根据各第一非结构元中的各电阻率和第一非结构元的空间范围,确定各第一非结构元的目标电阻率;根据目标电阻率和采集深度,确定地表电位测量值;若地表电位测量值未超出标准电位值范围,则将目标电阻率作为第一电阻率。Based on the above embodiments, the first inversion module is also used to determine the target resistivity of each first non-structure element according to the resistivity in each first non-structure element and the spatial range of the first non-structure element; determine the surface potential measurement value according to the target resistivity and the acquisition depth; if the surface potential measurement value does not exceed the standard potential value range, the target resistivity is used as the first resistivity.

在上述实施例的基础上,第一反演模块还包括根据第一电阻率和电阻率,确定目标差值;若目标差值大于预设差值,则将第一非结构元作为待调整非结构元;根据目标差值,确定修正系数,并根据修正系数和待调整非结构元,确定第二非结构元。Based on the above embodiment, the first inversion module also includes determining a target difference based on the first resistivity and the resistivity; if the target difference is greater than the preset difference, the first non-structural element is used as the non-structural element to be adjusted; according to the target difference, a correction coefficient is determined, and according to the correction coefficient and the non-structural element to be adjusted, a second non-structural element is determined.

在上述实施例的基础上,第一反演模块还包括根据目标差值,确定空间范围修正值;根据空间范围修正值和第一非结构元,生成第二非结构元。Based on the above embodiment, the first inversion module further includes determining a spatial range correction value according to the target difference; and generating a second non-structure element according to the spatial range correction value and the first non-structure element.

在上述实施例的基础上,图像生成模块还用于获取各采集层的地质信息;根据各第二非结构元的第二电阻率和采集深度,生成各采集层的电阻率分布子图像;根据各电阻率分布子图像和各地质信息,生成复杂结构地质条件的视电阻率层图像。Based on the above embodiments, the image generation module is also used to obtain geological information of each acquisition layer; generate a resistivity distribution sub-image of each acquisition layer according to the second resistivity and acquisition depth of each second non-structural element; and generate an apparent resistivity layer image of complex structural geological conditions according to each resistivity distribution sub-image and various geological information.

需要说明的是:上述实施例提供的装置在实现其功能时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的装置和方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。It should be noted that: when the device provided in the above embodiment realizes its function, only the division of the above functional modules is used as an example. In actual application, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the device and method embodiments provided in the above embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

在本申请所提供的几种实施方式中,应该理解到,所披露的装置,可通过其他的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些服务接口,装置或单元的间接耦合或通信连接,可以是电性或其他的形式。In the several implementation modes provided in this application, it should be understood that the disclosed devices can be implemented in other ways. For example, the device embodiments described above are only schematic, such as the division of units, which is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some service interfaces, and the indirect coupling or communication connection of devices or units can be electrical or other forms.

作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储器包括:U盘、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a memory and includes several instructions for a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the various embodiments of the present application. The aforementioned memory includes: various media that can store program codes, such as USB flash drives, mobile hard drives, magnetic disks or optical disks.

以上者,仅为本公开的示例性实施例,不能以此限定本公开的范围。即但凡依本公开教导所作的等效变化与修饰,皆仍属本公开涵盖的范围内。本领域技术人员在考虑说明书及实践真理的公开后,将容易想到本公开的其他实施方案。The above are only exemplary embodiments of the present disclosure and cannot be used to limit the scope of the present disclosure. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure. After considering the disclosure of the specification and the truth of practice, those skilled in the art will easily think of other embodiments of the present disclosure.

本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未记载的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的范围和精神由权利要求限定。This application is intended to cover any variation, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the art not described in the present disclosure. The description and examples are to be regarded as exemplary only, and the scope and spirit of the present disclosure are defined by the claims.

Claims (9)

1.一种空间自适应非结构元电阻率层析成像方法,其特征在于,包括:1. A spatially adaptive non-structured element resistivity tomography method, comprising: 获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各所述电阻率和各所述采集深度,生成各所述采集层对应的第一非结构元;Obtaining resistivity and acquisition depths of multiple acquisition layers under complex structural geological conditions, and generating first non-structural elements corresponding to each acquisition layer according to each resistivity and each acquisition depth; 对各所述第一非结构元进行反演操作,得到各所述第一非结构元的第一电阻率,并根据所述第一电阻率调整所述第一非结构元,得到第二非结构元;Performing an inversion operation on each of the first non-structure elements to obtain a first resistivity of each of the first non-structure elements, and adjusting the first non-structure element according to the first resistivity to obtain a second non-structure element; 对各所述第二非结构元进行反演操作,得到各所述第二非结构元的第二电阻率;Performing an inversion operation on each of the second non-structural elements to obtain a second resistivity of each of the second non-structural elements; 重复所述反演操作,直到所述反演操作达到预设迭代终止条件,根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成所述复杂结构地质条件的视电阻率层图像。The inversion operation is repeated until the inversion operation reaches a preset iteration termination condition, and an apparent resistivity layer image of the complex structural geological condition is generated according to the second resistivity of each of the second non-structural elements and the acquisition depth. 2.根据权利要求1所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述根据各所述电阻率和各所述采集深度,生成各所述采集层对应的第一非结构元,包括:2. The spatially adaptive non-structured element resistivity tomography method according to claim 1, characterized in that the step of generating the first non-structured element corresponding to each acquisition layer according to each resistivity and each acquisition depth comprises: 根据各所述第一电阻率和各所述采集深度,匹配各所述采集层的深度分辨率和水平分辨率;According to each of the first resistivities and each of the acquisition depths, matching the depth resolution and the horizontal resolution of each of the acquisition layers; 根据各所述第一电阻率在各所述采集层中的分布位置,确定网格形状;Determining a grid shape according to the distribution position of each of the first resistivities in each of the acquisition layers; 根据所述深度分辨率、所述水平分辨率和所述网格形状,生成各所述采集层对应的所述第一非结构元。The first non-structural elements corresponding to each of the acquisition layers are generated according to the depth resolution, the horizontal resolution and the grid shape. 3.根据权利要求2所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述根据所述深度分辨率、所述水平分辨率和所述网格形状,生成各所述采集层对应的所述第一非结构元之后,还包括:3. The spatially adaptive non-structured element resistivity tomography method according to claim 2, characterized in that after generating the first non-structured element corresponding to each acquisition layer according to the depth resolution, the horizontal resolution and the grid shape, the method further comprises: 获取各所述采集层中水平方向的初始采集点位置和初始深度位置;Acquire the initial acquisition point position and initial depth position in the horizontal direction of each acquisition layer; 根据各所述初始采集点位置和各所述分布位置中所述电阻率的水平位置,确定各所述电阻率的第一网格坐标;Determining first grid coordinates of each of the resistivities according to the positions of each of the initial acquisition points and the horizontal positions of the resistivities in each of the distribution positions; 根据各所述初始深度位置和各所述分布位置中所述电阻率的深度位置,确定各所述电阻率的第二网格坐标;Determining second grid coordinates of each of the resistivities according to each of the initial depth positions and the depth position of the resistivity in each of the distribution positions; 根据各所述第一网格坐标和各所述第二网格坐标,将各所述电阻率匹配至对应的所述第一非结构元。According to each of the first grid coordinates and each of the second grid coordinates, each of the resistivities is matched to the corresponding first non-structural element. 4.根据权利要求1所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述对各所述第一非结构元进行反演操作,得到各所述第一非结构元的第一电阻率,包括:4. The spatially adaptive non-structured element resistivity tomography method according to claim 1, wherein the performing an inversion operation on each of the first non-structured elements to obtain the first resistivity of each of the first non-structured elements comprises: 根据各所述第一非结构元中的各所述电阻率和所述第一非结构元的空间范围,确定各所述第一非结构元的目标电阻率;determining a target resistivity of each of the first non-structure elements according to each of the resistivities in each of the first non-structure elements and a spatial range of the first non-structure elements; 根据所述目标电阻率和所述采集深度,确定地表电位测量值;Determining a surface potential measurement value according to the target resistivity and the acquisition depth; 若所述地表电位测量值未超出标准电位值范围,则将所述目标电阻率作为所述第一电阻率。If the measured value of the ground surface potential does not exceed the standard potential value range, the target resistivity is used as the first resistivity. 5.根据权利要求1所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述根据所述第一电阻率调整所述第一非结构元,得到第二非结构元,包括:5. The spatially adaptive non-structured element resistivity tomography method according to claim 1, wherein the adjusting the first non-structured element according to the first resistivity to obtain the second non-structured element comprises: 根据所述第一电阻率和所述电阻率,确定目标差值;determining a target difference according to the first resistivity and the resistivity; 若所述目标差值大于预设差值,则将所述第一非结构元作为待调整非结构元;If the target difference is greater than the preset difference, the first non-structural element is used as the non-structural element to be adjusted; 根据所述目标差值,确定修正系数,并根据所述修正系数和所述待调整非结构元,确定所述第二非结构元。A correction coefficient is determined according to the target difference, and the second non-structural element is determined according to the correction coefficient and the non-structural element to be adjusted. 6.根据权利要求5所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述根据所述目标差值,确定修正系数,并根据所述修正系数和所述待调整非结构元,确定所述第二非结构元,包括:6. The spatially adaptive non-structure element resistivity tomography method according to claim 5, characterized in that the step of determining a correction coefficient according to the target difference, and determining the second non-structure element according to the correction coefficient and the non-structure element to be adjusted comprises: 根据所述目标差值,确定空间范围修正值;Determining a spatial range correction value according to the target difference; 根据所述空间范围修正值和所述第一非结构元,生成所述第二非结构元。The second non-structure element is generated according to the spatial range correction value and the first non-structure element. 7.根据权利要求1所述的空间自适应非结构元电阻率层析成像方法,其特征在于,所述根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成所述复杂结构地质条件的视电阻率层图像,包括:7. The spatially adaptive non-structured element resistivity tomography method according to claim 1, characterized in that the step of generating the apparent resistivity layer image of the complex structural geological condition according to the second resistivity of each second non-structured element and the acquisition depth comprises: 获取各所述采集层的地质信息;Acquiring geological information of each of the acquisition layers; 根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成各所述采集层的电阻率分布子图像;generating a resistivity distribution sub-image of each acquisition layer according to the second resistivity of each second non-structural element and the acquisition depth; 根据各所述电阻率分布子图像和各所述地质信息,生成所述复杂结构地质条件的视电阻率层图像。An apparent resistivity layer image of the complex structural geological condition is generated according to each of the resistivity distribution sub-images and each of the geological information. 8.一种空间自适应非结构元电阻率层析成像系统,其特征在于,所述系统包括:8. A spatially adaptive non-structured element resistivity tomography system, characterized in that the system comprises: 数据获取模块,用于获取复杂结构地质条件的在多个采集层的电阻率和采集深度,根据各所述电阻率和各所述采集深度,生成各所述采集层对应的第一非结构元;A data acquisition module, used to acquire the resistivity and acquisition depth of multiple acquisition layers under complex structural geological conditions, and generate a first non-structural element corresponding to each acquisition layer according to each resistivity and each acquisition depth; 第一反演模块,用于对各所述第一非结构元进行反演操作,得到各所述第一非结构元的第一电阻率,并根据所述第一电阻率调整所述第一非结构元,得到第二非结构元;A first inversion module is used to perform an inversion operation on each of the first non-structure elements to obtain a first resistivity of each of the first non-structure elements, and adjust the first non-structure element according to the first resistivity to obtain a second non-structure element; 第二反演模块,用于对各所述第二非结构元进行反演操作,得到各所述第二非结构元的第二电阻率;A second inversion module is used to perform an inversion operation on each of the second non-structure elements to obtain a second resistivity of each of the second non-structure elements; 图像生成模块,用于重复所述反演操作,直到所述反演操作达到预设迭代终止条件,根据各所述第二非结构元的所述第二电阻率和所述采集深度,生成所述复杂结构地质条件的视电阻率层图像。An image generation module is used to repeat the inversion operation until the inversion operation reaches a preset iteration termination condition, and generate an apparent resistivity layer image of the complex structural geological condition according to the second resistivity of each second non-structural element and the acquisition depth. 9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有指令,当所述指令被执行时,执行如权利要求1-7任意一项所述的空间自适应非结构元电阻率层析成像方法步骤。9. A computer-readable storage medium, characterized in that the computer-readable storage medium stores instructions, and when the instructions are executed, the steps of the spatially adaptive non-structured element resistivity tomography method according to any one of claims 1 to 7 are executed.
CN202410126376.0A 2024-01-30 2024-01-30 Space self-adaptive unstructured element resistivity tomography method and system Active CN117991348B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410126376.0A CN117991348B (en) 2024-01-30 2024-01-30 Space self-adaptive unstructured element resistivity tomography method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410126376.0A CN117991348B (en) 2024-01-30 2024-01-30 Space self-adaptive unstructured element resistivity tomography method and system

Publications (2)

Publication Number Publication Date
CN117991348A true CN117991348A (en) 2024-05-07
CN117991348B CN117991348B (en) 2025-02-21

Family

ID=90886692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410126376.0A Active CN117991348B (en) 2024-01-30 2024-01-30 Space self-adaptive unstructured element resistivity tomography method and system

Country Status (1)

Country Link
CN (1) CN117991348B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102419455A (en) * 2011-08-23 2012-04-18 安徽理工大学 Interwell parallel resistivity CT (computed tomography) testing method
CN102495428A (en) * 2011-12-12 2012-06-13 山东大学 Resistivity real-time imaging monitoring method and system for water-bursting geological disaster in construction period of underground engineering
US20160033668A1 (en) * 2014-07-30 2016-02-04 Battelle Memorial Institute Method of imaging the electrical conductivity distribution of a subsurface
CN108169801A (en) * 2018-01-16 2018-06-15 陕西铁道工程勘察有限公司 High-resolution earth resistivity fast imaging method
CN109919844A (en) * 2019-02-28 2019-06-21 河南师范大学 A High-Resolution Method for Conductivity Distribution Reconstruction in Electrical Tomography
CN112666612A (en) * 2020-11-02 2021-04-16 中国铁路设计集团有限公司 Magnetotelluric two-dimensional inversion method based on tabu search
US20220308251A1 (en) * 2020-06-03 2022-09-29 Shandong University Three-dimensional electrical resistivity tomography method and system
CN115166840A (en) * 2022-06-10 2022-10-11 阿里云计算有限公司 Method and device for realizing resistivity inversion
CN115238549A (en) * 2022-07-25 2022-10-25 中南大学 A method of using ERT to monitor the safety factor under the condition of landslide rainfall
CN115238550A (en) * 2022-07-25 2022-10-25 中南大学 Self-adaptive unstructured grid landslide rainfall geoelectric field numerical simulation calculation method
CN115689802A (en) * 2022-10-26 2023-02-03 重庆大学 A Multi-resolution HVDC Ground Current 3D Model Calculation Method
US20230119421A1 (en) * 2021-10-14 2023-04-20 Earthsystems Technologies, Inc. System and method for resistivity interpretation
CN116704147A (en) * 2023-07-03 2023-09-05 广西龙马高速公路有限公司 Electromagnetic exploration three-dimensional modeling method for complex geological structure

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102419455A (en) * 2011-08-23 2012-04-18 安徽理工大学 Interwell parallel resistivity CT (computed tomography) testing method
CN102495428A (en) * 2011-12-12 2012-06-13 山东大学 Resistivity real-time imaging monitoring method and system for water-bursting geological disaster in construction period of underground engineering
US20160033668A1 (en) * 2014-07-30 2016-02-04 Battelle Memorial Institute Method of imaging the electrical conductivity distribution of a subsurface
CN108169801A (en) * 2018-01-16 2018-06-15 陕西铁道工程勘察有限公司 High-resolution earth resistivity fast imaging method
CN109919844A (en) * 2019-02-28 2019-06-21 河南师范大学 A High-Resolution Method for Conductivity Distribution Reconstruction in Electrical Tomography
US20220308251A1 (en) * 2020-06-03 2022-09-29 Shandong University Three-dimensional electrical resistivity tomography method and system
CN112666612A (en) * 2020-11-02 2021-04-16 中国铁路设计集团有限公司 Magnetotelluric two-dimensional inversion method based on tabu search
US20230119421A1 (en) * 2021-10-14 2023-04-20 Earthsystems Technologies, Inc. System and method for resistivity interpretation
CN115166840A (en) * 2022-06-10 2022-10-11 阿里云计算有限公司 Method and device for realizing resistivity inversion
CN115238549A (en) * 2022-07-25 2022-10-25 中南大学 A method of using ERT to monitor the safety factor under the condition of landslide rainfall
CN115238550A (en) * 2022-07-25 2022-10-25 中南大学 Self-adaptive unstructured grid landslide rainfall geoelectric field numerical simulation calculation method
CN115689802A (en) * 2022-10-26 2023-02-03 重庆大学 A Multi-resolution HVDC Ground Current 3D Model Calculation Method
CN116704147A (en) * 2023-07-03 2023-09-05 广西龙马高速公路有限公司 Electromagnetic exploration three-dimensional modeling method for complex geological structure

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CAI YUMEI: "Forward Modeling Analysis of Cross-hole Resistivity Tomography (ERT) Using Unstructured Grids", JOURNAL OF PHYSICS: CONFERENCE SERIES, 25 January 2024 (2024-01-25), pages 012056 *
冯锐, 李智明, 李志武, 周燕云: "电阻率层析成像技术", 中国地震, no. 01, 25 March 2004 (2004-03-25) *
尚耀军: "高密度电法在贵州金沙某滑坡勘察中的应用", 工程地球物理学报, vol. 10, no. 6, 30 November 2013 (2013-11-30), pages 771 - 776 *
马鑫: "井间地震波层析成像与直流电阻率法联合反演", 西部探矿工程, vol. 31, no. 06, 14 June 2019 (2019-06-14), pages 107 - 112 *
高级;张海江;秦臻;: "不同观测方式下基于伴随矩阵的2.5D全通道电阻率反演研究", 地球物理学进展, no. 06, 15 December 2016 (2016-12-15) *

Also Published As

Publication number Publication date
CN117991348B (en) 2025-02-21

Similar Documents

Publication Publication Date Title
US11042676B2 (en) Representing structural uncertainty in a mesh representing a geological environment
RU2337404C1 (en) Computer method for modelling during drilling and stratified underground fraction visualisation
CA2920499C (en) Stratigraphic function
CN111612903B (en) Geological data visualization method based on mixed data model
EP2869096B1 (en) Systems and methods of multi-scale meshing for geologic time modeling
EP2810101B1 (en) Improving efficiency of pixel-based inversion algorithms
CN110823962B (en) Three-dimensional imaging method and system for landslide mass
US12282128B2 (en) System and method for quantitative seismic integration modeling workflow
CN102749647B (en) Quantitative method for space of holes in reservoir
US20160202390A1 (en) In-Situ Wellbore, Core and Cuttings Information System
NO20111053A1 (en) Simultaneous removal of structural inclination angle data
WO2017127058A1 (en) Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations
US20200379135A1 (en) Analytics and machine learning method for estimating petrophysical property values
CN119126212B (en) High-precision detection method, system, equipment and storage medium for urban underground space
CN115100373B (en) Three-dimensional integrated characterization method, device, medium and equipment for surface and underground environment
WO2014099201A1 (en) Geophysical modeling of subsurface volumes based on horizon extraction
WO2017131608A1 (en) Permeability anisotropy assessment in subsurface anisotropic formations
CN118444367B (en) A geological survey imaging method and medium based on cross-well acoustic CT
CN110333543A (en) Interpretation and Imaging Method and System for Low Resistivity Body Based on Reflection Coefficient Analysis
CN116071510A (en) Marine wind power engineering geological model based on geological survey data consistency and modeling method
CN107748393B (en) A kind of bearing calibration based on the stratigraphic dip of numerical simulation to resistivity effects
CN117991348A (en) Space self-adaptive unstructured element resistivity tomography method and system
CN111965707A (en) A seismic inversion reservoir prediction method for complex structures with overburden faults
CN116341232A (en) Method and system for automatically generating horizons in sea breeze geological modeling based on image recognition
US20140185413A1 (en) Dip seismic attributes

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