Method and system for converting geological model into seepage model

Abstract

The invention discloses a method and system for converting a geologic model into a seepage model, and relates to the technical field of geologic model conversion, and the method comprises the steps: building a research region three-dimensional terrain model and a research region three-dimensional geologic model; discretizing the three-dimensional terrain model and the three-dimensional geological model of the research area through data conversion, and layering according to different stratums to form an ASCII (American Standard Code for Information Interchange) file; uniformly discretized node coordinate data obtained after conversion is completed are imported into seepage numerical calculation in a layered mode, a hydrogeological model is constructed, and a seepage numerical calculation model is constructed by assigning hydrogeological parameters and source and sink items. It is ensured that the three-dimensional geologic model conforms to the actual situation, and the simulation matching degree is high; data conversion between a geologic model and an underground water numerical model is achieved through data conversion of the technology, and the defects of international mainstream software in the aspect of pretreatment work of mountainous area complex terrains, stratum lithology and the like are overcome; a large amount of pretreatment workload is reduced, and the efficiency and precision of mountain underground water numerical model construction are improved.

Description

Technical Field

[0001] This invention relates to the field of geological model conversion technology, and in particular to a geological model conversion method and system for seepage flow model. Background Technology

[0002] When performing numerical simulations of groundwater flow in mountainous areas, mainstream international groundwater numerical simulation software such as GMS, VisualModlow, and FEFLOW require importing topographic and stratigraphic lithology data into the model layer by layer according to the software's required file format to construct a 3D geological model structure. This has the following drawbacks: In complex terrain and deep canyon areas with nearly vertical rock strata, constructing a groundwater numerical model requires importing topographic and stratigraphic lithology data layer by layer according to the software's required format. This results in an exceptionally large amount of basic data input, increasing the difficulty of model construction and leading to long processing times and low efficiency. Currently, topographic and stratigraphic lithology data acquired during the survey and design phase of large-scale water conservancy and hydropower projects are mostly displayed in the form of 3D BIM models. There is a lack of readily available interactive conversion interfaces between BIM software and groundwater numerical simulation software, preventing direct conversion and resulting in a massive amount of data processing. Furthermore, during the construction of the groundwater numerical model, some strata need to be simplified to reduce the number of meshes in the model, leading to lower accuracy of the 3D geological structure model.

[0003] This invention leverages the advantages of BIM software in complex terrain and geological modeling to ensure a high degree of realism between the 3D geological model and the actual situation. It enables data conversion between the BIM geological model and the seepage numerical model, overcomes the shortcomings of mainstream international seepage calculation software in preprocessing work such as complex terrain and lithology in mountainous areas, improves the efficiency and accuracy of groundwater numerical model construction in mountainous areas, and provides a foundation for the evolution analysis of groundwater seepage fields in mountainous areas. Summary of the Invention

[0004] In view of the above-mentioned problems, the present invention is proposed.

[0005] Therefore, the problem to be solved by this invention is: how to accurately convert geological models and seepage numerical models.

[0006] To address the aforementioned technical problems, this invention provides the following technical solution: a geological model to seepage model conversion method, comprising: establishing a three-dimensional topographic model and a three-dimensional geological model of the study area; discretizing the three-dimensional topographic model and the three-dimensional geological model of the study area through data conversion, discretizing them into layers according to different strata to form ASCII files; importing the converted uniformly discretized node coordinate data into seepage numerical calculation layer by layer to construct a hydrogeological model; and constructing a seepage numerical calculation model by assigning values ​​to hydrogeological parameters and source-sink terms.

[0007] As a preferred embodiment of the geological model seepage model method described in this invention, the three-dimensional terrain model of the study area includes: establishing terrain surfaces based on contour lines and elevation points according to the topographic map of the study area; establishing terrain volumes based on the terrain surfaces and the modeling range; and establishing a three-dimensional terrain model of the study area. The three-dimensional geological model of the study area includes: calculating a preliminary interpolation result for each location (x, y, z) in the entire study area using the DSI interpolation algorithm; adjusting the interpolation result v(x, y, z) for each sample point according to fuzzy constraints; and applying precise constraints to certain specific sample points to form a three-dimensional geological model. The DSI interpolation algorithm is expressed as follows:

[0008] P = {p1, p2, p3, ..., p} n}

[0009]

[0010] Where P is the set of sample points in the borehole lithological stratification, n is the number of sample points, v(x, y, z) is the interpolation function, and v i Let p be the geological attribute value of the i-th sample point. i Let d(p) be the i-th sample point, α be the interpolation attenuation rate, and d(p) be the interpolation attenuation rate. i Let (x, y, z) be the point p. i The Euclidean distance to position (x, y, z); the fuzzy constraint is expressed as:

[0011] v i -δ i ≤v(x,y,z)≤v i +δ i

[0012]

[0013] Where, δ i Let β be the fuzzy range of the i-th sample point, and v be the degree of fuzziness. j Let j be the geological attribute value of the j-th sample point; the precise constraint is expressed as:

[0014] v(x j y j , z j ) = v j

[0015] Where, x j y j , z j The three-dimensional coordinates of the sample points for precise constraints.

[0016] As a preferred embodiment of the geological model transseepage model method described in this invention, the discretization is represented as follows:

[0017]

[0018]

[0019]

[0020]

[0021] Where D is the fractal dimension, ε is the mesh size, N(ε) is the number of meshes in the study area containing the model, and X u Y v Z k Let x be the coordinates of a discrete point. min y min z min Let x be the minimum coordinate of the study area. max y max z max Let u, v, and k be the minimum coordinates of the study area, where u, v, and k are integer indices.

[0022] As a preferred embodiment of the geological model transseepage model method described in this invention, the ASCII file is represented as follows:

[0023] f(X u Y v Z k ) = sin(X u )·cos(Y v )+tan(Z k )+ln(1+X u 2 +Y v 2 +Z k 2 )

[0024] Line u,v,k =str(f(X) u Y v Z k ))

[0025] Where, f(X) u Y v Z k ) is the encoding function, Line u,v,k This is a line in an ASCII file, representing the encoded value of a discrete point. str is the number converted to a string.

[0026] As a preferred embodiment of the geological model transformation and seepage modeling method of the present invention, the hydrogeological model includes defining the three-dimensional structure of the model using coded data in a generated ASCII file, taking advantage of the characteristics of modeling commands.

[0027] As a preferred embodiment of the geological model-to-seepage flow model method of the present invention, the hydrogeological parameters assigned include permeability, pore size, and water compressibility coefficient; the permeability includes, if f(X) u Y v Z k The penetration rate threshold T was not exceeded. s Then the baseline permeability value of the study area, if f(X) u Y v Z k Exceeding the penetration threshold T s The penetration rate is then expressed as:

[0028]

[0029] Where s is the penetration rate, s d T is the benchmark value for permeability. s σ is the penetration threshold. s 2 To control the width coefficient of the exponential function; if the current water level H exceeds the predetermined value of the groundwater level, the permeability threshold T is automatically adjusted. s , represented as:

[0030] T s ′=T s +γ(HH y )

[0031] Among them, T s ′ represents the adjusted permeability threshold, and H represents the current water level. y The groundwater level is a predetermined value, and γ is an adjustment coefficient; the pore size includes, if f(X) u Y v Z k Exceeding the pore gap threshold T n Then the reference value of the hole gap in the study area is used. If f(X) u Y v Z k The gap threshold T was not exceeded. n The hole gap is then expressed as:

[0032] n = n d +sin(T n -f(X u Y v Z k ))

[0033] Where n is the pore clearance, n d T is the reference value for hole clearance. n The pore space threshold is set; if the current soil saturation S is lower than the predetermined soil saturation value, the pore space threshold T is automatically adjusted. n , represented as:

[0034] T n ′=T n -δ(S y -S)

[0035] Among them, T n ′ represents the adjusted pore size threshold, and S represents the current soil saturation. y The soil saturation is a predetermined value, and δ is an adjustment coefficient; the compressibility coefficient of water includes, if f(X) u Y v Z k (Below the lower limit of the compressibility coefficient T) β1 The compression coefficient is expressed as:

[0036] β=β d -λ1×(T β1 -f(X u Y v Z k ))

[0037] If f(X) u Y v Z k (Higher than the upper limit of the compressibility coefficient T) β2 The compression coefficient is expressed as:

[0038] β=β d +λ2×(f(X u Y v Z k )-T β2 )

[0039] If f(X) u Y v Z k (Higher than the lower limit of compressibility T) β1 Below the upper limit of the compressibility coefficient T β2 The compression coefficient is expressed as:

[0040] β=β d +μ·cos(f(X u Y v Z k ))

[0041] Where β is the compressibility coefficient of water, β dHere, λ1 is the lower limit adjustment factor, λ2 is the upper limit adjustment factor, μ is the adjustment factor within the limit, and T is the reference value for the compressibility coefficient of water. β1 T is the lower limit of the compressibility factor. β2 This represents the upper limit of the compression factor.

[0042] As a preferred embodiment of the geological model transformation flow model method described in this invention, wherein: the source-sink term assignment includes, if f(X) u Y v Z k Exceeding the source-sink threshold T Q Then use the source-sink term baseline value of the study area, if f(X) u Y v Z k The threshold T for source and sink items has not been exceeded. Q Then the source and sink terms are represented as:

[0043] Q = Q d +log(1+|f(X u Y v Z k )-T Q |)

[0044] Where Q represents the source and sink terms, Q d T is the baseline value for source and sink terms. Q The source-sink threshold is set; if the current groundwater flow velocity V exceeds the predetermined value, the source-sink threshold T is automatically adjusted. Q , represented as:

[0045] T Q ′=T Q +τ(VV y )

[0046] Among them, T Q ′ represents the adjusted source-sink threshold, and V represents the current groundwater flow velocity. y The groundwater flow velocity is a predetermined value, and τ is an adjustment coefficient; the seepage numerical calculation model includes the construction of the seepage numerical calculation model based on the assignment of values ​​to hydrogeological parameters and source and sink terms.

[0047] Another objective of this invention is to provide a system for a geological model transseepage model method, which solves the problem of a geological model transseepage model by constructing a geological model transseepage model system.

[0048] To address the aforementioned technical problems, this invention provides the following technical solution: a geological model-to-seepage flow model system, comprising a data collection module, a three-dimensional model construction module, a data discretization module, and a hydrogeological model construction module; the data collection module collects topographic maps, contour lines, and elevation point information of the study area, as well as borehole lithological stratification, sample points, and geological attribute values; the three-dimensional model construction module establishes topographic surfaces using contour lines and elevation points based on the topographic data, and uses DSI interpolation algorithm and fuzzy constraints to form a three-dimensional geological model; the data discretization module discretizes the three-dimensional topography and geological model into data layered according to different strata, and outputs the discretized data as an ASCII file; the hydrogeological model construction module uses the encoded data in the ASCII file to define the three-dimensional structure of the model and assigns values ​​to hydrogeological parameters and source / sink terms.

[0049] A computer device includes a memory and a processor, the memory storing a computer program, characterized in that the processor executes the computer program to implement the steps of the geological model seepage flow model method as described above.

[0050] A computer-readable storage medium having a computer program stored thereon, characterized in that, when the computer program is executed by a processor, it implements the steps of the geological model-to-seepage flow model method as described above.

[0051] The beneficial effects of this invention are as follows: The geological model to seepage model method provided by this invention utilizes the advantages of BIM software in complex terrain and geological modeling to ensure that the three-dimensional geological model is highly consistent with the actual situation; through data conversion, this technology realizes the data conversion between the BIM geological model and the groundwater numerical model, overcoming the shortcomings of mainstream international software such as FEFLOW in pre-processing work such as complex terrain and lithology in mountainous areas; reducing a large amount of pre-processing workload and improving the efficiency and accuracy of groundwater numerical model construction in mountainous areas. Attached Figure Description

[0052] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. Wherein:

[0053] Figure 1 The first embodiment of the present invention provides an overall flowchart of a geological model transformation flow modeling method.

[0054] Figure 2 This is a structural diagram of a geological model transseepage model system provided in the second embodiment of the present invention.

[0055] Figure 3 This is a three-dimensional topographic model of the study area provided by a geological model transformation flow model method in the third embodiment of the present invention.

[0056] Figure 4 This is a three-dimensional BIM geological model map of the study area provided by a geological model transformation and seepage modeling method in the third embodiment of the present invention.

[0057] Figure 5 The study area grid model diagram of a geological model transformation flow model method provided in the third embodiment of the present invention.

[0058] Figure 6 The third embodiment of the present invention provides a groundwater seepage model diagram constructed using FEFLOW software for a geological model to seepage model method. Detailed Implementation

[0059] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0060] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0061] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in different places in this specification does not necessarily refer to the same embodiment, nor is it a single or selective embodiment that is mutually exclusive with other embodiments.

[0062] Example 1

[0063] Reference Figure 1 This is the first embodiment of the present invention, which provides a geological model to seepage model method, including: establishing a three-dimensional topographic model and a three-dimensional geological model of the study area; discretizing the three-dimensional topographic model and the three-dimensional geological model of the study area through data conversion, discretizing them into layers according to different strata, forming ASCII files; importing the converted uniformly discretized node coordinate data into seepage numerical calculation layer by layer to construct a hydrogeological model; and constructing a seepage numerical calculation model by assigning values ​​to hydrogeological parameters and source and sink terms.

[0064] Step 1: Construction of 3D terrain and geological models.

[0065] Based on the topographic map of the study area, topographic surfaces are established according to contour lines and elevation points. A "topographic volume" is then created based on the topographic surfaces and the modeling area, thus establishing a three-dimensional topographic model of the study area. Next, based on the lithological stratification of boreholes, the DSI interpolation algorithm is used for modeling. Through fuzzy and precise constraint techniques, it is ensured that the model accurately passes through known points and that the model's trend conforms to geological laws, forming a three-dimensional geological body model.

[0066] The DSI interpolation algorithm is used to calculate a preliminary interpolation result for each location (x, y, z) in the entire study area. For each sample point, the interpolation result v(x, y, z) is adjusted according to fuzzy constraints. For certain specific sample points, precise constraints are applied to form a three-dimensional geological model. The DSI interpolation algorithm is expressed as follows:

[0067] P = {p1, p2, p3, ..., p} n}

[0068]

[0069] Where P is the set of sample points in the borehole lithological stratification, n is the number of sample points, v(x, y, z) is the interpolation function, and v i Let p be the geological attribute value of the i-th sample point. i Let d(p) be the i-th sample point, α be the interpolation attenuation rate, and d(p) be the interpolation attenuation rate. i Let (x, y, z) be the point p. i The Euclidean distance to the position (x, y, z); the fuzzy constraint is expressed as:

[0070] v i -δ i ≤v(x,y,z)≤v i +δ i

[0071]

[0072] Where, δ i Let v be the fuzzy range of the i-th sample point, β be the fuzziness level, a constant between 0 and 1, and v j Let j be the geological attribute value of the j-th sample point; the precise constraint is expressed as:

[0073] v(x j y j , z j ) = v j

[0074] Where, x j y j , z j The three-dimensional coordinates of the sample points for precise constraints.

[0075] Step 2: Conversion between the three-dimensional geological model and the seepage calculation model.

[0076] After establishing a 3D BIM terrain model and geological model for the study area, the terrain and geological models were discretized into .dat files layered according to different strata through data conversion. The X, Y, and Z values ​​for each layer represent the coordinates X, Y, and elevation Z, respectively, forming readable ASCII files, thus enabling the conversion between the 3D BIM model and international seepage numerical simulation software. The discretization is represented as:

[0077]

[0078]

[0079]

[0080]

[0081] Where D is the fractal dimension, ε is the mesh size, N(ε) is the number of meshes in the study area containing the model, and X u Y v Z k Let x be the coordinates of a discrete point. min y min z min Let x be the minimum coordinate of the study area. max y max z max Here, u, v, and k represent the minimum coordinates of the study area, and are integer indices used to generate discrete points. For example, if 10 discrete points are desired in the X direction, the value of u can range from 0 to 9. The resulting readable ASCII file representation is as follows:

[0082] f(X u Y v Z k ) = sin(X u )·cos(Y v )+tan(Z k )+ln(1+X u 2 +Y v 2 +Z k 2 )

[0083] Line u,v,k =str(f(X) u Y v Z k ))

[0084] Where, f(X) u Yv Z k () is the encoding function that maps three-dimensional coordinates to a real number, Line. u,v,k This is a line in an ASCII file, representing the encoded value of a discrete point. str is the number converted to a string.

[0085] The transformed, uniformly discretized node coordinate data is imported layer by layer into the seepage numerical calculation software. The hydrogeological model is constructed by utilizing the features of the modeling command. Then, by assigning values ​​to the hydrogeological parameters and source-sink terms, the seepage numerical calculation model is completed.

[0086] Leveraging the features of the modeling commands, the 3D structure of the model is defined using encoded data from the generated ASCII file. Hydrogeological parameters are assigned, including permeability, pore size, and the compressibility of water.

[0087] If f(X) u Y v Z k The penetration rate threshold T was not exceeded. s Then the baseline permeability value of the study area, if f(X) u Y v Z k Exceeding the penetration threshold T s The penetration rate is then expressed as:

[0088]

[0089] Where s is the penetration rate, s d T is the benchmark value for permeability. s The penetration threshold, To control the width coefficient of the exponential function; if the current water level H exceeds the predetermined value of the groundwater level, the permeability threshold T is automatically adjusted. s , represented as:

[0090] T s ′=T s +γ(HH y )

[0091] Among them, T s ′ represents the adjusted permeability threshold, and H represents the current water level. y γ is the predetermined groundwater level, and γ is the adjustment coefficient.

[0092] If f(X) u Y v Z k Exceeding the pore gap threshold T n Then the reference value of the hole gap in the study area is used. If f(X) u Y v Z kThe gap threshold T was not exceeded. n The hole gap is then expressed as:

[0093] n = n d +sin(T n -f(X u Y v Z k ))

[0094] Where n is the pore clearance, n d T is the reference value for hole clearance. n The pore space threshold is set; if the current soil saturation S is lower than the predetermined soil saturation value, the pore space threshold T is automatically adjusted. n , represented as:

[0095] T n ′=T n -δ(S y -S)

[0096] Among them, T n ′ represents the adjusted pore size threshold, and S represents the current soil saturation. y δ is the predetermined value for soil saturation, and δ is the adjustment coefficient.

[0097] If f(X) u Y v Z k (Below the lower limit of the compressibility coefficient T) β1 The compression coefficient is expressed as:

[0098] β=β d -λ1×(T β1 -f(X u Y v Z k ))

[0099] If f(X) u Y v Z k (Higher than the upper limit of the compressibility coefficient T) β2 The compression coefficient is expressed as:

[0100] β=β d +λ2×(f(X u Y v Z k )-T β2 )

[0101] If f(X) u Y v Z k (Higher than the lower limit of compressibility T) β1 Below the upper limit of the compressibility coefficient T β2The compression coefficient is expressed as:

[0102] β=β d +μ·cos(f(X u Y v Z k ))

[0103] Where β is the compressibility coefficient of water, β d Here, λ1 is the lower limit adjustment factor, λ2 is the upper limit adjustment factor, μ is the adjustment factor within the limit, and T is the reference value for the compressibility coefficient of water. β1 T is the lower limit of the compressibility factor. β2 This represents the upper limit of the compression factor.

[0104] If f(X) u Y v Z k Exceeding the source-sink threshold T Q Then use the source-sink term baseline value of the study area, if f(X) u Y v Z k The threshold T for source and sink items has not been exceeded. Q Then the source and sink terms are represented as:

[0105] Q = Q d +log(1+|f(X u Y v Z k )-T Q |)

[0106] Where Q represents the source and sink terms, Q d T is the baseline value for source and sink terms. Q The source-sink threshold is set; if the current groundwater flow velocity V exceeds the predetermined value, the source-sink threshold T is automatically adjusted. Q , represented as:

[0107] T Q ′=T Q +τ(VV y )

[0108] Among them, T Q ′ represents the adjusted source-sink threshold, and V represents the current groundwater flow velocity. y The groundwater flow velocity is a predetermined value, and τ is an adjustment coefficient; the seepage numerical calculation model includes the construction of the seepage numerical calculation model based on the assignment of values ​​to hydrogeological parameters and source and sink terms.

[0109] Example 2

[0110] Reference Figure 2This is the second embodiment of the present invention, which differs from the previous embodiment in that it provides a geological model seepage flow model system, including: a data collection module, a three-dimensional model construction module, a data discretization module, and a hydrogeological model construction module.

[0111] The data acquisition module collects topographic maps, contour lines, and elevation points of the study area, as well as borehole lithological stratification, sample points, and geological attribute values.

[0112] The 3D model building module uses contour lines and elevation points to create terrain surfaces based on terrain data, and uses DSI interpolation algorithm and fuzzy constraints to form a 3D geological model.

[0113] The data discretization module discretizes the 3D terrain and geological model into data layered according to different strata, and outputs the discretized data as an ASCII file.

[0114] The hydrogeological model building module uses coded data in ASCII files to define the three-dimensional structure of the model and assign values ​​to hydrogeological parameters and source / sink terms.

[0115] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0116] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device.

[0117] More specific examples of computer-readable media (a non-exhaustive list) include: electrical connections (electronic devices) having one or more wires, portable computer disk drives (magnetic devices), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Furthermore, computer-readable media can even be paper or other suitable media on which the program can be printed, because the program can be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in computer memory.

[0118] It should be understood that various parts of the present invention can be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0119] Example 3

[0120] Reference Figures 3-6 This is the third embodiment of the present invention, which differs from the previous two embodiments in that it is used to verify and explain the technical effects adopted in the present invention, so as to verify the real effect of the method.

[0121] Step 1: Establish a 3D BIM terrain model and geological model of the study area, such as... Figure 3 and Figure 4 As shown.

[0122] Step 2: Discretize the terrain model and geological model, and save them as *.txt files.

[0123] Step 3: Import the model boundary extent .dxf file into Surfer software. Select "data" in the Grid panel, input the X, Y, and Z data, and then perform interpolation and encryption processing to form a *.grid file. This establishes a Grid model of the study area, such as... Figure 5 As shown.

[0124] Step 4: Save the Grid model as a *.dat file, and then convert the data to form a readable ASCII file.

[0125] Step 5: The saved ASCII code file, where I, J, and K (I represents the X direction, J represents the y direction, and K represents the Z direction) uniquely determines the spatial location of the node. I represents the X direction, J represents the y direction, and K represents the Z direction. x, y, and z represent the actual coordinates and elevation of the node in the three directions, respectively.

[0126] Step Six: Based on the FEFLOW modeling concept and its modeling command characteristics, convert the saved uniformly discretized node coordinate data (.dat file) into a file that FEFLOW software can recognize. Then, establish a groundwater numerical simulation model in the FEFLOW software, such as... Figure 6 As shown.

[0127] This embodiment utilizes both the conventional method and our invented method for conversion simultaneously, and the detection comparison results are shown in the table below:

[0128] Table 1 Comparison between traditional methods and our invented methods

[0129] Judgment category Conventional method Our invention method Model construction time 48h 24h Model accuracy error 12％ 5％ Data conversion workload 40h 10h Topography and geology model fitting degree 85％ 95％

[0130] The comparison results above show that the model building time of our invention is 24 hours, which is 24 hours less than the 48 hours of the traditional method; the model accuracy error is 5%, which is 7% less than the 12% of the traditional method; the data conversion workload is 10 hours, which is 30 hours less than the 40 hours of the traditional method; and the terrain and geological model fitting degree is 95%, which is 10% more than the 85% of the traditional method.

[0131] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A geological model transformation and seepage model method, characterized in that: include, Establish a three-dimensional topographic model and a three-dimensional geological model of the study area; The three-dimensional terrain model and the three-dimensional geological model of the study area were discretized through data conversion, and then discretized into ASCII files according to different strata. The transformed uniformly discretized node coordinate data is imported into the seepage numerical calculation layer by layer to construct a hydrogeological model. By assigning values ​​to hydrogeological parameters and source-sink terms, a seepage numerical calculation model is constructed.

2. The geological model transformation and seepage model method as described in claim 1, characterized in that: The three-dimensional terrain model of the study area includes: establishing terrain surfaces based on contour lines and elevation points according to the topographic map of the study area; establishing terrain volumes based on terrain surfaces and modeling range; and establishing a three-dimensional terrain model of the study area. The three-dimensional geological model of the study area includes: using the DSI interpolation algorithm to calculate a preliminary interpolation result for each location (x, y, z) in the entire study area; adjusting the interpolation result v(x, y, z) according to fuzzy constraints for each sample point; and applying precise constraints to certain specific sample points to form a three-dimensional geological model. The DSI interpolation algorithm is expressed as follows: P={p1,p2,p3,...,p n } Where P is the set of sample points in the borehole lithological stratification, n is the number of sample points, and v(x, y, z) is the interpolation function. i Let p be the geological attribute value of the i-th sample point. i Let d(p) be the i-th sample point, α be the interpolation attenuation rate, and d(p) be the interpolation attenuation rate. i Let (x, y, z) be the point p. i Euclidean distance to position (x, y, z); The fuzzy constraint is expressed as follows: v i -δ i ≤v(x,y,z)≤v i +δ i Where, δ i Let β be the fuzzy range of the i-th sample point, and v be the degree of fuzziness. j Let j be the geological attribute value of the j-th sample point; The precise constraint is expressed as follows: v(x j ,y j ,z j )=v j Where, x j y j , z j The three-dimensional coordinates of the sample points for precise constraints.

3. The geological model transformation and seepage model method as described in claim 2, characterized in that: The discretization is expressed as follows: Where D is the fractal dimension, ε is the mesh size, N(ε) is the number of meshes in the study area containing the model, and X u Y v Z k Let x be the coordinates of a discrete point. min y min z min Let x be the minimum coordinate of the study area. max y max z max Let u, v, and k be the minimum coordinates of the study area, where u, v, and k are integer indices.

4. The geological model transformation flow model method as described in claim 3, characterized in that: The ASCII file is represented as follows: Line u,v,k =str(f(X u ,Y v ,Z k )) Where, f(X) u Y v Z k ) is the encoding function, Line u,v,k This is a line in an ASCII file, representing the encoded value of a discrete point. str is the number converted to a string.

5. The geological model to seepage flow model method as described in claim 4, characterized in that: The hydrogeological model includes defining the three-dimensional structure of the model using encoded data in a generated ASCII file, taking advantage of the features of the modeling command.

6. The geological model transformation and seepage model method as described in claim 5, characterized in that: The hydrogeological parameters assigned include permeability, pore size, and water compressibility. The permeability includes, if f(X) u Y v Z k The penetration rate threshold T was not exceeded. s Then the baseline permeability value of the study area, if f(X) u Y v Z k Exceeding the penetration threshold T s The penetration rate is then expressed as, Where s is the penetration rate, s d T is the benchmark value for permeability. s The penetration threshold, To control the width coefficient of the exponential function; If the current water level H exceeds the predetermined value of the groundwater level, the permeability threshold T will be automatically adjusted. s , is represented as , T s ′=T s +γ(H-H y ) Among them, T s ′ represents the adjusted permeability threshold, and H represents the current water level. y The groundwater level is a predetermined value, and γ is an adjustment coefficient. The aperture gap includes, if f(X) u Y v Z k Exceeding the pore gap threshold T n Then the reference value of the hole gap in the study area is used. If f(X) u Y v Z k The gap threshold T was not exceeded. n The hole gap is then expressed as, n=n d +sin(T n -f(X u ,Y v ,Z k )) Where n is the pore clearance, n d T is the reference value for hole clearance. n The aperture gap threshold; If the current soil saturation S is lower than the predetermined value for soil saturation, the pore size threshold T will be automatically adjusted. n , is represented as , T n ′=T n -δ(S y -S) Among them, T n ′ represents the adjusted pore size threshold, and S represents the current soil saturation. y δ is the predetermined value for soil saturation, and δ is the adjustment coefficient. The compressibility coefficient of water includes, if f(X) u Y v Z k (Below the lower limit of the compressibility coefficient T) β1 Then the compression coefficient is expressed as, β=β d -λ1×(T β1 -f(X u ,Y v ,Z k )) If f(X) u Y v Z k (Higher than the upper limit of the compressibility coefficient T) β2 Then the compression coefficient is expressed as, β=β d +λ2×(f(X u ,Y v ,Z k )-T β2 ) If f(X) u Y v Z k (Higher than the lower limit of compressibility T) β1 Below the upper limit of the compressibility coefficient T β2 Then the compression coefficient is expressed as, β=β d +μ·cos(f(X u ,Y v ,Z k )) Where β is the compressibility coefficient of water, β d Here, λ1 is the lower limit adjustment factor, λ2 is the upper limit adjustment factor, μ is the adjustment factor within the limit, and T is the reference value for the compressibility coefficient of water. β1 T is the lower limit of the compressibility factor. β2 This represents the upper limit of the compression factor.

7. The geological model transformation and seepage model method as described in claim 6, characterized in that: The source and sink term assignments include, if f(X) u Y v Z k Exceeding the source-sink threshold T Q Then use the source-sink term baseline value of the study area, if f(X) u Y v Z k The threshold T for source and sink items has not been exceeded. Q Then the source and sink terms are represented as follows: Q=Q d +log(1+|f(X u ,Y v ,Z k )-T Q |) Where Q represents the source and sink terms, Q d T is the baseline value for source and sink terms. Q The threshold for source and sink terms; If the current groundwater flow velocity V exceeds the predetermined value, the source and sink threshold T will be automatically adjusted. Q , is represented as , T Q ′=T Q +τ(V-V y ) Among them, T Q ′ represents the adjusted source-sink threshold, and V represents the current groundwater flow velocity. y τ is a predetermined value for groundwater flow velocity, and τ is an adjustment coefficient. The seepage numerical calculation model includes the construction of the seepage numerical calculation model based on the assignment of values ​​to hydrogeological parameters and source-sink terms.

8. A system employing a geological model to seepage model method as described in any one of claims 1 to 7, characterized in that: It includes a data collection module, a 3D model building module, a data discretization module, and a hydrogeological model building module; The data acquisition module collects topographic maps, contour lines, and elevation points of the study area, as well as borehole lithological stratification, sample points, and geological attribute values. The 3D model building module establishes terrain surfaces using contour lines and elevation points based on terrain data, and uses DSI interpolation algorithm and fuzzy constraints to form a 3D geological model. The data discretization module discretizes the three-dimensional terrain and geological model into data layered according to different strata, and outputs the discretized data as an ASCII file; The hydrogeological model construction module uses coded data in ASCII files to define the three-dimensional structure of the model and assigns values ​​to hydrogeological parameters and source-sink terms.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that: When the processor executes the computer program, it implements the steps of the geological model seepage model method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by the processor, it implements the steps of the geological model seepage model method according to any one of claims 1 to 7.

Images