Method, device and equipment for calculating recoverable resource quantity of salt lake brine and storage medium

Abstract

The invention provides a salt lake brine recoverable resource quantity calculation method, device and equipment and a storage medium, and relates to the technical field of salt lake resource evaluation and development, the method comprises the steps that the storage resource quantity of salt lake brine in a target area is calculated, and the storage resource quantity comprises the specific yield resource quantity and the elastic resource quantity; calculating the supply resource quantity of the salt lake brine in the target area; and adding the specific yield resource quantity, the elastic resource quantity and the replenishment resource quantity in an arithmetic manner to obtain the recoverable resource quantity of the salt lake brine in the target area. The method has the technical effect that the recoverable resource quantity of the salt lake brine is accurately evaluated.

Description

Technical Field

[0001] The present application relates to the technical field of salt lake resource evaluation and development, and specifically to a method, device, equipment and storage medium for calculating the recoverable amount of salt lake brine resources. Background Art

[0002] Salt lake brine is an important mineral resource, rich in potassium, lithium, boron, magnesium and other elements. Accurately assessing the recoverable resources of salt lake brine is the premise and basis for the rational development and utilization of this resource.

[0003] At present, a common method is to directly use the drilling fluid to calculate the static reserves of salt lake brine by area extrapolation or volume estimation. This method is simple to operate, but due to the dynamic replenishment process of salt lake brine, the above method often underestimates the recoverable potential of salt lake brine resources, resulting in inaccurate assessment results of recoverable resources of salt lake brine. Summary of the invention

[0004] The present application provides a method, device, equipment and storage medium for calculating the recoverable resources of salt lake brine, which are used to accurately evaluate the recoverable resources of salt lake brine.

[0005] In the first aspect, the present application provides a method for calculating the recoverable resources of salt lake brine, the method comprising: calculating the storage resources of salt lake brine in a target area, the storage resources comprising water supply resources and elastic resources; calculating the replenishment resources of salt lake brine in the target area; arithmetically adding the water supply resources, the elastic resources and the replenishment resources to obtain the recoverable resources of salt lake brine in the target area.

[0006] By adopting the above technical solution, the storage resources and replenishment resources of the salt lake brine in the target area are calculated, and the recoverable resources of the salt lake brine are obtained by adding the two together, thereby realizing a comprehensive assessment of the salt lake brine resources. Among them, the storage resources take into account the static water resources, i.e., the water supply resources, and the dynamic water resources, i.e., the elastic resources. Through the calculation of static water resources, dynamic water resources, and replenishment resources, the total resources of the salt lake brine body can be more comprehensively reflected, so that the calculated value of the recoverable resources of the salt lake brine can be obtained, providing a basis for the development and utilization of the salt lake brine. This calculation method takes into account both static water and dynamic water resources, as well as the replenishment effect, and accurately evaluates the recoverable resources of the salt lake brine.

[0007] Optionally, the calculation of the storage resource volume of salt lake brine in the target area includes: obtaining the area of ​​salt lake brine in the target area, the porosity of the rock of the salt lake brine reservoir, the thickness of the salt lake brine layer, the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine, and the top elevation of the pressure brine reservoir of the salt lake brine; calculating the water supply resource volume according to the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir, and the thickness of the brine layer; calculating the elastic resource volume according to the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine, and the top elevation of the pressure brine reservoir of the salt lake brine; and arithmetically adding the water supply resource volume and the elastic resource volume to obtain the storage resource volume of salt lake brine in the target area.

[0008] By adopting the above technical solution, by obtaining the relevant parameters of salt lake brine, including area, reservoir parameters, head parameters, etc., according to the water supply resource calculation formula and elastic resource calculation formula, the water supply resource and elastic resource are calculated respectively, and the two are added to obtain the storage resource of salt lake brine. Among them, the water supply resource calculation takes into account the area, reservoir rock porosity and layer thickness of the salt lake brine, reflecting the static water resource of the salt lake brine body; the elastic resource calculation takes into account the elastic water release characteristics and head parameters of the salt lake brine, reflecting the dynamic water resource of the salt lake brine body. This calculation method can more comprehensively reflect the total storage resource of the salt lake brine body by reasonably obtaining parameters, using formulas to calculate the water supply and elastic resources, and adding the two together, and can more accurately evaluate the total reserves of salt lake brine, providing important parameter basis for subsequent development and utilization.

[0009] Optionally, the calculating of the water yield resource quantity according to the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir and the thickness of the brine layer comprises: substituting the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir and the thickness of the brine layer into a first formula to calculate the water yield resource reserves; wherein the first formula is: In the formula, Q1 is the water resource reserve, A is the area of ​​the salt lake brine, is the porosity of the rock of the salt lake brine reservoir, and M is the thickness of the salt lake brine layer.

[0010] By adopting the above technical solution, a water supply resource calculation formula is introduced when calculating the water supply resource. When calculating, the relevant parameters of the salt lake brine are substituted into the formula for calculation, and the water supply resource reserves can be directly obtained. This formula fully considers the key factors that determine the water supply of salt lake brine, including three parameters: salt lake area, reservoir porosity and layer thickness. Using this formula can avoid tedious calculation processes and directly solve the water supply resource based on key parameters, thereby achieving the purpose of simplifying and quickly calculating the water supply resource of salt lake brine. It provides convenience for subsequent salt lake brine resource assessment.

[0011] Optionally, the elastic resource volume is calculated according to the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine, including: substituting the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine into a second formula to calculate the elastic resource reserves; wherein the second formula is: Q2=u * (hH)A; wherein Q2 is the elastic resource volume, A is the elastic water release coefficient of the salt lake brine, h is the average pressure head elevation of the salt lake brine, H is the top elevation of the pressure brine reservoir of the salt lake brine, and A is the area of ​​the salt lake brine.

[0012] By adopting the above technical solution, the elastic resource calculation formula is introduced when calculating the elastic resource volume. When calculating, the relevant parameters of the salt lake brine are substituted into the formula for calculation, and the elastic resource reserves can be directly obtained. This formula fully considers the key factors that determine the elastic water release of salt lake brine, including the three parameters of elastic water release coefficient, head difference and area. The use of this formula can avoid the complicated calculation process and directly solve the elastic resource volume based on the key parameters, thereby achieving the purpose of simplifying and quickly calculating the elastic resource volume of salt lake brine. It provides convenience for the subsequent evaluation of salt lake brine resources.

[0013] Optionally, the calculation of the replenishment resource amount of salt lake brine in the target area includes: establishing a hydrogeological conceptual model of salt lake brine in the target area; based on the hydrogeological conceptual model, using visualization software as a water flow simulation platform for salt lake brine in the target area, dividing the target area into a regular rectangular grid of several rows and columns on a plane, wherein the regular rectangular grid includes a plurality of rectangular units; using the row and column numbers of the regular rectangular grid as the numbers of the rectangular units, and establishing a three-dimensional finite difference model; calculating the underground horizontal measurement of each of the rectangular units through the three-dimensional finite difference model, and calculating the replenishment of salt lake brine in the target area based on the underground horizontal measurement of each of the rectangular units, and using the replenishment as the replenishment resource amount of salt lake brine in the target area.

[0014] By adopting the above technical scheme, when calculating the replenishment resources, the hydrogeological conceptual model of salt lake brine is first established, and then the visualization software is used as a water flow simulation platform to divide the target area into regular rectangular grids and establish a three-dimensional finite difference model. The underground horizontal measurement of each unit is calculated in the model, and the replenishment of salt lake brine is calculated as the replenishment resources. This technical scheme fully considers the hydrogeological process of salt lake brine replenishment, uses numerical simulation methods, simulates the replenishment mechanism of salt lake brine based on the hydrogeological model, and realizes the calculation of replenishment resources. At the same time, the refinement of replenishment simulation is realized by using visualization and three-dimensional finite difference numerical calculation methods. This technical scheme can effectively simulate the complex salt lake recharge process, and the calculation results are more in line with the actual situation, providing important parameter support for the evaluation of salt lake brine resources.

[0015] Optionally, the establishment of a hydrogeological conceptual model of salt lake brine in the target area includes: Where Ω is the seepage area, ε is the source-sink coefficient of the aquifer, h is the groundwater level, h0 is the initial water level distribution of the groundwater system, and K x , K y are the horizontal permeability coefficients of the water-containing medium in the x and y directions, K z is the vertical z permeability coefficient of the water-bearing medium, S is the water storage rate of the aquifer below the free surface, μ is the gravity water supply of the phreatic aquifer, Γ0 is the upper boundary of the seepage area, that is, the free surface of the groundwater, p is the evaporation and precipitation infiltration intensity of the phreatic surface, h1 is the known boundary water level value, Γ1 is the known water level boundary, Γ2 is the flow boundary of the seepage area, is the normal direction of the boundary surface, K n is the permeability coefficient in the normal direction of the boundary surface, q is the flow rate per unit area of ​​the flow boundary, inflow is positive, outflow is negative, and it is 0 when it is isolated from water, and t is the research time.

[0016] By adopting the above technical scheme, when establishing the hydrogeological conceptual model, the groundwater flow control equation was used to describe the hydrogeological process of salt lake brine. The control equation fully considers the various hydrological and geological factors that affect the recharge of salt lake brine, including seepage area, aquifer parameters, source and sink terms, boundary conditions, etc. The use of this control equation to establish a hydrogeological model can accurately describe the recharge mechanism of salt lake brine and lay a theoretical foundation for subsequent numerical simulation calculations. The form of this control equation comprehensively and systematically considers various hydrological and geological factors that regulate recharge, and the established hydrogeological model can truly reflect the law of salt lake recharge. This technical means can effectively improve the scientificity and accuracy of the model, provide reliable model support for the subsequent calculation of recharge volume, and make the recharge resource assessment results more valuable for reference.

[0017] Optionally, the calculating the underground horizontal measurement of each of the rectangular units through the three-dimensional finite difference model includes: setting hydrogeological parameters in the three-dimensional finite difference model, the hydrogeological parameters including the spatial position, lithological parameters and boundary conditions of the aquifer; and numerically simulating the groundwater flow of the salt lake brine according to the spatial position, lithological parameters and boundary conditions of the aquifer to obtain the underground horizontal measurement of each of the rectangular units.

[0018] By adopting the above technical scheme, when calculating the underground level measurement, the hydrogeological parameters are first set in the three-dimensional finite difference model, including the location of the aquifer, lithology and boundary conditions. Then, according to the set hydrogeological parameters, the three-dimensional finite difference model is used to numerically simulate the groundwater flow process of the salt lake brine to obtain the underground level measurement of each rectangular unit. This technical scheme uses the three-dimensional finite difference method to numerically calculate the complex salt lake brine recharge process, which can avoid the empirical assumptions in the conceptual model and realize the objectification of the calculation of the recharge amount. At the same time, setting the hydrogeological parameters provides important constraints for the model calculation. This technical scheme realizes the effective combination of the model parameterization process and the numerical calculation process, so that the recharge calculation is more in line with the actual hydrogeological conditions of the salt lake brine. The calculation results can provide more reliable technical support for the subsequent salt lake brine resource assessment.

[0019] In the second aspect, the present application provides a device for calculating the recoverable resources of salt lake brine, the device comprising: a first calculation module, a second calculation module and a generation module; wherein the first calculation module is used to calculate the storage resources of salt lake brine in the target area, the storage resources including the water supply resources and the elastic resources; the second calculation module is used to calculate the replenishment resources of salt lake brine in the target area; the generation module is used to arithmetically add the water supply resources, the elastic resources and the replenishment resources to obtain the recoverable resources of salt lake brine in the target area.

[0020] In the third aspect, the present application provides an electronic device, which adopts the following technical solution: it includes a processor, a memory, a user interface and a network interface, the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes a computer program such as any of the above-mentioned methods for calculating the recoverable resources of salt lake brine.

[0021] In a fourth aspect, the present application provides a computer-readable storage medium, which adopts the following technical solution: a computer program that can be loaded by a processor and execute any of the above-mentioned methods for calculating the recoverable resources of salt lake brine is stored.

[0022] In summary, the present application includes at least one of the following beneficial technical effects: By calculating the storage resources and replenishment resources of salt lake brine in the target area, and adding the two together to obtain the recoverable resources of salt lake brine, a comprehensive assessment of the salt lake brine resources is achieved. Among them, the storage resources take into account the static water resources, i.e., the water supply resources, and the dynamic water resources, i.e., the elastic resources. Through the calculation of static water resources, dynamic water resources, and replenishment resources, the total resources of the salt lake brine body can be more comprehensively reflected, so that the calculated value of the recoverable resources of salt lake brine can be obtained, providing a basis for the development and utilization of salt lake brine. This calculation method takes into account both static water and dynamic water resources, as well as the replenishment effect, and accurately evaluates the recoverable resources of salt lake brine. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 It is a flow chart of a method for calculating the recoverable resources of salt lake brine provided in an embodiment of the present application; Figure 2 It is a schematic diagram of the structure of a device for calculating the recoverable resources of salt lake brine provided in an embodiment of the present application; Figure 3 It is a structural schematic diagram of an electronic device provided in an embodiment of the present application.

[0024] Description of reference numerals: 1000, electronic device; 1001, processor; 1002, communication bus; 1003, user interface; 1004, network interface; 1005, memory. DETAILED DESCRIPTION

[0025] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below in conjunction with the drawings in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments.

[0026] In the description of the embodiments of the present application, words such as "illustrative", "for example" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "illustrative", "for example" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "illustrative", "for example" or "for example" is intended to present related concepts in a concrete way.

[0027] Figure 1 Schematic diagram of a method for calculating the recoverable resources of salt lake brine provided in the embodiment of the present application. Figure 1 As shown, the method includes S101-S103.

[0028] S101, calculating the storage resource volume of salt lake brine in the target area, where the storage resource volume includes water supply resource volume and elastic resource volume.

[0029] In one example, the present invention provides a method for calculating the recoverable resources of salt lake brine. First, the storage resources of the salt lake brine in the target area are calculated, and the storage resources include two parts: water yield resources and elastic resources. The reason for calculating the storage resources is that it is an important component of the recoverable resources of the salt lake brine, and it characterizes the resource reserves of the brine deposit under static conditions. Among them, the water yield resources reflect the amount of brine that can be released from the brine reservoir under the action of gravity, while the elastic resources reflect the amount of brine that can be discharged from the brine reservoir under the action of elastic pressure. In order to calculate the water yield resources, it is necessary to obtain parameters such as the salt lake brine area (A), the porosity of the brine reservoir rock (φ), and the brine layer thickness (M) in the target area, and then substitute them into the first formula Calculate, where Q1 is the water resource reserve, A is the area of ​​salt lake brine, is the porosity of the salt lake brine reservoir rock, and M is the thickness of the salt lake brine layer.

[0030] The brine area, rock porosity and layer thickness jointly determine the water release capacity of the brine reservoir under gravity. Similarly, to calculate the elastic resource volume, it is necessary to obtain the elastic water release coefficient (u), the average pressure head elevation (h) and the top elevation (H) of the brine in the salt lake, and substitute them into the second formula Q2 = u * (hH)A is used for calculation, where Q2 is the elastic resource volume, A is the elastic water release coefficient of the salt lake brine, h is the average pressure head elevation of the salt lake brine, H is the top elevation of the pressure brine reservoir of the salt lake brine, and A is the area of ​​the salt lake brine.

[0031] The elastic water release coefficient reflects the compressibility of rock and brine, the difference between the average pressure head elevation and the brine reservoir roof elevation reflects the pressure on the pressure brine, and the area parameter reflects the spatial range of the pressure. The three together determine the water release capacity of the brine reservoir under elastic pressure. It should be pointed out that the above two formulas are derived on the basis of continuous medium mechanics and seepage mechanics, and have clear physical concepts and applicable conditions.

[0032] Through the above calculations, the water supply resource and elastic resource of the salt lake brine in the target area can be obtained, and the storage resource can be obtained by adding the two. This resource can more comprehensively reflect the resource occurrence of brine deposits under the dual effects of gravity and elastic pressure, and is an important basis for the subsequent calculation of recoverable resources. At the same time, by distinguishing between gravity resources and elastic resources, the spatial distribution characteristics of brine resources can also be revealed, enrichment areas can be identified, and a basis for the reasonable deployment of mining projects can be provided. In summary, the calculation of storage resources is a key link in the evaluation of recoverable resources of salt lake brine. The calculation idea of ​​gravity water supply plus elastic water release is more scientific and reasonable, which can fully consider the geological conditions and dynamic factors of brine occurrence. The calculation results are highly reliable and have important value in guiding the development and utilization of mining areas.

[0033] S102, calculating the replenishment resource amount of salt lake brine in the target area.

[0034] In one example, after calculating the storage resources, the replenishment resources of the salt lake brine in the target area are further calculated. The reason for calculating the replenishment resources is that salt lake brine, as a special underground fluid, has obvious fluidity and dynamic change characteristics. During the mining process, as the brine level drops, the brine in the surrounding area will be replenished to the mining center driven by the hydraulic gradient, forming a certain amount of lateral runoff replenishment. This part of the replenishment can be used as the newly recovered resources under dynamic mining conditions, but the traditional static reserve calculation method often ignores this phenomenon. In order to scientifically evaluate the recoverable potential of salt lake brine resources, it is necessary to quantitatively calculate and predict the brine replenishment resources.

[0035] The present invention adopts a numerical simulation method to achieve the purpose. First, according to the hydrogeological conditions of the mining area, a three-dimensional seepage conceptual model of salt lake brine is established. On the basis of the conceptual model, mature numerical simulation software such as GMS, Visual Modflow, etc. are used to divide the target area into regular rectangular grids on the plane, and the row and column numbers are used as the numbers of the grid units, thereby establishing a three-dimensional finite difference model. In the finite difference model, it is necessary to set corresponding hydrogeological parameters, including lithological parameters such as the spatial position, permeability, water supply degree, and storage coefficient of the aquifer, as well as external conditions such as the initial water level and boundary conditions. These parameters can be obtained by means of indoor tests, pumping tests, tracer tests, etc., or by referring to regional empirical values. Based on the three-dimensional finite difference model, the groundwater flow of brine in the target area is numerically simulated, and the underground horizontal measurement of each grid unit can be calculated by solving a group of nonlinear partial differential equations, mainly including the lateral runoff recharge. The recharge of each unit is accumulated to obtain the total recharge of brine in the target area, which is included in the recoverable resource as the recharge resource.

[0036] The advantage of the numerical simulation method is that it can accurately depict the seepage and migration process of brine in porous media, reasonably reflect the dynamic changes caused by mining, and provide a reliable technical means for quantitative evaluation of replenishment resources. By conducting sensitivity analysis and uncertainty analysis, the uncertainty of model parameters and their impact on simulation results can also be diagnosed, providing improvement directions for subsequent work.

[0037] In addition, based on the model verification, this method can also be used to predict the changes in brine recharge under future mining plans, so as to dynamically optimize the mining deployment and realize the rational development and utilization of brine resources in the mining area. In summary, the calculation of recharge resources is an indispensable and important part of the evaluation of recoverable resources of salt lake brine. The use of numerical simulation methods can fully consider the dynamic effects of brine mining, improve the scientificity and accuracy of resource estimation, and is of great significance for guiding the scientific management of mining areas.

[0038] Based on the above embodiment, as an optional implementation, in S102, calculating the replenishment resource amount of salt lake brine in the target area specifically includes: S201, establish a hydrogeological conceptual model of salt lake brine in the target area.

[0039] In one example, the hydrogeological conceptual model is an important basis for conducting salt lake brine resource evaluation. It is guided by hydrogeological theory and comprehensively considers multiple factors such as the origin of salt lakes, tectonic background, stratigraphic lithology, and groundwater dynamic conditions. It provides a generalized description of the spatial framework of salt lake brine storage, recharge runoff discharge processes, and hydraulic connections.

[0040] The specific approach to constructing a hydrogeological conceptual model is to first collect and organize regional geological and hydrogeological data, and through indoor and outdoor surveys, remote sensing interpretation and other means, identify the basic geological elements such as stratum lithology, geological structure, brine-bearing layer structure, and boundary conditions in the region; secondly, use borehole pumping tests, tracer tests, isotope dating and other methods to analyze and study the hydrogeological parameters such as the permeability, water supply, hydrochemical characteristics, and brine recharge runoff conditions of the brine-bearing layer; thirdly, explore the control of regional tectonic evolution, climate and environmental changes, magmatic hydrothermal activities and other geological effects on brine enrichment; finally, based on the principles of system analysis and pattern recognition, summarize the basic characteristics of salt lake brine storage and the key processes of migration and transformation from three levels: spatial structure, field problems, and boundary conditions, and express them intuitively with graphic symbols. Through the above process, a relatively complete hydrogeological conceptual model can be established.

[0041] On the one hand, this model can explain the controlling factors of brine enrichment from the perspective of genetic mechanism, identify brine-rich basins, find brine-rich strata, and reveal the temporal and spatial differentiation of brine mineralization. On the other hand, it can provide important geological basis and boundary constraints for subsequent three-dimensional geological modeling, numerical simulation, mining planning and other work, and is the logical starting point for quantitative evaluation of salt lake brine resources. It should be pointed out that due to the complex and changeable occurrence environment of salt lake brine, the principle of adapting measures to local conditions and highlighting key points should be followed when constructing conceptual models, and appropriate model types should be selected for different types of salt lakes. For closed salt lakes, the focus should be on brine vertical recharge and brine circulation at the basin scale; for semi-closed salt lakes, the focus should be on peripheral recharge and regional brine migration; and for open salt lakes, the focus should be on atmospheric precipitation, surface water and its transformation relationship with underground brine.

[0042] Establishing the hydrogeological conceptual model of salt lake brine in the target area specifically includes: Where Ω is the seepage area, ε is the source-sink coefficient of the aquifer, h is the groundwater level, h0 is the initial water level distribution of the groundwater system, and K x , K y are the horizontal permeability coefficients of the water-containing medium in the x and y directions, K z is the vertical z permeability coefficient of the water-bearing medium, S is the water storage rate of the aquifer below the free surface, μ is the gravity water supply of the phreatic aquifer, Γ0 is the upper boundary of the seepage area, that is, the free surface of the groundwater, p is the evaporation and precipitation infiltration intensity of the phreatic surface, h1 is the known boundary water level value, Γ1 is the known water level boundary, Γ2 is the flow boundary of the seepage area, is the normal direction of the boundary surface, K n is the permeability coefficient in the normal direction of the boundary surface, q is the flow rate per unit area of ​​the flow boundary, inflow is positive, outflow is negative, and it is 0 when it is isolated from water, and t is the research time.

[0043] In the conceptual model, the seepage area is represented by Ω, which represents the spatial range of brine. The driving forces of brine seepage mainly include gravitational potential energy and pressure potential energy, which are represented by water level h and pressure head respectively. The seepage capacity of brine in porous media is measured by the permeability coefficient K. Considering the anisotropy of salt lake sediments, the permeability coefficient may be different in the x, y, and z directions, so K is used. x , K y , K z Respectively.

[0044] In addition, the conceptual model also considers the water storage performance of the aquifer, which is characterized by the water supply degree μ and the water storage rate S. During the seepage process, brine may also exchange with the outside world, such as precipitation infiltration recharge, evaporation discharge, etc. These factors can be described by the source and sink term ε. Based on the above parameters, the conceptual model establishes the control equation of brine seepage, that is, the partial differential equation. The left-hand side of the equation represents the storage change of brine, which is equal to the net inflow flux of brine, and the right-hand side represents the net inflow of brine through the seepage flux. In summary, the equation reflects the mass balance relationship in the brine migration process. At the same time, the conceptual model also needs to give appropriate initial conditions and boundary conditions to reflect the hydrogeological conditions of the study area. Among them, the initial condition gives the initial groundwater level distribution h0 of the region, and the boundary conditions include the free surface Γ0 of the groundwater, the known water level boundary Γ1 and the flow boundary Γ2. On the free surface, the brine contacts the atmosphere, evaporation or precipitation infiltration occurs, which is controlled by the parameter p. On the known water level boundary, the water level value h1 is known and can be given directly. On the flow boundary, the flow of brine obeys Darcy's law, which is determined by the normal permeability coefficient K n and source-sink terms q. Since most salt lake brine is stored in unconfined aquifers, the conceptual model distinguishes between unconfined aquifers and confined aquifers. The upper boundary of the unconfined aquifer is the free surface, which is characterized by the gravity water supply μ, while the upper boundary of the confined aquifer is the impermeable roof, which is characterized by the pressure head. By reasonably selecting the model parameters and giving appropriate initial and boundary conditions, the seepage and migration process of brine in porous media can be quantitatively described.

[0045] S202, based on the hydrogeological conceptual model, using visualization software as a water flow simulation platform for salt lake brine in the target area, the target area is divided into a regular rectangular grid of several rows and columns on a plane, and the regular rectangular grid includes a plurality of rectangular units.

[0046] In one example, the seepage migration of brine in porous media is a complex three-dimensional nonlinear process that requires the use of numerical simulation. Meshing is the basis of numerical simulation methods such as finite difference. By discretizing the continuous seepage area into a finite number of grid cells, a water balance equation can be established on each cell, thereby achieving an approximate simulation of the brine flow field in the entire area.

[0047] The specific method of grid generation is to first determine the appropriate grid size according to the spatial scale and calculation accuracy requirements of the study area. The grid size should be small enough to ensure the accuracy of the calculation results, but not too fine to increase the amount of calculation. Under normal circumstances, for salt lake mining areas with a horizontal scale of tens of square kilometers, the grid size can be selected from tens to hundreds of meters. After determining the grid size, the grid generation tools provided by the visualization software, such as Map->2DGrid in GMS or Grid->Define Uniform Grid in Visual Modflow, are used to divide the target area into regular rectangular grids in the x and y horizontal directions according to fixed row and column spacing. Each grid cell is represented by its row and column number, and the row and column numbers are numbered starting from 1. When grid generation, it should be noted that the grid lines should be aligned with the key geological bodies such as the boundaries, faults, and impermeable layers of the study area as much as possible to improve the accuracy of the simulation. If necessary, local encryption can be performed in key areas, that is, a smaller grid size is used, and a larger grid size is used in secondary areas to take into account both calculation accuracy and efficiency. After the grid division is completed, each grid unit represents a discrete seepage sub-area, and its hydrogeological parameters such as permeability coefficient, water supply degree, source and sink terms can be regarded as the average value or representative value within the unit.

[0048] S203, using the row and column numbers of the regular rectangular grid as numbers of the rectangular units to establish a three-dimensional finite difference model.

[0049] In one example, after the grid division has been completed, each grid cell is first numbered. Since a regular rectangular grid is used, each cell can be uniquely identified by its row number i and column number j, that is, the cell in the i-th row and j-th column is recorded as (i, j). In the horizontal direction, the row numbers are 1, 2, ..., m from top to bottom, and the column numbers are 1, 2, ..., n from left to right. This numbering method is simple and intuitive, and is easy for computers to recognize and process.

[0050] After the numbering is completed, considering that salt lake brine is usually stored in the phreatic aquifer, vertical seepage plays an important role in its recharge, so it is necessary to divide the aquifer into layers in the vertical direction. Assuming that the aquifer is divided into l layers in the vertical direction, the entire seepage area can be divided into m×n×l three-dimensional grid units, and each unit can be uniquely identified by its row number i, column number j and layer number k, recorded as (i, j, k). In this way, the precise positioning and description of each grid unit can be achieved through the combination of row number, column number and layer number.

[0051] After establishing the three-dimensional grid system, the partial differential equations in the hydrogeological conceptual model are discretized on each grid unit. Specifically, the partial derivative terms are replaced by the water level difference quotient of adjacent grid points using the central difference format in the finite difference method, and the time step △t and the space step △x, △y, △z are introduced, and the partial differential equations are finally converted into algebraic equations with the grid node water level as the unknown quantity. Since each grid unit satisfies the law of conservation of mass, the water balance equation can be established at the center of the unit, and the balance equations of each unit are combined to form a large sparse matrix. This matrix reflects the water balance relationship of the entire seepage area, and the water level value of each grid node can be obtained by solving the matrix. In the solution process, it is also necessary to combine the given initial conditions and boundary conditions to select an appropriate solution algorithm, such as the catch-up method, the strong implicit process method, etc. Through iterative solution, the dynamic distribution of important parameters such as brine groundwater level, flow velocity, flux, etc. in the study area in three dimensions of space and time can be simulated and calculated.

[0052] Establishing a three-dimensional finite difference model is a key step in carrying out numerical simulation of salt lake brine resources. By discretely solving the hydrogeological conceptual model, the transition from qualitative description to quantitative calculation can be achieved. Using row and column numbers of regular rectangular grids to identify units has the advantages of simple numbering, unique position, and ease of calculation. It is a common method for finite difference modeling. The finite difference model established based on row and column numbering can reasonably describe the three-dimensional seepage and migration process of brine in porous media. The simulation results can intuitively reflect the spatiotemporal distribution law of the dynamic resources of brine, providing an important basis for the rational exploitation of salt lake brine.

[0053] S204, calculating the underground horizontal measurement of each rectangular unit through a three-dimensional finite difference model, and calculating the recharge amount of salt lake brine in the target area according to the underground horizontal measurement of each rectangular unit, and taking the recharge amount as the recharge resource amount of salt lake brine in the target area.

[0054] In one example, underground level measurement reflects the dynamic changes of brine within each unit and is an important indicator for evaluating brine resources, while brine recharge directly determines the exploitable potential of brine resources and is an important basis for formulating mining plans.

[0055] Specifically, the calculation of underground level measurement is based on the solution results of the three-dimensional finite difference model. Through numerical simulation, the groundwater level value of each grid node at different time steps can be obtained. On this basis, Darcy's law is used to calculate the water flux of each grid unit through the six surfaces, that is, the volume flow of brine, and the water balance equation of each unit is established in combination with the source and sink terms of the unit. In a complete calculation period, the water balance equation of each unit is integrated to obtain the underground level measurement of the unit. The positive and negative values ​​of the underground level measurement reflect the increase or decrease of the brine reserves in the unit. When it is positive, it indicates that the brine reserves have increased, and vice versa, it indicates that the brine reserves have decreased. This indicator intuitively reflects the enrichment or loss of brine in each unit and is an important basis for evaluating the distribution of brine resources. Based on the underground level measurement of each unit, the balance of the entire target area is accumulated to estimate the brine recharge on a regional scale. Brine recharge reflects the total amount of brine supplied to the target area through atmospheric precipitation infiltration, river and lake lateral seepage, underground lateral runoff, etc. within a certain period of time. It represents the natural recharge capacity of brine resources and their recovery potential after mining, and is a key indicator for evaluating the exploitability of brine resources.

[0056] Taking the brine replenishment as the brine replenishment resource can provide an important reference for the rational formulation of mining plans. Under normal circumstances, the annual brine mining volume should not exceed its multi-year average replenishment volume to achieve sustainable utilization of brine resources.

[0057] Based on the above embodiment, as an optional implementation, in S204, calculating the underground level of each rectangular unit by using a three-dimensional finite difference model specifically includes: S401, setting hydrogeological parameters in the three-dimensional finite difference model, the hydrogeological parameters including the spatial position, lithology parameters and boundary conditions of the aquifer.

[0058] S402, based on the spatial position, lithological parameters and boundary conditions of the aquifer, the groundwater flow of the salt lake brine is numerically simulated to obtain the underground horizontal measurement of each rectangular unit.

[0059] In one example, the seepage and migration of brine in porous media is a complex three-dimensional nonlinear process, which is affected by many factors such as the lithology of the aquifer and boundary conditions. It is difficult to solve it directly through analytical methods, but requires the help of numerical simulation technology. Numerical simulation can fully consider the spatial heterogeneity of the aquifer and the spatiotemporal changes of boundary conditions on the basis of a three-dimensional finite difference model, quantitatively describe the spatiotemporal evolution of the brine seepage field, predict the dynamic changes of brine under different mining schemes, and provide an important basis for the rational exploitation of salt lake brine.

[0060] In the specific implementation, it is first necessary to accurately describe the spatial position and lithological parameters of the aquifer in the three-dimensional finite difference model. This requires, on the basis of the previous work, integrating the data of borehole cores, geophysical exploration, pumping tests, etc., to interpolate or invert the key parameters such as the spatial distribution, thickness, permeability, and water supply of the aquifer, and assign the results to each grid unit of the three-dimensional model. Since salt lake sediments usually have strong spatial heterogeneity, the aquifer parameters may show large differences between adjacent units. Therefore, when assigning parameters, geostatistical methods such as Kriging should be used to accurately characterize the spatial variation characteristics of the parameters. While accurately describing the properties of the aquifer, it is also necessary to reasonably set the boundary conditions of the model.

[0061] For salt lake brine, the main boundary types include: constant head boundaries, such as areas connected to surface water bodies such as lakes and rivers; constant flow boundaries, such as precipitation infiltration, evaporation and discharge; and water-blocking boundaries, such as zero flux boundaries connected to bedrock. After the parameter assignment and boundary setting are completed, the three-dimensional finite difference model can be used for numerical solution. By discretely solving the partial differential equations in the hydrogeological conceptual model, the groundwater level value of each grid node at any time can be obtained under the premise of satisfying the initial conditions and boundary conditions. On this basis, Darcy's law is used to calculate the actual water flux of each grid unit, and the water balance equation of each unit is time-integrated in combination with the source and sink terms of the unit, and finally the underground level measurement of each unit within the set time period can be obtained.

[0062] S103, the water supply resource volume, the elastic resource volume and the replenishment resource volume are arithmetically added together to obtain the recoverable resource volume of the salt lake brine in the target area.

[0063] In one example, after calculating the storage resources (including water supply resources and elastic resources) and the recharge resources, the three are arithmetically summed to finally obtain the recoverable resources of salt lake brine in the target area.

[0064] The above three resource quantities reflect different mechanisms of brine resource formation. Among them, the water supply resource quantity corresponds to the release process of brine under the action of gravity, the elastic resource quantity corresponds to the discharge process of brine under the action of elastic pressure, and the recharge resource quantity corresponds to the lateral migration and recharge process of brine under dynamic mining conditions. The formation mechanisms of the three resources are relatively independent, so they can be superimposed. Secondly, the pressure resource method mainly considers elastic resources, while the hydrodynamic method mainly considers recharge resources.

[0065] The arithmetic summation of the water-feeding resource volume, the elastic resource volume and the recharge resource volume is a reasonable solution for calculating the recoverable resource volume. In specific implementation, the water-feeding resource volume and the elastic resource volume of the target area can be calculated according to the first formula and the second formula respectively, and then the recharge resource volume can be calculated by the three-dimensional finite difference numerical simulation method, and finally the three resource volumes can be added together to obtain the total recoverable resource volume.

[0066] Based on the above method, the present application also discloses a device for calculating the recoverable resources of salt lake brine, such as Figure 2 As shown, Figure 2 It is a structural schematic diagram of a device for calculating the recoverable resources of salt lake brine provided in an embodiment of the present application.

[0067] A device for calculating the recoverable resource amount of salt lake brine, comprising: a first calculation module, a second calculation module and a generation module; wherein the first calculation module is used to calculate the storage resource amount of salt lake brine in a target area, and the storage resource amount includes the water supply resource amount and the elastic resource amount; the second calculation module is used to calculate the replenishment resource amount of salt lake brine in the target area; the generation module is used to arithmetically add the water supply resource amount, the elastic resource amount and the replenishment resource amount to obtain the recoverable resource amount of salt lake brine in the target area.

[0068] It should be noted that: when the device provided in the above embodiment realizes its function, only the division of the above functional modules is used as an example. In actual application, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the device and method embodiments provided in the above embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

[0069] See also Figure 3 , is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Figure 3 As shown, the electronic device 1000 may include: at least one processor 1001 , at least one network interface 1004 , a user interface 1003 , a memory 1005 , and at least one communication bus 1002 .

[0070] The communication bus 1002 is used to realize the connection and communication between these components.

[0071] The user interface 1003 may include a display screen (Display) and a camera (Camera), and the optional user interface 1003 may also include a standard wired interface and a wireless interface.

[0072] The network interface 1004 may optionally include a standard wired interface or a wireless interface (such as a WI-FI interface).

[0073] Among them, the processor 1001 may include one or more processing cores. The processor 1001 uses various interfaces and lines to connect various parts in the entire server, and executes various functions of the server and processes data by running or executing instructions, programs, code sets or instruction sets stored in the memory 1005, and calling data stored in the memory 1005. Optionally, the processor 1001 can be implemented in at least one hardware form of digital signal processing (Digital, Signal, Processing, DSP), field programmable gate array (Field-Programmable, Gate, Array, FPGA), and programmable logic array (Programmable, Logic, Array, PLA). The processor 1001 can integrate one or more combinations of central processing unit (Central, Processing, Unit, CPU), image processor (Graphics, Processing, Unit, GPU) and modem. Among them, the CPU mainly processes the operating system, user interface and application program, etc.; the GPU is responsible for rendering and drawing the content to be displayed on the display screen; the modem is used to process wireless communication. It can be understood that the above-mentioned modem may not be integrated into the processor 1001, and it can be implemented by a single chip.

[0074] Among them, the memory 1005 may include a random access memory (Random, Access, Memory, RAM) and may also include a read-only memory (Read-Only, Memory). Optionally, the memory 1005 includes a non-transitory, computer-readable, storage, medium. The memory 1005 can be used to store instructions, programs, codes, code sets or instruction sets. The memory 1005 may include a program storage area and a data storage area, wherein the program storage area may store instructions for implementing an operating system, instructions for at least one function (such as a touch function, a sound playback function, an image playback function, etc.), instructions for implementing the above-mentioned various method embodiments, etc.; the data storage area may store data involved in the above-mentioned various method embodiments, etc. The memory 1005 may optionally be at least one storage device located away from the aforementioned processor 1001. As Figure 3 As shown, the memory 1005 as a computer storage medium may include an operating system, a network communication module, a user interface module, and an application program for a method for calculating the recoverable resources of salt lake brine.

[0075] exist Figure 3In the electronic device 1000 shown, the user interface 1003 is mainly used to provide an input interface for the user and obtain data input by the user; and the processor 1001 can be used to call an application program stored in the memory 1005 for a method for calculating the recoverable resources of brine in a salt lake. When executed by one or more processors, the electronic device executes one or more of the methods described in the above embodiments.

[0076] An electronic device readable storage medium stores instructions, which, when executed by one or more processors, enable the electronic device to execute one or more of the methods described in the above embodiments.

[0077] It should be noted that, for the aforementioned method embodiments, for the sake of simplicity, they are all described as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the order of the actions described, because according to the present application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required for the present application.

[0078] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0079] In the several embodiments provided in the present application, it should be understood that the disclosed devices can be implemented in other ways. For example, the device embodiments described above are only schematic, such as the division of the units, which is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some service interfaces, and the indirect coupling or communication connection of devices or units can be electrical or other forms.

[0080] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0081] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0082] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a memory and includes several instructions for a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present application. The aforementioned memory includes: various media that can store program codes, such as USB flash drives, mobile hard drives, magnetic disks or optical disks.

[0083] The above is only an exemplary embodiment of the present disclosure, and the scope of the present disclosure cannot be limited thereto. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure. After considering the specification and practicing the disclosure here, those skilled in the art will easily think of other embodiments of the present disclosure. This application is intended to cover any modification, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the technical field not recorded in the present disclosure. The description and examples are only regarded as exemplary, and the scope and spirit of the present disclosure are defined by the claims.

Claims

1. A method for calculating the recoverable resources of salt lake brine, characterized in that: The method comprises: Calculate the storage resource volume of salt lake brine in the target area, wherein the storage resource volume includes the water supply resource volume and the elastic resource volume; calculate the replenishment resource volume of salt lake brine in the target area; The water supply resource volume, the elastic resource volume and the replenishment resource volume are arithmetically added together to obtain the recoverable resource volume of the salt lake brine in the target area.

2. The method for calculating the recoverable resources of salt lake brine according to claim 1, characterized in that: The calculation of the storage resource volume of salt lake brine in the target area includes: Obtain the area of ​​salt lake brine in the target area, the porosity of the rock of the salt lake brine reservoir, the thickness of the salt lake brine layer, the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine; calculate the water supply resource according to the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir and the thickness of the brine layer; Calculate the elastic resource volume according to the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine; The water supply resource amount and the elastic resource amount are arithmetically added to obtain the storage resource amount of salt lake brine in the target area.

3. The method for calculating the recoverable resources of salt lake brine according to claim 2, characterized in that: The calculation of the water supply resource quantity according to the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir and the thickness of the brine layer includes: Substitute the area of ​​the salt lake brine, the porosity of the rock of the salt lake brine reservoir and the thickness of the brine layer into the first formula to calculate the water supply resource reserves; wherein, The first formula is: In the formula, Q1 is the water resource reserve, A is the area of ​​the salt lake brine, is the porosity of the rock of the salt lake brine reservoir, and M is the thickness of the salt lake brine layer.

4. The method for calculating the recoverable resources of salt lake brine according to claim 2, characterized in that: The elastic resource volume is calculated according to the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine, including: Substitute the elastic water release coefficient of the salt lake brine, the average pressure head elevation of the salt lake brine and the top elevation of the pressure brine reservoir of the salt lake brine into the second formula to calculate the elastic resource reserves; wherein, The second formula is: Q2 = u * (hH)A; In the formula, Q2 is the elastic resource volume, A is the elastic water release coefficient of the salt lake brine, h is the average pressure head elevation of the salt lake brine, H is the top elevation of the pressure brine reservoir of the salt lake brine, and A is the area of ​​the salt lake brine.

5. The method for calculating the recoverable resources of salt lake brine according to claim 1, characterized in that: The calculation of the replenishment resource amount of salt lake brine in the target area includes: Establish a hydrogeological conceptual model of salt lake brine in the target area; According to the hydrogeological conceptual model, the target area is divided into a regular rectangular grid of several rows and columns on a plane using visualization software as a water flow simulation platform for salt lake brine in the target area, wherein the regular rectangular grid includes a plurality of rectangular units; Using the row and column numbers of the regular rectangular grid as numbers of rectangular units to establish a three-dimensional finite difference model; The underground level measurement of each rectangular unit is calculated by the three-dimensional finite difference model, and the replenishment amount of the salt lake brine in the target area is calculated according to the underground level measurement of each rectangular unit, and the replenishment amount is used as the replenishment resource amount of the salt lake brine in the target area.

6. The method for calculating the recoverable resources of salt lake brine according to claim 1, characterized in that: The establishment of a hydrogeological conceptual model of salt lake brine in the target area includes: Where Ω is the seepage area, ε is the source-sink coefficient of the aquifer, h is the groundwater level, h0 is the initial water level distribution of the groundwater system, and K x , K y are the horizontal permeability coefficients of the water-containing medium in the x and y directions, K z is the vertical z permeability coefficient of the water-bearing medium, S is the water storage rate of the aquifer below the free surface, μ is the gravity water supply of the phreatic aquifer, Γ0 is the upper boundary of the seepage area, that is, the free surface of the groundwater, p is the evaporation and precipitation infiltration intensity of the phreatic surface, h1 is the known boundary water level value, Γ1 is the known water level boundary, Γ2 is the flow boundary of the seepage area, is the normal direction of the boundary surface, K n is the permeability coefficient in the normal direction of the boundary surface, q is the flow rate per unit area of ​​the flow boundary, inflow is positive, outflow is negative, and it is 0 when it is isolated from water, and t is the research time.

7. The method for calculating the recoverable resources of salt lake brine according to claim 5, characterized in that: The calculating the underground level measurement of each rectangular unit by the three-dimensional finite difference model includes: Setting hydrogeological parameters in the three-dimensional finite difference model, wherein the hydrogeological parameters include the spatial position, lithological parameters and boundary conditions of the aquifer; According to the spatial position, lithological parameters and boundary conditions of the aquifer, the groundwater flow of the salt lake brine is numerically simulated to obtain the underground horizontal measurement of each rectangular unit.

8. A device for calculating the recoverable resources of salt lake brine, characterized in that: The device comprises: a first calculation module, a second calculation module and a generation module; wherein, The first calculation module is used to calculate the storage resource volume of salt lake brine in the target area, and the storage resource volume includes the water supply resource volume and the elastic resource volume; The second calculation module is used to calculate the replenishment resource amount of salt lake brine in the target area; The generation module is used to arithmetically add the water supply resource volume, the elastic resource volume and the replenishment resource volume to obtain the recoverable resource volume of the salt lake brine in the target area.

9. An electronic device, characterized in that: It includes a processor, a memory, a user interface and a network interface, the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes the method as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: A computer program is stored which can be loaded by a processor and execute the method according to any one of claims 1 to 7.

Images