Automatic modeling method and system for irregular interface structure and water body nesting

Abstract

The present invention belongs to the field of geophysical exploration technology. To address the significant deviation between current geological models and actual working conditions, a method and system for automatically modeling irregular interface structures and water body nesting is provided. The method includes constructing spatial correlation rules for stratum position-amplitude fluctuation, spatial correlation rules for stratum position-water-bearing caves, and a table of correspondences between lithology-surrounding rock classification-physical property parameters. This method then generates three-dimensional irregular stratum models and irregular water-bearing structure models for different undulating interfaces. Under the same grid coordinates, the physical property parameters of the irregular water-bearing structure model are directly mapped to the corresponding positions in the three-dimensional irregular stratum model, forming a composite model of stratum interface fluctuation and water body nesting. The resulting geological model better reflects actual working conditions, providing effective data support for the prevention and control of geological disasters.

Description

Technical Field

[0001] The present invention belongs to the technical field of geophysical exploration, and in particular relates to an automatic modeling method and system for irregular interface structure and water body nesting. Background Art

[0002] The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.

[0003] As tunnel construction extends into geologically complex areas, the contradiction between complex geological conditions and the need for accurate forecasting becomes increasingly prominent, and the need to construct complex geological models becomes increasingly urgent. In actual geological environments, the interweaving of different geological interfaces, fault structures, and diverse water bodies creates an extremely complex geological landscape, which increases the difficulty of detecting hazard sources. Current geological preliminary survey data is limited. Developing geological models solely based on this limited data results in significant deviations from actual operating conditions, making it impossible to provide effective data support for the prevention and control of geological disasters. Summary of the Invention

[0004] In order to solve the technical problems existing in the above-mentioned background technology, the present invention provides a method and system for automatic modeling of irregular interface structure and water body nesting, which can realize automatic modeling of irregular interface structure and water body nesting, so that the generated geological model is more in line with the actual working conditions, and provides effective data support for the prevention and control of geological disasters.

[0005] In order to achieve the above object, the present invention adopts the following technical solutions:

[0006] A first aspect of the present invention provides a method for automatically modeling irregular interface structures and water body nesting.

[0007] An automatic modeling method for irregular interface structure and water body nesting, comprising:

[0008] Acquire multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This allows for the construction of spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters.

[0009] Based on the spatial correlation rules between stratum position and undulation amplitude, a three-dimensional irregular stratum model with different undulation interfaces was constructed. Based on the spatial correlation rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model was constructed by combining multiple ellipsoids.

[0010] Under the same grid coordinates, the physical parameters of the irregular water-bearing structure model are directly mapped to the corresponding positions of the three-dimensional irregular stratum model, forming a composite model of stratum interface fluctuations and water body nesting.

[0011] As an implementation method, the process of constructing the spatial association rules between stratum position and fluctuation amplitude is as follows:

[0012] Statistical analysis of stratum depth, interface fluctuation distribution and fluctuation amplitude from drill core data;

[0013] At the fault location, the inverse distance weighted method of spatial interpolation is used to generate a continuous undulation surface, and then the probability density distribution of the stratum interface undulation and the undulation wavelength range are calculated, and the spatial correlation rules between stratum position and undulation amplitude are obtained.

[0014] As an implementation method, the process of generating a continuous relief surface using the inverse distance weighted method of spatial interpolation is as follows:

[0015] According to each drilling The coordinates of , ) and stratum depth, calculate the average depth of the area, and record the depth of each borehole Fluctuation relative to the average depth ;

[0016] Taking the regional average depth as the baseline, calculate the points on the regional average depth line ( , ) to drill hole distance : ;

[0017] Calculate interpolation ; ; ;in, is the total weight; is the sub-weight;

[0018] Based on drilling Belonging coordinate system, based on drilling Fluctuation relative to the average depth For values, a continuous undulating surface is generated.

[0019] As an embodiment, the water-bearing structural features include attribute information and location information of water-bearing caves. The attribute information of water-bearing caves is stored in the form of an attribute table and associated with the spatial location information to obtain the spatial association rules between the formation location and the water-bearing caves.

[0020] As an implementation method, the process of calibrating the physical parameter threshold is as follows:

[0021] Classify lithology and surrounding rock according to physical parameters;

[0022] Calculate the mean and standard deviation of physical property parameters of similar lithologies, and then determine the normal distribution range;

[0023] The physical property boundaries of different lithologies are defined by the normal distribution range, and the physical property parameter thresholds are calibrated by combining the measured data with engineering experience data.

[0024] As an implementation method, the process of constructing a three-dimensional irregular stratum model with different undulating interfaces is as follows:

[0025] According to the detection resolution requirements, the grid size is preset and a two-dimensional uniform initial model is formed accordingly;

[0026] According to the preset number of strata and the height of each stratum, the two-dimensional uniform initial model is divided into a number of rectangular blocks in the vertical direction; wherein the number of rectangular blocks is the same as the number of strata;

[0027] Configure the corresponding physical parameters for each rectangular block, copy the two-dimensional uniform initial model along the set direction, and obtain a corresponding three-dimensional array;

[0028] For each replicated two-dimensional uniform initial model, the maximum disturbance value and the corresponding random distribution function of each interface are set according to the spatial correlation rule between the formation position and the fluctuation amplitude;

[0029] A random disturbance value is added to the upper and lower interfaces of each rectangular block in the vertical direction. According to the position of the interface after disturbance, the attribute value of the corresponding position in the three-dimensional array is updated to form a three-dimensional irregular stratum model with different undulating interfaces.

[0030] As an implementation method, the process of constructing an irregular water-bearing structure model by combining multiple ellipsoids is as follows:

[0031] Based on the spatial distribution density, morphological parameters, and connectivity characteristics of water-bearing structures, and in accordance with the spatial correlation rules between stratum position and water-bearing caves, the center coordinate positions, rotation angles, and semi-axis length ratios of multiple ellipsoids were randomly generated within a preset range on both sides of the fault in the area ahead of the tunnel face to simulate inclined water channels and beaded caves.

[0032] According to the central coordinate positions of multiple ellipsoids and the corresponding rotation angles and semi-axis length ratios, the corresponding ellipsoid models are constructed, and these ellipsoid models are combined and remapped to the original coordinate system to form an irregular water-bearing structure model.

[0033] A second aspect of the present invention provides an automatic modeling system for irregular interface structure and water body nesting.

[0034] An irregular interface structure and water body nesting automatic modeling system, comprising:

[0035] The data processing module is used to obtain multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This module then constructs spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters.

[0036] The model construction module is used to construct a three-dimensional irregular stratum model with different undulating interfaces based on the spatial association rules between stratum position and undulation amplitude; based on the spatial association rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model is constructed by combining multiple ellipsoids;

[0037] The model fusion module is used to directly map the physical parameters of the irregular water-bearing structure model to the corresponding positions of the three-dimensional irregular stratum model under the same grid coordinates, forming a composite model with stratum interface fluctuations and water body nesting.

[0038] A third aspect of the present invention provides a computer-readable storage medium.

[0039] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps in the above-mentioned method for automatic modeling of irregular interface construction and water body nesting.

[0040] A fourth aspect of the present invention provides a computer device.

[0041] A computer device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the above-mentioned method for automatic modeling of irregular interface construction and water body nesting are implemented.

[0042] Compared with the prior art, the present invention has the following beneficial effects:

[0043] The present invention uses multi-source data from tunnels to construct spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical property parameters. Based on this, a three-dimensional irregular stratum model with different fluctuation interfaces is constructed, and an irregular water-bearing structure model is constructed by combining multiple ellipsoids. Finally, the physical property parameters of the irregular water-bearing structure model are directly mapped to the corresponding positions of the three-dimensional irregular stratum model, forming a composite model with stratum interface fluctuations and water bodies nested. This model can more accurately and comprehensively simulate geological models and comprehensively and realistically reflect actual complex geological conditions, providing accurate model support for geological research and disaster prevention.

[0044] Advantages of additional aspects of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] The accompanying drawings, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0046] Figure 1 This is a general flow chart of the automatic modeling of irregular interface structure and water body nesting in an embodiment of the present invention;

[0047] Figure 2 This is a process for constructing a three-dimensional irregular stratum model with different undulating interfaces in an embodiment of the present invention;

[0048] Figure 3 The water-bearing structure modeling process in the embodiment of the present invention;

[0049] Figure 4 Schematic diagram of the three-dimensional irregular stratum model result in an embodiment of the present invention;

[0050] Figure 5 This is a schematic diagram of the randomly generated complex model results in an embodiment of the present invention. DETAILED DESCRIPTION

[0051] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0052] It should be noted that the following detailed descriptions are illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which the present invention belongs.

[0053] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.

[0054] Example 1

[0055] Combine Figure 1 The embodiment of the present invention provides an automatic modeling method for irregular interface structure and water body nesting, which includes:

[0056] S101: Acquire multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This allows for the construction of spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters.

[0057] In the specific implementation process of step S101, the process of constructing the spatial association rule between the stratum position and the fluctuation amplitude is as follows:

[0058] Statistical analysis of stratum depth, interface fluctuation distribution and fluctuation amplitude from drill core data;

[0059] At the fault location, the inverse distance weighted method of spatial interpolation is used to generate a continuous undulation surface, and then the probability density distribution of the stratum interface undulation and the undulation wavelength range are calculated, and the spatial correlation rules between stratum position and undulation amplitude are obtained.

[0060] This embodiment uses the inverse distance weighted method of spatial interpolation to generate a continuous relief surface, which can solve the problem of data discontinuity at the fault location and achieve the purpose of improving the accuracy of subsequent model construction based on limited data.

[0061] The process of generating a continuous undulating surface using the inverse distance weighted method of spatial interpolation is as follows:

[0062] According to each drilling The coordinates of , ) and stratum depth, calculate the average depth of the area, and record the depth of each borehole Fluctuation relative to the average depth ;

[0063] Taking the regional average depth as the baseline, calculate the points on the regional average depth line ( , ) to drill hole distance : ;

[0064] Calculate interpolation ; ; ;in, is the total weight; is the sub-weight;

[0065] Based on drilling Belonging coordinate system, based on drilling Fluctuation relative to the average depth For values, a continuous undulating surface is generated.

[0066] Among them, all interpolation points are grouped and the probability density is obtained through the histogram. The distance between adjacent wave peaks can be measured on the undulation surface map to estimate the wavelength.

[0067] It should be noted that, in other embodiments, in addition to the spatial interpolation inverse distance weighted method to generate the continuous undulating surface, other existing interpolation methods can also be used to achieve the generation of the continuous undulating surface, which will not be described in detail here.

[0068] During the implementation process, geological radar (GPR) emits high-frequency electromagnetic waves into the ground. Based on the reflection characteristics of these waves in different media, the location, size, and shape of underground water-bearing caves are detected. Once areas potentially containing water-bearing caves are initially detected using GPR and other methods, drilling verification is performed. Drilling allows for direct observation of the cave's filling (whether it contains water or mud, etc.), measurement of parameters such as its height and diameter, and acquisition of information on the lithology and thickness of the stratum in which the cave resides, enabling differentiation between water-filled caves, water-conducting faults, and fractured water zones. The GPR detection data and drilling verification data are integrated to create a spatial database. This database records the precise location, size, filling, stratum name, and depth of each water-bearing cave.

[0069] Among them, the water-bearing structural characteristics include the attribute information and location information of water-bearing caves. The attribute information of water-bearing caves is stored in the form of an attribute table and associated with the spatial location information, thereby obtaining the spatial association rules between the stratigraphic location and the water-bearing caves, providing constraints for multi-interface nested modeling.

[0070] In the specific implementation process, the process of calibrating the threshold value of physical parameters is as follows:

[0071] Classify lithology and surrounding rock according to physical parameters;

[0072] Calculate the mean and standard deviation of physical property parameters of similar lithologies, and then determine the normal distribution range;

[0073] The physical property boundaries of different lithologies are defined by the normal distribution range, and the physical property parameter thresholds are calibrated by combining the measured data with engineering experience data.

[0074] By calibrating the threshold values ​​of physical property parameters, the present invention can accurately construct a correspondence table of lithology, surrounding rock classification and physical property parameters, thereby improving the accurate data basis for the construction of irregular water-bearing structure models.

[0075] S102: Based on the spatial association rules between stratum position and undulation amplitude, a three-dimensional irregular stratum model with different undulation interfaces is constructed; based on the spatial association rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model is constructed by combining multiple ellipsoids.

[0076] In the specific implementation process of step S102, combined with Figure 2 and Figure 4 , the process of constructing a three-dimensional irregular stratum model with different undulating interfaces is:

[0077] According to the detection resolution requirements, the grid size is preset and a two-dimensional uniform initial model is formed accordingly; for example, a rectangular area with a length of L and a width of R. Assume that the initial bottom layer is uniform in this two-dimensional area and has an initial property (such as resistivity, wave velocity, dielectric constant, etc.), which is recorded as P0. This two-dimensional area can be represented as a two-dimensional array , where each element The initial value is P0, i represents the row index, j represents the column index, .

[0078] According to the preset number of strata and the height of each stratum, the two-dimensional uniform initial model is divided into a number of rectangular blocks in the vertical direction; wherein the number of rectangular blocks is the same as the number of strata; for example, based on the preliminary geological data, such as electromagnetic, drilling, etc., the preliminary number of strata n and the height of each stratum h are obtained. i The height of each rectangular block is h i , ; Where H is the total height of the two-dimensional area in the vertical direction.

[0079] Configure the corresponding physical parameters for each rectangular block, copy the two-dimensional uniform initial model along the set direction (for example, the x direction) to obtain a corresponding three-dimensional array; for example, the number of copies is D, and a three-dimensional array is obtained. ;

[0080] For each replicated two-dimensional uniform initial model, the maximum disturbance value and the corresponding random distribution function of each interface are set according to the spatial correlation rule between the formation position and the fluctuation amplitude;

[0081] A random disturbance value is added to the upper and lower interfaces of each rectangular block in the vertical direction. According to the position of the interface after disturbance, the attribute value of the corresponding position in the three-dimensional array is updated to form a three-dimensional irregular stratum model with different undulating interfaces.

[0082] For example, set the maximum disturbance value of each interface and the corresponding random distribution function , add random perturbation values ​​in the vertical direction to the upper and lower interfaces of each rectangular block (except the upper interface of the top layer and the lower interface of the bottom layer). Let the upper interface of the k-th layer rectangular block be , the lower interface For each point on the interface , random perturbation values ​​are generated by a random function, such as The upper interface position after disturbance is , the lower interface position is . Update the two-dimensional array according to the interface position after disturbance If the interface position changes, the attribute values ​​of the elements above and below the new interface position are updated through linear interpolation according to the distance from the interface, thus forming a three-dimensional irregular stratum model.

[0083] In this embodiment, if Figure 3 As shown in Figure 2, the process of constructing an irregular water-bearing structure model by combining multiple ellipsoids is as follows:

[0084] Based on the spatial distribution density, morphological parameters, and connectivity characteristics of water-bearing structures, and in accordance with the spatial correlation rules between stratum position and water-bearing caves, the center coordinate positions, rotation angles, and semi-axis length ratios of multiple ellipsoids were randomly generated within a preset range on both sides of the fault in the area ahead of the tunnel face to simulate inclined water channels and beaded caves.

[0085] According to the central coordinate positions of multiple ellipsoids and the corresponding rotation angles and semi-axis length ratios, the corresponding ellipsoid models are constructed, and these ellipsoid models are combined and remapped to the original coordinate system to form an irregular water-bearing structure model.

[0086] This embodiment uses an ellipsoid modeling method to achieve accurate simulation of water-containing conditions such as caves and water-conducting channels in tunnels.

[0087] For example, based on the grid division of the three-dimensional complex stratigraphic model, three-dimensional grid coordinates x, y, and z are created according to the given grid shape (such as the common area of ​​the joint inversion studied is 30×30×30), and the commonly used meshgrid function can be selected based on the Python language to generate multidimensional grid coordinates.

[0088] Based on the spatial distribution density, morphological parameters, and connectivity characteristics of water-bearing structures obtained through initial data structured processing, and in accordance with the spatial association rule of "stratum position-water-bearing caves", the center coordinate positions of multiple ellipsoids were randomly generated within a certain range on both sides of the fault within a 30-meter area in front of the tunnel face, along with the corresponding rotation angles and semi-axis length ratios, to simulate inclined water-conducting channels and beaded caves.

[0089] The ellipsoid calculation formula is: ;

[0090] in,( x 0 ,y 0 ,z0) is the center coordinate of the ellipsoid, and a, b, and c are the lengths of the semi-axes along the x, y, and z axes, respectively. This formula describes the geometric conditions satisfied by points inside the ellipsoid, namely the relationship between the distance from the center of the ellipsoid and the length of the semi-axes.

[0091] In the process of creating the rotation matrix rotation ellipsoid, first, the single-axis rotation matrix is ​​set as:

[0092] ;

[0093] ;

[0094] ;

[0095] 、 and are the rotation angles along the x, y, and z axes respectively; the final matrix Obtained by multiplying the single-axis matrices in sequence, for example, rotating in the order Z→Y→X:

[0096] ;

[0097] The complete final matrix after expansion is:

[0098] ;

[0099] Set the parameter value of the internal area of ​​the ellipsoid to the specified physical parameter, that is, the left side of the ellipsoid formula is <1.

[0100] In this embodiment, the process of remapping the multi-ellipsoid combination to the original coordinate system is as follows:

[0101] First, the matrix coordinates are multiplied by the rotation matrix to obtain the rotated coordinates. After obtaining the rotated coordinates, the relative coordinates are obtained by dividing each coordinate point by the semi-axis length of the ellipsoid to normalize the ellipsoid, that is, each point is converted into a relative position relative to the ellipsoid.

[0102] Then calculate the sum of the squares of these relative coordinate points to the origin, and then take the square root to get the distance after rotation;

[0103] When performing distance calculations on a grid, the rotated coordinate points are resized to match the shape of the grid. This allows for a one-to-one correspondence between the rotated coordinate points and the grid points, so that the position of each grid point is consistent with the rotated coordinate point in space.

[0104] S103: Under the same grid coordinates, the physical property parameters of the irregular water-bearing structure model are directly mapped to corresponding positions of the three-dimensional irregular stratum model to form a composite model of stratum interface fluctuations and water body nesting.

[0105] Figure 5 For the same geological body, three corresponding complex model results of resistivity, wave velocity and dielectric constant are given, where the numerical range of resistivity is set to 0-1000Ω·m, the range of wave velocity is set to 0-3000m / s, and the range of dielectric constant is set to 4-10. Figure 5 It can be seen that such a model can reflect complex geological conditions more comprehensively and realistically, and can provide accurate model support for geological research and disaster prevention and control.

[0106] Due to the lack of an automated modeling framework, which has led to a data supply gap, deep learning model training requires a sufficient number of real-world datasets. However, the existing process has the following fundamental bottlenecks: manual modeling is inefficient and relies on engineers' experience to manually adjust parameters. Building a single complex model takes several hours, which makes it difficult to meet the deep learning requirements for tens of thousands of samples. In addition, the scene generalization ability is insufficient, and the model generation logic lacks initial data constraints, resulting in a large number of invalid samples of "simple bottom layer + regular water body" in the dataset, and the model has poor generalization ability in the measured data after training.

[0107] In order to solve the above problems, the present invention also uses the number of strata and the physical properties of the strata obtained during the tunnel construction process to adjust the relevant information in the preliminary survey data, compare it with the corresponding data in the preliminary survey data, calculate the deviation size, and correct the interface undulation parameters of the three-dimensional stratum model and the ellipsoid parameters of the water-bearing structure model in real time to make the composite model more in line with the actual address situation, dynamically adjust the composite model parameters to improve the modeling accuracy of hidden disaster sources, and construct a complete technical chain from data processing to model application, which significantly improves the modeling and accuracy of hidden disaster sources under complex geological conditions, and provides key technical support for intelligent disaster prevention in tunnel engineering.

[0108] Example 2

[0109] An embodiment of the present invention provides an automatic modeling system for irregular interface structure and water body nesting, which includes:

[0110] The data processing module is used to obtain multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This module then constructs spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters.

[0111] The model construction module is used to construct a three-dimensional irregular stratum model with different undulating interfaces based on the spatial association rules between stratum position and undulation amplitude; based on the spatial association rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model is constructed by combining multiple ellipsoids;

[0112] The model fusion module is used to directly map the physical parameters of the irregular water-bearing structure model to the corresponding positions of the three-dimensional irregular stratum model under the same grid coordinates, forming a composite model with stratum interface fluctuations and water body nesting.

[0113] It should be noted here that the various modules in the embodiment of the present invention correspond one-to-one to the various steps in the above embodiment, and their specific implementation processes are the same, which will not be described in detail here.

[0114] Example 3

[0115] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the steps in the above-mentioned method for automatic modeling of irregular interface construction and water body nesting are implemented.

[0116] Example 4

[0117] This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for automatic modeling of irregular interface construction and water body nesting as described above are implemented.

[0118] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Thus, the present invention may take the form of hardware embodiments, software embodiments, or embodiments combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage and optical storage) containing computer-usable program code.

[0119] The present invention is described with reference to flowcharts and / or block diagrams of computer program products of methods, devices (systems) according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1A device that provides the functions specified in a block or multiple blocks.

[0120] Those skilled in the art will appreciate that all or part of the processes in the above-described method embodiments can be implemented by instructing the relevant hardware through a computer program. The program can be stored in a computer-readable storage medium, and when executed, the program can include the processes in the above-described method embodiments. The storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM).

[0121] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection of the present invention.

Claims

1. A method for automatic modeling of irregular interface structure and water body nesting, characterized in that: include: Acquire multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This allows for the construction of spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters. Based on the spatial correlation rules between stratum position and undulation amplitude, a three-dimensional irregular stratum model with different undulation interfaces was constructed. Based on the spatial correlation rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model was constructed by combining multiple ellipsoids. Under the same grid coordinates, the physical parameters of the irregular water-bearing structure model are directly mapped to the corresponding positions of the 3D irregular stratum model, forming a composite model of stratum interface fluctuations and water body nesting. The process of constructing the spatial association rules between stratum position and fluctuation amplitude is as follows: Statistical analysis of stratum depth, interface fluctuation distribution and fluctuation amplitude from drill core data; At the fault location, the inverse distance weighted method of spatial interpolation is used to generate a continuous undulation surface, and then the probability density distribution of the stratum interface undulation and the undulation wavelength range are calculated, and the spatial correlation rules between stratum position and undulation amplitude are obtained.

2. The method for automatically modeling irregular interface structures and water body nesting according to claim 1, characterized in that: The process of generating a continuous undulating surface using the inverse distance weighted method of spatial interpolation is as follows: According to each drilling The coordinates of , ) and stratum depth, calculate the average depth of the area, and record the depth of each borehole Fluctuation relative to the average depth ; Taking the regional average depth as the baseline, calculate the points on the regional average depth line ( , ) to drill hole distance : ; Calculate interpolation ; ; ;in, is the total weight; is the sub-weight; Based on drilling Belonging coordinate system, based on drilling Fluctuation relative to the average depth For values, a continuous undulating surface is generated.

3. The method for automatically modeling irregular interface structures and water body nesting according to claim 1, characterized in that: The water-bearing structural features include attribute information and location information of water-bearing caves. The attribute information of water-bearing caves is stored in the form of an attribute table and associated with spatial location information to obtain spatial association rules between stratum location and water-bearing caves.

4. The method for automatically modeling irregular interface structures and water body nesting according to claim 1, wherein: The process of calibrating the threshold value of physical parameters is as follows: Classify lithology and surrounding rock according to physical parameters; Calculate the mean and standard deviation of physical property parameters of similar lithologies, and then determine the normal distribution range; The physical property boundaries of different lithologies are defined by the normal distribution range, and the physical property parameter thresholds are calibrated by combining the measured data with engineering experience data.

5. The method for automatically modeling irregular interface structure and water body nesting according to claim 1, characterized in that: The process of constructing a three-dimensional irregular stratum model with different undulating interfaces is as follows: According to the detection resolution requirements, the grid size is preset and a two-dimensional uniform initial model is formed accordingly; According to the preset number of strata and the height of each stratum, the two-dimensional uniform initial model is divided into a number of rectangular blocks in the vertical direction; wherein the number of rectangular blocks is the same as the number of strata; Configure the corresponding physical parameters for each rectangular block, copy the two-dimensional uniform initial model along the set direction, and obtain a corresponding three-dimensional array; For each replicated two-dimensional uniform initial model, the maximum disturbance value and the corresponding random distribution function of each interface are set according to the spatial correlation rule between the formation position and the fluctuation amplitude; A random disturbance value is added to the upper and lower interfaces of each rectangular block in the vertical direction. According to the position of the interface after disturbance, the attribute value of the corresponding position in the three-dimensional array is updated to form a three-dimensional irregular stratum model with different undulating interfaces.

6. The method for automatically modeling irregular interface structure and water body nesting according to claim 1, characterized in that: The process of constructing an irregular water-bearing structure model by combining multiple ellipsoids is as follows: Based on the spatial distribution density, morphological parameters, and connectivity characteristics of water-bearing structures, and in accordance with the spatial correlation rules between stratum position and water-bearing caves, the center coordinate positions, rotation angles, and semi-axis length ratios of multiple ellipsoids were randomly generated within a preset range on both sides of the fault in the area ahead of the tunnel face to simulate inclined water channels and beaded caves. According to the central coordinate positions of multiple ellipsoids and the corresponding rotation angles and semi-axis length ratios, the corresponding ellipsoid models are constructed, and these ellipsoid models are combined and remapped to the original coordinate system to form an irregular water-bearing structure model.

7. An automatic modeling system for irregular interface structure and water body nesting, characterized in that: include: The data processing module is used to obtain multi-source tunnel data, extract geological parameter interfaces, analyze water-bearing structural characteristics, and calibrate physical parameter thresholds. This module then constructs spatial association rules between stratum position and fluctuation amplitude, spatial association rules between stratum position and water-bearing caves, and a correspondence table between lithology, surrounding rock classification, and physical parameter parameters. The model construction module is used to construct a three-dimensional irregular stratum model with different undulating interfaces based on the spatial association rules between stratum position and undulation amplitude; based on the spatial association rules between stratum position and water-bearing caves and the correspondence table between lithology, surrounding rock classification and physical property parameters, an irregular water-bearing structure model is constructed by combining multiple ellipsoids; The model fusion module is used to directly map the physical parameters of the irregular water-bearing structure model to the corresponding positions of the 3D irregular stratum model under the same grid coordinates, forming a composite model with the undulation of the stratum interface and the nesting of water bodies; The process of constructing the spatial association rules between stratum position and fluctuation amplitude is as follows: Statistical analysis of stratum depth, interface fluctuation distribution and fluctuation amplitude from drill core data; At the fault location, the inverse distance weighted method of spatial interpolation is used to generate a continuous undulation surface, and then the probability density distribution of the stratum interface undulation and the undulation wavelength range are calculated, and the spatial correlation rules between stratum position and undulation amplitude are obtained.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps in the method for automatic modeling of irregular interface structure and water body nesting as described in any one of claims 1 to 6 are implemented.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps in the method for automatic modeling of irregular interface structure and water body nesting as described in any one of claims 1 to 6 are implemented.

Images