Method and system for automatic generation of parameterized three-dimensional geological and structural models for geotechnical engineering

Abstract

This invention discloses a method and system for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering, belonging to the field of information-based modeling technology for geotechnical engineering. The method includes the following steps: acquiring borehole data of the geotechnical engineering area and extracting corresponding geological parameters; applying parameter perturbations of a preset amplitude to the target geological parameters and generating candidate three-dimensional stratigraphic geometry results; determining the effective area of ​​influence of the target geological parameters based on the three-dimensional stratigraphic geometry results; establishing a monotonic constraint relationship between the change in target geological parameters and the change in stratigraphic geometry within the effective area, and generating a three-dimensional stratigraphic model; and combining the three-dimensional stratigraphic model with the engineering structural model to obtain a parametric three-dimensional geological and structural model. This invention solves the problems of unpredictable stratigraphic geometry and unstable structure-stratigraphic relationships caused by parameter changes in existing three-dimensional geological modeling by analyzing the stratigraphic geometric spatial response caused by geological parameter perturbations and introducing a monotonic constraint mechanism.

Description

Technical Field

[0001] This invention relates to the field of information modeling technology in geotechnical engineering, and more specifically, to a method and system for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering. Background Technology

[0002] With the increasing demand for digital and refined design in geotechnical engineering, parametric 3D geological and structural models have become an important foundation for engineering investigation and analysis, structural design, and construction decisions. Current technologies typically use borehole data as a basis, taking geological parameters such as stratum thickness, interface elevation, dip angle, and strike as input. A 3D geological model is then generated through methods such as surface interpolation, implicit functions, or voxel reconstruction, and further combined with the engineering structural model.

[0003] However, existing parametric modeling methods generally assume a stable correspondence between geological parameters and the generated 3D stratigraphic geometry, failing to adequately consider the impact of factors such as uneven borehole spatial distribution, variations in data density, and differences in interpolation algorithms on the model's geometric results. In practical applications, even with identical geological parameter values ​​or minor adjustments, the generated 3D stratigraphic geometry may still exhibit significant differences under different data conditions, resulting in a lack of a stable and monotonic response relationship between parameter changes and stratigraphic geometric changes.

[0004] The aforementioned problems make it difficult to predict the response of three-dimensional geological models to parameter adjustments. Engineers find it difficult to effectively control, quickly correct, or conduct multi-scheme comparative analysis based on the parameters, thus limiting the practicality and reliability of parametric three-dimensional geological and structural models in geotechnical engineering design and analysis.

[0005] The aforementioned disclosed technical solutions suffer from at least the following technical problems: Existing methods for automatically generating parametric 3D geological and structural models for geotechnical engineering typically use geological parameters such as stratum thickness, dip angle, and strike as input, generating corresponding 3D geological models through surface interpolation or implicit functions. However, under different borehole spatial distributions or data densities, the same set of geological parameters often generates 3D models with significantly different geometric shapes. This results in a lack of a stable and monotonic correspondence between parameter changes and the generated 3D geometric results, making the model's response to parameter adjustments unpredictable. Consequently, it is difficult to effectively control, quickly correct, or conduct comparative analysis of the model based on parameters, reducing the practicality of parametric modeling in geotechnical engineering design and analysis. Summary of the Invention

[0006] To overcome the aforementioned deficiencies of the prior art, embodiments of the present invention provide a method and system for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering. By analyzing the stratigraphic geometric spatial response caused by geological parameter disturbances and introducing a monotonic constraint mechanism, the method achieves controllable and consistent generation of stratigraphic models and engineering structural models, thus solving the problems of unpredictable stratigraphic geometry and unstable structure-stratigraphic relationships caused by parameter changes in existing three-dimensional geological modeling.

[0007] To achieve the above objectives, the present invention provides the following technical solution: On the one hand, the automatic generation method for parametric three-dimensional geological and structural models for geotechnical engineering includes the following steps: acquiring borehole data of the geotechnical engineering area and extracting the corresponding geological parameters; applying parameter perturbation of a preset amplitude to the target geological parameters and generating corresponding candidate three-dimensional stratigraphic geometry results; based on the three-dimensional stratigraphic geometry results, calculating the spatial response relationship between the stratigraphic geometry change and the target geological parameter change, and determining the effective area of ​​the target geological parameter; within the effective area, establishing a monotonic constraint relationship between the target geological parameter change and the stratigraphic geometry change, and generating a three-dimensional stratigraphic model; and combining the three-dimensional stratigraphic model with the engineering structural model to obtain a parametric three-dimensional geological and structural model.

[0008] In a preferred embodiment, the step of acquiring borehole data in the geotechnical engineering area and extracting corresponding geological parameters includes: analyzing the borehole data in the geotechnical engineering area to obtain the spatial location and stratification record of each borehole; identifying the boundary position between different strata in each borehole based on the stratification record; extracting geological parameters describing the spatial morphology of the stratum by using the boundary position of the same stratum in different boreholes as the correlation object, and establishing the correlation between the geological parameters and the spatial location of the borehole.

[0009] In a preferred embodiment, identifying the boundary locations between different strata in each borehole based on the layered records includes: sorting the layered records of a single borehole along the borehole depth direction to form a continuous stratigraphic sequence; comparing adjacent stratigraphic records one by one in the stratigraphic sequence to identify locations where stratigraphic properties change; determining the borehole depth value corresponding to the location where stratigraphic properties change as a candidate stratigraphic boundary point; verifying the validity of the stratigraphic segments corresponding to the candidate stratigraphic boundary points, merging segments that do not meet the validity conditions, and correcting the corresponding candidate stratigraphic boundary points accordingly; and determining the corrected candidate stratigraphic boundary points as the stratigraphic boundary points in that borehole.

[0010] In a preferred embodiment, the step of applying a preset amplitude of parameter perturbation to the target geological parameter and generating corresponding candidate three-dimensional stratigraphic geometry results includes: determining the baseline parameter distribution of the target geological parameter at each borehole based on the original values ​​of the target geological parameter in the borehole data; applying a uniform perturbation coefficient to the baseline parameter distribution while keeping other geological parameters unchanged to obtain a forward perturbation parameter distribution and a reverse perturbation parameter distribution respectively; correcting the position of the stratigraphic interface control points corresponding to each borehole according to the forward perturbation parameter distribution and the reverse perturbation parameter distribution to form a perturbed stratigraphic interface control point set; and generating corresponding candidate three-dimensional stratigraphic geometry results based on the perturbed stratigraphic interface control point set while keeping the stratigraphic topology unchanged.

[0011] In a preferred embodiment, determining the baseline parameter distribution of the target geological parameter at each borehole based on the original values ​​of the target geological parameter in the borehole data includes: analyzing the geological records of each borehole within the geotechnical engineering area and extracting the parameter values ​​corresponding to the target geological parameter from each borehole; associating the parameter values ​​of the target geological parameter with the spatial location of the corresponding borehole to form a parameter set organized according to the spatial location of the borehole; verifying the validity of the parameter set and eliminating abnormal values ​​that do not meet the engineering rationality; and constructing the baseline parameter distribution of the target geological parameter at each borehole based on the verified parameter set.

[0012] In a preferred embodiment, the step of calculating the spatial response relationship between the stratigraphic geometric change and the change in the target geological parameter based on the three-dimensional stratigraphic geometry results includes: obtaining the three-dimensional stratigraphic geometry results corresponding to the baseline value of the target geological parameter, and the candidate three-dimensional stratigraphic geometry results corresponding to the positive and negative perturbations of the target geological parameter; spatially aligning the baseline three-dimensional stratigraphic geometry results and the candidate three-dimensional stratigraphic geometry results to establish a unified spatial sampling baseline; calculating the stratigraphic geometric change of the candidate three-dimensional stratigraphic geometry results relative to the baseline three-dimensional stratigraphic geometry results at each spatial location under the unified spatial sampling baseline; normalizing the stratigraphic geometric change based on the parameter change corresponding to the positive and negative perturbations of the target geological parameter to obtain the response value of the stratigraphic geometric change to the change in the target geological parameter at each spatial location; and organizing the response values ​​at each spatial location according to spatial location to form the spatial response distribution of the target geological parameter.

[0013] In a preferred embodiment, determining the effective area of ​​action of the target geological parameter includes: acquiring the stratigraphic geometric response of each spatial location under positive and negative disturbances of the target geological parameter based on the spatial response distribution of the target geological parameter; determining the stability of the stratigraphic geometry response to changes in the target geological parameter at each spatial location based on the symmetry index between the positive and negative disturbance responses; identifying spatial locations within the spatial range where the response stability satisfies the opposite directions of positive and negative disturbances and the amplitude of change is consistent, as stable response locations; performing spatial connectivity analysis on the stable response locations to extract continuous spatial regions that satisfy stratigraphic continuity constraints; and determining the continuous spatial regions as the effective area of ​​action of the target geological parameter.

[0014] In a preferred embodiment, establishing a monotonic constraint relationship between the change in target geological parameters and the change in stratigraphic geometry within the effective area includes: determining the correspondence between the direction of change of target geological parameters and the direction of change of stratigraphic geometry within the effective area based on the stratigraphic geometric response results under positive and negative disturbances of target geological parameters; identifying a stable correspondence between the direction of parameter change and the direction of change of stratigraphic geometry within the effective area as a monotonic correspondence between target geological parameters and stratigraphic geometry; and constraining the direction of change of stratigraphic geometry caused by the change in target geological parameters during subsequent stratigraphic geometry generation or parameter update, allowing only the direction to be consistent with the monotonic correspondence.

[0015] In a preferred embodiment, generating a three-dimensional stratigraphic model includes: within the effective action area, acquiring the change in target geological parameters that satisfies the monotonic constraint relationship; based on the change in target geological parameters, adjusting the positions of stratigraphic interface control points within the effective action area so that the direction of position adjustment conforms to the direction of stratigraphic geometric change defined by the monotonic constraint relationship; after completing the position adjustment of the stratigraphic interface control points, reconstructing the stratigraphic interface based on the adjusted stratigraphic interface control points to form a continuous three-dimensional stratigraphic interface; and generating a three-dimensional stratigraphic model that satisfies the monotonic constraint relationship based on the continuous three-dimensional stratigraphic interface.

[0016] On the other hand, the automatic generation system for parametric 3D geological and structural models for geotechnical engineering includes the following modules: a parameter acquisition module for acquiring borehole data of the geotechnical engineering area and extracting corresponding geological parameters; a candidate stratum generation module for applying parameter perturbations of a preset amplitude to the target geological parameters and generating corresponding candidate 3D stratum geometry results; an effective action area identification module for calculating the spatial response relationship between the stratum geometry change and the target geological parameter change based on the 3D stratum geometry results and determining the effective action area of ​​the target geological parameters; a monotonic constraint modeling module for establishing a monotonic constraint relationship between the target geological parameter change and the stratum geometry change within the effective action area and generating a 3D stratum model; and a model combination module for combining the 3D stratum model with the engineering structure model to obtain a parametric 3D geological and structural model.

[0017] The technical effects and advantages of the present invention regarding the automatic generation method and system for parametric three-dimensional geological and structural models in geotechnical engineering are as follows: This invention introduces a constraint relationship between changes in target geological parameters and changes in stratigraphic geometry, enabling the three-dimensional stratigraphic model to maintain consistency and stability in the direction of geometric changes during parameter adjustment, thereby avoiding unreasonable fluctuations in stratigraphic morphology as parameters change. Based on this, the generated three-dimensional stratigraphic model can be effectively combined with engineering structural models to form parameterized geological and structural models, improving the controllability of geotechnical engineering modeling and its engineering adaptability to parameter changes. Attached Figure Description

[0018] Figure 1 This is a flowchart illustrating the automatic generation method for parametric three-dimensional geological and structural models in geotechnical engineering according to the present invention. Figure 2 This is a schematic diagram of the parametric three-dimensional geological and structural model automatic generation system for geotechnical engineering according to the present invention. Detailed Implementation

[0019] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0020] Example 1, Figure 1 This invention presents a method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering, comprising the following steps: S1. Obtain borehole data for the geotechnical engineering area and extract the corresponding geological parameters; In this embodiment, acquiring borehole data of the geotechnical engineering area and extracting corresponding geological parameters includes: The original borehole data within the geotechnical engineering area is parsed in a unified format to obtain the spatial coordinates, depth, and layer records of each borehole. Based on borehole stratification records, the corresponding stratigraphic boundary points in each borehole are identified according to the depth of the stratigraphic interface. Using the boundary points of the same stratum in different boreholes as the associated objects, geological parameters describing the spatial morphology of the stratum are extracted. The geological parameters include at least one or more of the following: stratum interface elevation, stratum thickness, dip angle, and dip angle. The geological parameters are then linked to the spatial location of the corresponding borehole to form a parameterized geological input set for subsequent parameter disturbance and spatial response analysis.

[0021] Furthermore, the step of identifying the corresponding formation boundary points in each borehole based on the borehole layering record and according to the depth of formation interface occurrence includes: The stratigraphic records of a single borehole are sorted according to the borehole depth direction to form a continuous stratigraphic sequence along the borehole depth direction; In a continuous stratigraphic sequence, traverse adjacent stratigraphic records to determine whether the lithological class or engineering classification of two adjacent strata has changed. When the lithology or engineering classification of two adjacent strata changes, the borehole depth value corresponding to the location of the change is extracted as a candidate stratum boundary point. Stratigraphic segments with a thickness less than the preset minimum effective thickness are merged, and the positions of the candidate stratigraphic boundaries are corrected accordingly. The revised candidate stratigraphic boundary point is determined as the stratigraphic boundary point in the borehole.

[0022] In a continuous stratigraphic sequence, traversing adjacent stratigraphic records to determine whether the lithological type or engineering classification of two adjacent strata has changed specifically involves: Compare the lithological codes or engineering classification identifiers corresponding to two adjacent strata; When the lithology codes are different or the engineering classification identifiers are inconsistent, it is determined that there is a stratigraphic interface change between the two adjacent strata.

[0023] S2. Apply parameter perturbation of preset amplitude to the target geological parameters and generate corresponding candidate three-dimensional stratigraphic geometry results respectively; In this embodiment, for each target geological parameter, the parameter perturbation, spatial response analysis, effective area determination, and monotonic constraint relationship establishment processes are performed independently.

[0024] In this embodiment, applying a preset amplitude of parameter perturbation to the target geological parameters and generating corresponding candidate three-dimensional stratigraphic geometric results includes: For at least one target geological parameter, based on the original value of the target geological parameter in the borehole data, determine its baseline parameter distribution at each borehole. While keeping other geological parameters unchanged, a disturbance coefficient is uniformly applied to the baseline parameter distribution of the target geological parameters at each borehole to obtain the forward disturbance parameter distribution and the reverse disturbance parameter distribution. The disturbance coefficient is limited to a preset disturbance amplitude range and is kept consistent at each borehole to eliminate the influence of uneven borehole distribution on geometric response analysis. Based on the forward and reverse disturbance parameter distributions, the positions of the formation interface control points corresponding to each borehole are corrected to form a disturbed set of formation interface control points. Without altering the stratigraphic topology and interpolation rules, the perturbed stratigraphic interface control point set is used to generate corresponding candidate three-dimensional stratigraphic geometric results. The candidate three-dimensional stratigraphic geometry results are then stored in a one-to-one correspondence with the three-dimensional stratigraphic geometry results generated by the baseline parameter distribution for subsequent calculation of stratigraphic geometry changes.

[0025] It should be noted that determining the baseline parameter distribution at each borehole based on the original values ​​of the target geological parameter in the borehole data includes: The original geological records of each borehole in the geotechnical engineering area are analyzed, and the parameter values ​​corresponding to the target geological parameters in each borehole are extracted. The values ​​of the target geological parameters are associated with the spatial coordinates of the corresponding boreholes to form a set of parameter values ​​indexed by the spatial location of the boreholes. The set of parameter values ​​is checked for consistency, and missing values ​​or abnormal values ​​that exceed the reasonable engineering range are removed. The verified parameter values ​​are organized according to the spatial location of the borehole to construct the baseline parameter distribution of the target geological parameters at each borehole.

[0026] The disturbance coefficient is determined based on the statistical discrete characteristics of the target geological parameters in the borehole data, specifically: Calculate the range or standard deviation of the target geological parameters at each borehole; The disturbance coefficient is limited to a value that does not exceed the range of values ​​or a preset proportion of the standard deviation.

[0027] The process involves uniformly applying a disturbance coefficient to the baseline parameter distribution of the target geological parameters at each borehole to obtain the positive and negative disturbance parameter distributions, including: Based on the aforementioned baseline parameter distribution, the baseline parameter values ​​for the target geological parameters at each borehole are determined. Based on the statistical discrete characteristics of the target geological parameters in the borehole data, determine the perturbation coefficients used for parameter perturbation; Using the same disturbance coefficient, the reference parameter values ​​at each borehole are subjected to forward and reverse disturbances respectively to obtain the corresponding forward and reverse disturbance parameter values; The values ​​of the forward and reverse disturbance parameters at each borehole are organized according to the spatial location of the borehole, forming forward disturbance parameter distributions and reverse disturbance parameter distributions respectively.

[0028] Based on the disturbed target geological parameters, the positions of the corresponding stratigraphic interface control points at each borehole are corrected according to the changes in the target geological parameters, specifically as follows: Convert the changes in the target geological parameters into the corresponding interface elevation offsets; The position of the formation interface control point is then shifted and corrected in the vertical direction of the borehole.

[0029] S3. Based on the three-dimensional stratigraphic geometry results, calculate the spatial response relationship between the stratigraphic geometric changes and the changes in the target geological parameters, and determine the effective area of ​​the target geological parameters. The effective area can be understood as the spatial region in which changes in target geological parameters have a dominant influence on the geomorphology of the strata.

[0030] In this embodiment, calculating the spatial response relationship between the stratigraphic geometric change and the target geological parameter change based on the three-dimensional stratigraphic geometry results includes: Obtain the three-dimensional stratigraphic geometry results corresponding to the benchmark geological parameters, as well as the candidate three-dimensional stratigraphic geometry results corresponding to the positive and negative perturbation parameter distributions; The baseline three-dimensional stratigraphic geometry results are spatially aligned with each candidate three-dimensional stratigraphic geometry result to establish a unified spatial sampling coordinate system; Under the unified spatial sampling coordinate system, the stratigraphic geometric change of the candidate three-dimensional stratigraphic geometric result relative to the reference three-dimensional stratigraphic geometric result is calculated for each spatial sampling position; the spatial sampling position is a discrete spatial sampling point generated within the spatial range of the three-dimensional stratigraphic model according to a preset spatial resolution or stratigraphic structural features.

[0031] Based on the changes in the positive and negative disturbance parameters, the response value of the formation geometric change at each spatial sampling location relative to the change in the target geological parameter is calculated; the response value is the ratio of the formation geometric change to the change in the target geological parameter or a normalized function value. The response values ​​at each spatial sampling location are organized according to their spatial location to form the spatial response distribution of the target geological parameter.

[0032] By calculating the spatial response relationship between the changes in stratigraphic geometry and the changes in target geological parameters based on candidate three-dimensional stratigraphic geometry results, this invention can quantitatively characterize the actual influence of changes in target geological parameters on stratigraphic geometry on a spatial scale, thereby transforming the originally unpredictable and uncontrollable parameter-geometric relationship into a response field with clear spatial distribution characteristics.

[0033] Furthermore, by comparing and calculating the three-dimensional stratigraphic geometry results corresponding to forward and reverse perturbations under a unified spatial sampling coordinate system, the interference introduced by differences in model discretization methods or local interpolation instability can be effectively reduced, resulting in better stability and repeatability of the spatial response distribution.

[0034] Based on the spatial response distribution, the effective influence of target geological parameters in different spatial regions can be accurately identified, thus providing a reliable data foundation for determining the effective influence area of ​​target geological parameters and establishing a stable correspondence between parameter changes and stratigraphic geometric changes, thereby improving the controllability and engineering applicability of parametric three-dimensional geological models.

[0035] The effective area for determining the target geological parameters includes: Based on the spatial response distribution of the target geological parameters, the stratigraphic geometric changes at each spatial location under positive and negative parameter perturbation conditions are obtained. For each spatial location, a symmetry index for the positive and negative perturbation responses is constructed to characterize the stability of the response of the stratigraphic geometry to changes in the target geological parameters at that spatial location. Based on the aforementioned symmetry index, stable response locations that satisfy the opposite directions of positive and negative disturbances and have consistent amplitude changes can be identified within the spatial range. In the stable response locations, the response intensity is further normalized based on the intensity of the stratigraphic geometric change caused by the unit parameter change at each spatial location, and spatial locations with significant response characteristics to the change of target geological parameters are screened out. Spatial connectivity analysis is performed on the selected spatial locations to extract continuous spatial regions that satisfy the stratigraphic continuity constraint; The continuous spatial region is defined as the effective area of ​​the target geological parameters.

[0036] The specific formula for calculating the symmetry index is as follows:

[0037] in, This is an index of symmetry at spatial location x. For positive disturbance response, This is the response to the reverse perturbation.

[0038] when A value of 0 indicates that at this spatial location, the stratigraphic geometric change caused by a positive change in the target geological parameter is opposite in direction to the stratigraphic geometric change caused by a negative change in the parameter, and the two changes have good consistency in amplitude.

[0039] This indicates that the stratigraphic geometry at this spatial location is primarily controlled by the target geological parameter; the stratigraphic geometry's response to this parameter exhibits good stability and reversibility; and a clear and predictable correspondence exists between parameter changes and geometric changes. Therefore, when When the value is 0, this spatial location can be regarded as the stable location of the target geological parameters.

[0040] when When the value is 1, it indicates that at this spatial location: when the target geological parameter changes in the positive and negative directions, the resulting stratigraphic geometric changes are in the same direction, or one of the changes produces almost no geometric response; the stratigraphic geometric changes do not reverse accordingly with the direction of parameter changes.

[0041] This indicates that the stratigraphic geometry at this spatial location is not stably controlled by the target geological parameters; the stratigraphic geometry may be primarily influenced by other geological parameters, geometric constraints, or numerical interpolation factors; and there is no stable, reversible correspondence between parameter changes and geometric changes. Therefore, when When the value is 1, the spatial location is not within the effective range of the target geological parameter.

[0042] when A value between 0 and 1 indicates that at this spatial location, both positive and negative changes in the target geological parameters will cause changes in the stratigraphic geometry; however, the two changes are not completely symmetrical in terms of amplitude or direction.

[0043] This indicates that the spatial location is somewhat sensitive to the target geological parameters; however, this sensitivity is affected by the coupling of other geological parameters, spatial constraints, or uncertainties in local data; and the correspondence between parameter changes and geometric changes is not stable enough. Therefore, this spatial location should not be directly used as the stable area of ​​the target geological parameters, but needs to be further screened in conjunction with other criteria.

[0044] S4. Within the effective area, establish a monotonic constraint relationship between the changes in target geological parameters and the changes in stratigraphic geometry, and generate a three-dimensional stratigraphic model. In this embodiment, establishing a monotonic constraint relationship between the change in target geological parameters and the change in stratigraphic geometry within the effective area includes: Within the effective area, the directions of the stratigraphic geometric changes generated when the target geological parameters undergo positive and negative changes are recorded respectively. Based on the direction of the stratigraphic geometric change, determine the consistency relationship between the direction of the change of the target geological parameters and the direction of the stratigraphic geometric change. The direction of parameter change that satisfies the consistency relationship within the effective area is determined to be a monotonic correspondence with the direction of formation geometric change. In subsequent stratigraphic geometry generation or parameter update processes, the direction of stratigraphic geometry changes caused by changes in target geological parameters is constrained, and only those consistent with the monotonic correspondence are allowed.

[0045] It should be noted that determining the consistency relationship between the direction of change of the target geological parameters and the direction of change of the stratigraphic geometry, based on the direction of the stratigraphic geometry change, includes: Within the effective area, the results of stratigraphic geometric changes corresponding to positive and negative changes in the target geological parameters are obtained respectively; By comparing the results of formation geometric changes under the positive and negative change conditions, it is determined whether the directions of formation geometric changes are opposite in the two cases; If, within the effective area, the direction of change of the stratigraphic geometry under positive change conditions is opposite to the direction of change under negative change conditions, then it is determined that there is a consistency relationship between the direction of change of the target geological parameters and the direction of change of the stratigraphic geometry. If, within the effective area, the direction of change of the stratigraphic geometry under both positive and negative change conditions is the same or unclear, then it is determined that there is no consistent relationship between the direction of change of the target geological parameters and the direction of change of the stratigraphic geometry.

[0046] Furthermore, the generation of the three-dimensional stratigraphic model includes: Within the effective area, the change in the target geological parameters that satisfy the monotonic constraint relationship is obtained; Based on the change in the target geological parameters, the position of the stratigraphic interface control points in the effective area is updated so that the movement direction of the stratigraphic interface control points is consistent with the direction of stratigraphic geometric change defined by the monotonic constraint relationship. After updating the positions of the stratigraphic interface control points, spatial interpolation or surface reconstruction is performed on the updated stratigraphic interface control points to generate a continuous three-dimensional stratigraphic interface. Based on the three-dimensional stratigraphic interface, a three-dimensional stratigraphic model that satisfies the monotonic constraint relationship is generated.

[0047] Furthermore, after updating the positions of the stratigraphic interface control points, the updated stratigraphic interface control points are subjected to spatial interpolation or surface reconstruction to generate a continuous three-dimensional stratigraphic interface, including: The updated stratigraphic interface control points are classified according to the corresponding stratigraphic interfaces, and control point sets for each stratigraphic interface are constructed. Based on each set of control points, spatial adjacency relationships between control points are established in three-dimensional space to form a spatial constraint structure for surface generation; Under the condition of satisfying the spatial adjacency relationship of control points, spatial interpolation or surface fitting is performed on each set of control points to generate the corresponding continuous stratigraphic interface surface. The continuity of the generated stratigraphic interface surface is checked to ensure that there are no intersections, breaks or discontinuities between adjacent stratigraphic interfaces. The stratigraphic interface surface that has passed the continuity check is used as the stratigraphic interface in the three-dimensional stratigraphic model.

[0048] S5. Combine the three-dimensional stratigraphic model with the engineering structural model to obtain a parametric three-dimensional geological and structural model.

[0049] The engineering structure model is a three-dimensional engineering structure model constructed based on engineering design data or existing structural models.

[0050] In this embodiment, the combination of the three-dimensional stratigraphic model and the engineering structure model to obtain a parametric three-dimensional geological and structural model includes: Obtain the engineering structure model corresponding to the three-dimensional stratigraphic model, and import the engineering structure model into a spatial coordinate system consistent with the three-dimensional stratigraphic model; Under the condition that the engineering structure model and the three-dimensional stratigraphic model are in a unified spatial coordinate system, the engineering structure model is spatially positioned and mapped to the spatial range of the three-dimensional stratigraphic model; Based on the spatial positioning results, the relative positional relationship between the engineering structure model and various stratigraphic interfaces and stratigraphic units in space is identified, and the spatial interaction area between the engineering structure and different stratigraphic units is determined. Within the spatial interaction area, the consistency of the contact relationship between the engineering structure model and the stratigraphic interface is checked based on the monotonic constraint relationship of the target geological parameters. If the contact relationship between the engineering structure model and the stratum interface is detected to not satisfy the monotonic constraint relationship, then the stratum geometry or the spatial position of the engineering structure in the corresponding spatial interaction area will be subject to restricted adjustment. After completing the consistency check and constraint adjustment, the three-dimensional stratigraphic model and the engineering structure model are combined to generate a three-dimensional geological and structural model that meets the parameterized constraints.

[0051] The spatial positioning of the engineering structure model, mapping the engineering structure model to the spatial range of the three-dimensional geological model, includes: After importing the engineering structure model into a spatial coordinate system consistent with the three-dimensional geological model, the initial spatial position of the engineering structure model in the spatial coordinate system is obtained; Obtain the overall spatial extent of the three-dimensional stratigraphic model in the spatial coordinate system; Based on the initial spatial position of the engineering structure model and the overall spatial range of the three-dimensional geological model, the engineering structure model is translated and / or rotated so that the spatial position of the engineering structure model corresponds to the spatial range of the three-dimensional geological model. After completing the translation and / or rotation, the spatial positioning result of the engineering structure model within the spatial range of the three-dimensional geological model is determined.

[0052] Based on the spatial positioning results, the process of identifying the relative spatial positions of the engineering structure model with various stratigraphic interfaces and stratigraphic units, and determining the spatial interaction areas between the engineering structure and different stratigraphic units, includes: Based on the spatial positioning results, the actual spatial range of each structural unit in the engineering structural model in three-dimensional space is obtained; Based on the three-dimensional stratigraphic model, the spatial extent of each stratigraphic unit enclosed by adjacent stratigraphic interfaces is obtained; The spatial range of the structural unit is compared with the spatial range of each stratigraphic unit to determine whether the structural unit overlaps with the corresponding stratigraphic unit in space. For structural units and stratigraphic units that overlap in space, the corresponding overlapping spatial region is defined as the spatial interaction region between the engineering structure and the stratigraphic unit.

[0053] The restricted adjustment of the spatial position of the stratigraphic geometry or engineering structure within the corresponding spatial interaction area includes: Within the spatial interaction area, identify local regions that cause the contact relationship between the engineering structure model and the stratum interface to not satisfy the monotonic constraint relationship; Based on the monotonic constraint relationship, the allowable direction of geological geometric change or the direction of spatial position adjustment of engineering structure within the local area are determined. While keeping the spatial position of the strata geometry and engineering structure outside the effective area unchanged, only the position of the strata interface control points or related parts of the engineering structure in the local area is corrected along the allowed adjustment direction; After completing the position correction, the contact relationship between the engineering structure model and the geological interface is rechecked until the contact relationship satisfies the monotonic constraint relationship.

[0054] Through the above technical solution, this invention transforms the uncertainties of geological parameters that are difficult to quantify and control in traditional 3D geological modeling into parameter-geometric response relationships with clear spatial distribution characteristics and monotonic constraint relationships. This enables the predictable and controllable generation of 3D stratigraphic models under parameter perturbation conditions. Furthermore, it achieves consistent coupling between engineering structure models and stratigraphic models, effectively avoiding distortion of structure-stratigraphic relationships or local geometric conflicts, and improving the reliability and engineering applicability of parametric 3D geological and structural models in geotechnical engineering analysis and design.

[0055] The method described in this invention is not limited to a single stratum or a single engineering structure type, and can be applied to multi-layered geological structures and different types of engineering structure models according to actual engineering needs.

[0056] Example 2, Figure 2 The present invention provides a system for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering, characterized by comprising the following modules: Parameter acquisition module: used to acquire borehole data of the geotechnical engineering area and extract the corresponding geological parameters; Candidate stratigraphic generation module: used to apply parameter perturbation of preset amplitude to the target geological parameters and generate corresponding candidate three-dimensional stratigraphic geometric results; Effective Area Identification Module: Based on three-dimensional stratigraphic geometry results, this module calculates the spatial response relationship between stratigraphic geometric changes and target geological parameter changes, and determines the effective area of ​​action of the target geological parameters. Monotonic constraint modeling module: Used to establish a monotonic constraint relationship between the change of target geological parameters and the change of stratigraphic geometry within the effective area, and generate a three-dimensional stratigraphic model; Model Combination Module: Used to combine 3D geological models with engineering structure models to obtain parametric 3D geological and structural models.

[0057] The above formulas are all dimensionless calculations. The formulas are derived from software simulations based on a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation.

[0058] The above embodiments can be implemented, in whole or in part, by software, hardware, firmware, or any other combination thereof. When implemented using software, the above embodiments can be implemented, in whole or in part, in the form of a computer program product.

[0059] Those skilled in the art will recognize that the modules and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0060] In addition, the functional modules in the various embodiments of this application can be integrated into one processing module, or each module can exist physically separately, or two or more modules can be integrated into one module.

[0061] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

[0062] In conclusion, the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering, characterized in that, Includes the following steps: Obtain borehole data for the geotechnical engineering area and extract the corresponding geological parameters; Apply parameter perturbations of a preset amplitude to the target geological parameters and generate corresponding candidate three-dimensional stratigraphic geometry results; Based on the three-dimensional stratigraphic geometry results, the spatial response relationship between stratigraphic geometric changes and target geological parameter changes is calculated, and the effective area of ​​action of the target geological parameters is determined. Within the effective area, a monotonic constraint relationship is established between the changes in target geological parameters and the changes in stratigraphic geometry, and a three-dimensional stratigraphic model is generated. By combining the three-dimensional stratigraphic model with the engineering structural model, a parametric three-dimensional geological and structural model is obtained.

2. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 1, characterized in that, The acquisition of borehole data for the geotechnical engineering area and the extraction of corresponding geological parameters include: The borehole data within the geotechnical engineering area are analyzed to obtain the spatial location and stratification records of each borehole. Based on the layered records, the boundary locations between different formations in each borehole are identified; Using the boundary locations of the same stratum in different boreholes as the correlation objects, geological parameters describing the spatial morphology of the stratum are extracted, and the correlation between the geological parameters and the spatial location of the borehole is established.

3. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 2, characterized in that, The method of identifying the boundary locations between different formations in each borehole based on layered records includes: The stratigraphic records of a single borehole are sorted along the borehole depth direction to form a continuous stratigraphic sequence; In the stratigraphic sequence, adjacent stratigraphic records are compared one by one to identify the locations where stratigraphic properties change; The borehole depth value corresponding to the location where the stratigraphic properties change is determined as the candidate stratigraphic boundary point; The validity of the stratigraphic segments corresponding to the candidate stratigraphic boundaries is checked, and segments that do not meet the validity conditions are merged, and the corresponding candidate stratigraphic boundaries are corrected accordingly. The revised candidate stratigraphic boundary point is determined as the stratigraphic boundary point in the borehole.

4. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 3, characterized in that, The process of applying a preset amplitude of parameter perturbation to the target geological parameters and generating corresponding candidate three-dimensional stratigraphic geometric results includes: Based on the original values ​​of the target geological parameter in the borehole data, determine the baseline parameter distribution of the target geological parameter at each borehole. While keeping other geological parameters unchanged, a uniform perturbation coefficient is applied to the baseline parameter distribution to obtain the positive perturbation parameter distribution and the negative perturbation parameter distribution, respectively. Based on the distribution of the forward disturbance parameters and the distribution of the reverse disturbance parameters, the positions of the formation interface control points corresponding to each borehole are corrected to form a set of disturbed formation interface control points. While keeping the stratigraphic topology unchanged, corresponding candidate three-dimensional stratigraphic geometry results are generated based on the perturbed stratigraphic interface control point set.

5. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 4, characterized in that, The determination of the baseline parameter distribution at each borehole based on the original values ​​of the target geological parameter in the borehole data includes: Analyze the geological records of each borehole in the geotechnical engineering area and extract the parameter values ​​corresponding to the target geological parameters from each borehole; The parameter values ​​of the target geological parameters are associated with the spatial location of the corresponding borehole to form a parameter set organized according to the spatial location of the borehole. The validity of the parameter set is checked, and abnormal values ​​that do not meet the engineering rationality are eliminated; Based on the verified parameter set, the baseline parameter distribution of the target geological parameters at each borehole is constructed.

6. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 5, characterized in that, The calculation of the spatial response relationship between stratigraphic geometric changes and target geological parameter changes based on three-dimensional stratigraphic geometry results includes: Obtain the three-dimensional stratigraphic geometry corresponding to the baseline values ​​of the target geological parameter, as well as the candidate three-dimensional stratigraphic geometry corresponding to the forward and reverse perturbations of the target geological parameter; Spatially align the baseline 3D stratigraphic geometry with the candidate 3D stratigraphic geometry to establish a unified spatial sampling baseline; Under the unified spatial sampling benchmark, the stratigraphic geometric change of the candidate three-dimensional stratigraphic geometric results at each spatial location relative to the benchmark three-dimensional stratigraphic geometric results is calculated; Based on the parameter changes corresponding to the positive and negative disturbances of the target geological parameters, the stratigraphic geometric changes are normalized to obtain the response values ​​of stratigraphic geometric changes to changes in target geological parameters at each spatial location. The response values ​​at each spatial location are organized according to spatial location to form the spatial response distribution of the target geological parameters.

7. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 6, characterized in that, The effective area for determining the target geological parameters includes: Based on the spatial response distribution of the target geological parameters, the stratigraphic geometric response of each spatial location under the conditions of positive and negative disturbance of the target geological parameters is obtained. For each spatial location, the stability of the stratigraphic geometry at that location in response to changes in target geological parameters is determined based on the symmetry index between the forward and reverse disturbance responses. Within the spatial range, identify spatial locations where the response stability satisfies the condition that the positive and negative disturbances have opposite directions and consistent amplitudes, and use these locations as stable response locations. Spatial connectivity analysis is performed on the stable response locations to extract continuous spatial regions that satisfy the formation continuity constraint; The continuous spatial region is defined as the effective area of ​​the target geological parameters.

8. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 7, characterized in that, The establishment of a monotonic constraint relationship between the changes in target geological parameters and the changes in stratigraphic geometry within the effective area includes: Within the effective area, based on the stratigraphic geometric response results under positive and negative disturbance conditions of the target geological parameters, the correspondence between the direction of change of the target geological parameters and the direction of change of stratigraphic geometry is determined. Within the effective area, the direction of parameter change and the direction of stratigraphic geometry change will be stably correlated, which will be identified as the monotonic correspondence between the target geological parameters and stratigraphic geometry. In subsequent stratigraphic geometry generation or parameter update processes, the direction of stratigraphic geometry changes caused by changes in target geological parameters is constrained, and only those consistent with monotonic correspondences are allowed.

9. The method for automatically generating parametric three-dimensional geological and structural models for geotechnical engineering according to claim 8, characterized in that, The generation of the three-dimensional stratigraphic model includes: Within the effective area of ​​action, obtain the change in target geological parameters that satisfy the monotonic constraint relationship; Based on the changes in the target geological parameters, the positions of the stratigraphic interface control points within the effective area are adjusted so that the direction of the position adjustment conforms to the direction of stratigraphic geometric change defined by the monotonic constraint relationship. After adjusting the positions of the stratigraphic interface control points, the stratigraphic interface is reconstructed based on the adjusted stratigraphic interface control points to form a continuous three-dimensional stratigraphic interface. Based on the continuous three-dimensional stratigraphic interface, a three-dimensional stratigraphic model that satisfies the monotonic constraint relationship is generated.

10. A system using the automatic generation method for parametric three-dimensional geological and structural models for geotechnical engineering as described in any one of claims 1-9, characterized in that, Includes the following modules: Parameter acquisition module: used to acquire borehole data of the geotechnical engineering area and extract the corresponding geological parameters; Candidate stratigraphic generation module: used to apply parameter perturbation of preset amplitude to the target geological parameters and generate corresponding candidate three-dimensional stratigraphic geometric results; Effective Area Identification Module: Based on three-dimensional stratigraphic geometry results, this module calculates the spatial response relationship between stratigraphic geometric changes and target geological parameter changes, and determines the effective area of ​​action of the target geological parameters. Monotonic constraint modeling module: Used to establish a monotonic constraint relationship between the change of target geological parameters and the change of stratigraphic geometry within the effective area, and generate a three-dimensional stratigraphic model; Model Combination Module: Used to combine 3D geological models with engineering structure models to obtain parametric 3D geological and structural models.

Images