A method, device, equipment and medium for calculating reservoir seepage

Abstract

The embodiment of the present invention provides a method, device, equipment and medium for calculating the seepage volume of a reservoir. The method includes: performing a three-dimensional finite element discretization operation on a three-dimensional model of a reservoir to be calculated to obtain a discretized three-dimensional model; using a Gaussian point in the discretized three-dimensional model as a first Gaussian point, and obtaining the coordinates and flow velocity of the first Gaussian point; using a Gaussian point in a surface model of a section to be calculated as a second Gaussian point, and obtaining the flow velocity of the second Gaussian point based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point; and obtaining the seepage volume of the section to be calculated based on the surface model and the flow velocity of the second Gaussian point. In this way, the seepage volume of the section to be calculated of the reservoir can be calculated by converting the Gaussian point flow velocity, so as to improve the calculation accuracy of the seepage volume.

Description

Technical Field

[0001] The invention belongs to the technical field of reservoir engineering, and in particular relates to a method, device, equipment and medium for calculating reservoir seepage. Background Art

[0002] With the development of reservoir engineering technology, the internal head of the reservoir is usually higher than the external head. Therefore, under the influence of the head difference, reservoir seepage will occur, that is, seepage from the inside of the reservoir to the outside of the reservoir.

[0003] However, reservoir seepage has certain safety hazards, and it is necessary to grasp the reservoir seepage situation in time to determine whether anti-seepage measures need to be taken. Therefore, how to calculate the reservoir seepage has become an urgent problem to be solved. Summary of the invention

[0004] The present invention provides a method, a device, an electronic device and a readable storage medium for calculating the seepage volume of a reservoir, so as to solve the problem of how to calculate the seepage volume of a reservoir.

[0005] In order to solve the above-mentioned technical problems, the present invention is achieved as follows:

[0006] In a first aspect, the present invention provides a method for calculating reservoir seepage, characterized in that the method comprises:

[0007] Perform three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model;

[0008] Taking the Gaussian point in the discretized three-dimensional model as a first Gaussian point, obtaining the flow velocity of the first Gaussian point, and obtaining the coordinates of the first Gaussian point;

[0009] Taking a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and acquiring the flow velocity of the second Gauss point based on the coordinates and flow velocity of the first Gauss point and the coordinates of the second Gauss point;

[0010] Based on the surface model and the flow velocity of the second Gaussian point, the seepage amount of the cross section to be calculated is obtained.

[0011] Optionally, acquiring the flow velocity of the second Gaussian point based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point includes:

[0012] An interpolation calculation is performed based on the flow velocity of the first Gaussian point, the coordinates of the first Gaussian point in the three-dimensional model, and the coordinates of the second Gaussian point in the three-dimensional model to obtain the flow velocity of the second Gaussian point.

[0013] Optionally, before acquiring the flow velocity of the second Gaussian point based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point, the method further includes:

[0014] Establishing a surface model of the cross section to be calculated in the discretized three-dimensional model;

[0015] Finite element division is performed on the surface model, and based on the divided surface model, coordinates of second Gauss points in the three-dimensional model are obtained, as well as Jacobian determinant values ​​of each of the second Gauss points.

[0016] Optionally, obtaining the seepage volume of the cross section to be calculated based on the surface model and the flow velocity of the second Gauss point includes:

[0017] Determining a unit normal vector of a unit surface in the surface model based on the node coordinates of each unit in the surface model;

[0018] For any second Gauss point, project the flow velocity of the second Gauss point based on the unit normal vector to obtain the flow velocity of the second Gauss point in the direction of the unit normal vector as the normal flow velocity;

[0019] Determine the seepage amount of each unit in the surface model based on the normal flow velocity of each second Gaussian point contained in each unit in the surface model, the Jacobian determinant value of the second Gaussian point and the weight of each second Gaussian point;

[0020] Based on the seepage amount of each unit in the surface model, the seepage amount of the cross section to be calculated is obtained.

[0021] Optionally, taking the Gaussian point in the discretized three-dimensional model as a first Gaussian point and acquiring the flow velocity of the first Gaussian point includes:

[0022] Obtaining the water head value of each node in the discretized three-dimensional model;

[0023] Based on the hydraulic head value of each node, the permeability matrix and the unit shape function, the flow velocity of the first Gaussian point is obtained.

[0024] Optionally, the reservoir to be calculated has a plurality of sections to be calculated, and the method further includes:

[0025] The seepage volumes of the multiple sections to be calculated are added together to obtain the total seepage volume of the reservoir to be calculated through the multiple sections to be calculated.

[0026] In a second aspect, the present invention provides a device for calculating reservoir seepage, characterized in that the device comprises:

[0027] A discrete module is used to perform a three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model;

[0028] A first acquisition module, configured to use the Gaussian point in the discretized three-dimensional model as a first Gaussian point, acquire the flow velocity of the first Gaussian point, and acquire the coordinates of the first Gaussian point;

[0029] A second acquisition module, configured to use a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and acquire the flow velocity of the second Gauss point based on the coordinates and flow velocity of the first Gauss point and the coordinates of the second Gauss point;

[0030] The third acquisition module is used to acquire the seepage volume of the cross section to be calculated based on the surface model and the flow velocity of the second Gaussian point.

[0031] In a third aspect, the present invention provides an electronic device comprising: a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method described in the first aspect when executing the program.

[0032] In a fourth aspect, the present invention provides a readable storage medium, which, when instructions in the storage medium are executed by a processor of an electronic device, enables the electronic device to execute the method described in the first aspect above.

[0033] The method for calculating the reservoir seepage provided by the embodiment of the present invention performs a three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model; uses the Gaussian point in the discretized three-dimensional model as the first Gaussian point, and obtains the coordinates and flow velocity of the first Gaussian point; uses the Gaussian point in the surface model of the section to be calculated as the second Gaussian point, and obtains the flow velocity of the second Gaussian point based on the flow velocity of the first Gaussian point and the coordinates of the second Gaussian point; obtains the seepage of the section to be calculated based on the surface model and the flow velocity of the second Gaussian point. In this way, the seepage of the section to be calculated of the reservoir can be calculated by converting the Gaussian point flow velocity, which can improve the calculation accuracy; and by establishing a surface model, it is not restricted by the grid of the original three-dimensional model of the reservoir, and the density of the interpolation points and the interpolated points can be improved, thereby improving the interpolation accuracy and further ensuring the accuracy of the seepage calculation. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0035] Figure 1 It is a flow chart of the steps of a method for calculating reservoir seepage provided by an embodiment of the present invention;

[0036] Figure 2 is a flow chart of another method for calculating reservoir seepage provided by an embodiment of the present invention;

[0037] Figure 3 The present invention provides a flow chart of a method for calculating reservoir seepage volume.

[0038] Figure 4 It is a structural diagram of a reservoir seepage calculation device provided by an embodiment of the present invention;

[0039] Figure 5 It is a structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0040] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0041] Figure 1 is a flow chart of the steps of a method for calculating reservoir seepage provided by an embodiment of the present invention, such as Figure 1 As shown, the method may include the following steps:

[0042] Step 101: Perform a three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model.

[0043] Step 102: taking the Gaussian point in the discretized three-dimensional model as the first Gaussian point, obtaining the flow velocity of the first Gaussian point, and obtaining the coordinates of the first Gaussian point.

[0044] Step 103: taking a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and obtaining the flow velocity of the second Gauss point based on the coordinates and flow velocity of the first Gauss point and the coordinates of the second Gauss point.

[0045] Step 104: Based on the surface model and the flow velocity of the second Gaussian point, obtain the seepage volume of the cross section to be calculated.

[0046] With respect to the above steps 101 to 104, the embodiment of the present invention can be applied to any reservoir that has a need for seepage calculation, and a reservoir can be selected as the reservoir to be calculated according to actual needs, and the embodiment of the present invention does not limit this. Among them, the seepage refers to the amount of water lost by the reservoir due to leakage, usually refers to the seepage flowing through a building or a dam foundation, which has an important impact on the water storage of the reservoir and is an important indicator of seepage control. Therefore, the analysis and prediction of the seepage is an important operation of reservoir management.

[0047] Among them, the above-mentioned three-dimensional model is built based on the reservoir to be calculated. The collected reservoir size, topography, geology, and material properties can be modeled, meshed, and assigned materials through finite element software to obtain a three-dimensional model corresponding to the reservoir to be calculated, which is convenient for data calculation and processing.

[0048] Furthermore, the above three-dimensional finite element discretization operation can be implemented by conventional finite element modeling methods, that is, the three-dimensional model is divided into a finite number of units, each unit is a simple part, and each unit is composed of nodes and connections between nodes. The discretized three-dimensional model contains multiple units and nodes.

[0049] Among them, the above-mentioned Gaussian points refer to Gaussian integral points. Each unit in the discretized three-dimensional model contains a number of Gaussian points. Gaussian points are a kind of numerical integration points. By numerically integrating the integrand at these points, an approximate integral value can be obtained. Specifically, the number of Gaussian points contained in each unit can be set by itself, and the embodiment of the present invention does not limit this. It should be noted that the more Gaussian points each unit has, the higher the accuracy of data calculation.

[0050] Furthermore, in the embodiment of the present invention, the Gaussian point in the discretized three-dimensional model can be used as the first Gaussian point, and the flow velocity of the first Gaussian point, that is, the flow velocity of the water in the reservoir at the first Gaussian point, can be obtained. And the coordinates of the first Gaussian point, that is, the coordinates of the first Gaussian point in the three-dimensional model, can be obtained.

[0051] Optionally, the flow velocity of the first Gaussian point can be calculated by conventional finite element calculation method, according to the head value of each node in the three-dimensional model, the permeability matrix and the unit shape function, to solve the flow velocity of each Gaussian point in each unit. It should be noted that the flow velocity of the first Gaussian point is a vector at this time. Among them, the above-mentioned permeability matrix is ​​related to the material properties, and the unit shape function is related to the unit shape. After the model mesh is completed and the material of each unit is assigned, it can be obtained by calculation.

[0052] The cross section to be calculated may be any reservoir cross section selected by oneself, may be an equipotential surface, and of course may be a non-equipotential surface, which is not limited in the embodiment of the present invention.

[0053] Furthermore, the surface model is established in the three-dimensional model according to the cross section to be calculated, that is, according to the position of the cross section to be calculated in the reservoir, a corresponding surface model is established in the three-dimensional model and meshing is performed. This step can also be achieved by three-dimensional software.

[0054] Furthermore, after the surface model is obtained, the flow velocity of the second Gauss point can be obtained according to the coordinates of the Gauss points of each unit in the surface model and the flow velocity of the first Gauss point.

[0055] Optionally, acquiring the flow velocity of the second Gaussian point based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point includes:

[0056] An interpolation calculation is performed based on the flow velocity of the first Gaussian point, the coordinates of the first Gaussian point in the three-dimensional model, and the coordinates of the second Gaussian point in the three-dimensional model to obtain the flow velocity of the second Gaussian point.

[0057] Specifically, the above interpolation calculation can be implemented by any interpolation algorithm, that is, based on the relationship between the coordinates of the first Gauss point and the coordinates of the second Gauss point, the flow velocity of the first Gauss point is used for interpolation to obtain the flow velocity corresponding to the position of the second Gauss point. Exemplarily, a linear interpolation algorithm can be used for implementation, and the embodiment of the present invention does not limit the interpolation algorithm used.

[0058] In summary, in the embodiment of the present invention, a three-dimensional finite element discretization operation is performed on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model; the Gaussian point in the discretized three-dimensional model is used as the first Gaussian point, and the coordinates and flow velocity of the first Gaussian point are obtained; the Gaussian point in the surface model of the section to be calculated is used as the second Gaussian point, and the flow velocity of the second Gaussian point is obtained based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point; the seepage volume of the section to be calculated is obtained based on the surface model and the flow velocity of the second Gaussian point. In this way, the seepage volume of the section to be calculated of the reservoir is calculated by converting the Gaussian point flow velocity, which can improve the calculation accuracy; and by establishing a surface model, it is not restricted by the grid of the original three-dimensional model of the reservoir, and the density of the interpolation points and the interpolated points can be improved, thereby improving the interpolation accuracy and further ensuring the accuracy of the seepage volume calculation. The seepage volume can be calculated for any section, which improves the flexibility of the calculation.

[0059] Optionally, before acquiring the flow velocity of the second Gaussian point based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point, the method further includes:

[0060] Establishing a surface model of the cross section to be calculated in the discretized three-dimensional model;

[0061] Finite element division is performed on the surface model, and based on the divided surface model, coordinates of second Gauss points in the three-dimensional model are obtained, as well as Jacobian determinant values ​​of each of the second Gauss points.

[0062] The Jacobian determinant value refers to the determinant value corresponding to the Gaussian point in the Jacobian determinant. Specifically, the corresponding Jacobian matrix can be obtained through the surface model, and then the value corresponding to the second Gaussian point can be obtained from the Jacobian matrix to calculate the Jacobian determinant value.

[0063] Optionally, obtaining the seepage volume of the cross section to be calculated based on the surface model and the flow velocity of the second Gauss point includes:

[0064] Determining a unit normal vector of a unit surface in the surface model based on the node coordinates of each unit in the surface model;

[0065] For any second Gauss point, project the flow velocity of the second Gauss point based on the unit normal vector to obtain the flow velocity of the second Gauss point in the direction of the unit normal vector as the normal flow velocity;

[0066] Determine the seepage amount of each unit in the surface model based on the normal flow velocity of each second Gaussian point contained in each unit in the surface model, the Jacobian determinant value of the second Gaussian point and the weight of each second Gaussian point;

[0067] Based on the seepage amount of each unit in the surface model, the seepage amount of the cross section to be calculated is obtained.

[0068] Optionally, taking the Gaussian point in the discretized three-dimensional model as a first Gaussian point and acquiring the flow velocity of the first Gaussian point includes:

[0069] Obtaining the water head value of each node in the discretized three-dimensional model;

[0070] Based on the hydraulic head value of each node, the permeability matrix and the unit shape function, the flow velocity of the first Gaussian point is obtained.

[0071] The head value mentioned above refers to the total head value, which is used to characterize the energy of the fluid. Specifically, the head value of each node mentioned above can be calculated by the basic equation of stable seepage. Furthermore, the flow velocity of the first Gaussian point can be calculated based on the head value of each node, the permeability matrix and the unit shape function by conventional finite element calculation method.

[0072] Optionally, the reservoir to be calculated has a plurality of sections to be calculated, and the method further includes:

[0073] The seepage volumes of the multiple sections to be calculated are added together to obtain the total seepage volume of the reservoir to be calculated through the multiple sections to be calculated.

[0074] In this way, through the combination of any cross sections, the total seepage volume of the reservoir passing through any combination of cross sections can be calculated, which can improve the practicability of the embodiment of the present invention in engineering.

[0075] Specifically, Figure 2 is a flow chart of another method for calculating reservoir seepage provided by an embodiment of the present invention, such as Figure 2 As shown:

[0076] (1) Solving the three-dimensional seepage field (performing three-dimensional finite element discretization on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model): conventional finite element calculation methods are used to discretize the model (the three-dimensional model of the reservoir to be calculated) in three dimensions, and the head value of each node is solved according to the basic equation of stable seepage.

[0077] (2) Solve the flow velocity of each Gaussian point in the three-dimensional unit: Use the conventional finite element calculation method to solve the flow velocity of each Gaussian point i in each unit based on the water head of each node, the permeability matrix and the unit shape function. , the flow velocity is a vector.

[0078] (3) Determine the calculation section (establish a surface model of the section to be calculated in the discretized three-dimensional model; perform finite element division on the surface model): select any section in the three-dimensional model for which the seepage volume needs to be calculated, establish a surface model of the section in the original model, and perform mesh division.

[0079] (4) Calculating Gaussian point information in the surface model (obtaining the coordinates of the second Gaussian point in the three-dimensional model based on the divided surface model, and obtaining the Jacobian determinant value of each second Gaussian point): Calculating the coordinates of the Gaussian point ig of each unit in the surface model of any cross-section in the three-dimensional model coordinate system and the Jacobian determinant value of each Gaussian point .

[0080] (5) Calculate the velocity of the Gaussian point in the surface model (based on the coordinates and velocity of the first Gaussian point and the coordinates of the second Gaussian point, obtain the velocity of the second Gaussian point): According to the coordinates and velocity of the Gaussian point i in the original model, and the coordinates of the Gaussian point ig in the cross-sectional surface model, and interpolate the flow velocity of the Gaussian point ig in the cross-sectional surface model , the flow velocity is a vector.

[0081] (6) Calculate the total seepage of the surface model (based on the surface model and the flow velocity of the second Gaussian point, obtain the seepage of the section to be calculated; add the seepage of multiple sections to be calculated to obtain the total leakage of the reservoir to be calculated through multiple sections to be calculated): calculate the unit normal vector of the unit surface according to the node coordinates of each unit in the surface model ; Project the flow velocity of each Gaussian point in each unit in the surface model to the normal direction of the unit surface to obtain the flow velocity of each Gaussian point along the normal direction of the unit , ; According to the flow velocity of each Gauss point along the unit normal direction , Jacobian values ​​of each Gaussian point and the weight of each Gaussian point Calculate the seepage flow through each unit , , where ng is the total number of Gaussian points in a unit; the sum of the seepage of all units is the total seepage in the surface model , , which is the total seepage flow through any cross section.

[0082] As another example, Figure 3 FIG. 1 is a flow chart of a method for calculating reservoir seepage volume provided by an embodiment of the present invention. Figure 3 As shown:

[0083] (1) Solving the three-dimensional seepage field: The conventional finite element calculation method is used to discretize the model in three dimensions, and the water head value of each node is solved according to the basic equation of stable seepage.

[0084] (2) Solve the flow velocity of each Gaussian point in the three-dimensional unit: Use the conventional finite element calculation method to solve the flow velocity of each Gaussian point i in each unit based on the water head of each node, the permeability matrix and the unit shape function. , the flow velocity is a vector.

[0085] (3) Select calculation sections: In order to calculate the total seepage of the reservoir from the surrounding banks and the bottom basin to the outside in multiple directions, multiple calculation sections need to be selected. The calculation section can form a semi-enclosed body with an open top and closed sides and bottom.

[0086] (4) Section meshing: Select a calculation section, create a surface model of the section in the original model, and perform meshing.

[0087] (5) Calculate the Gaussian point information in the surface model: Calculate the coordinates of the Gaussian points ig of each unit in the surface model of any cross section in the three-dimensional model coordinate system and the Jacobian determinant value of each Gaussian point .

[0088] (6) Calculate the velocity of the Gaussian point in the surface model: According to the coordinates and velocity of the Gaussian point i in the original model and the coordinates of the Gaussian point ig in the cross-sectional surface model, and interpolate the flow velocity of the Gaussian point ig in the cross-sectional surface model , the flow velocity is a vector.

[0089] (7) Calculate the total seepage of the surface model: Calculate the unit normal vector of the element surface based on the node coordinates of each element in the surface model ; Project the flow velocity of each Gaussian point in each unit in the surface model to the normal direction of the unit surface to obtain the flow velocity of each Gaussian point along the normal direction of the unit , ; According to the flow velocity of each Gauss point along the unit normal direction , Jacobian values ​​of each Gaussian point and the weight of each Gaussian point Calculate the seepage flow through each unit , , where ng is the total number of Gaussian points in a unit; the sum of the seepage of all units is the total seepage in the surface model , , which is the total seepage flow through any cross section.

[0090] (8) Repeat steps (4) to (7) for each calculation section to obtain the seepage volume of each calculation section. The total seepage volume of the reservoir flowing through the semi-enclosed volume can be obtained by summing them up.

[0091] It should be noted that in the related art, the seepage flow is usually calculated by the cross section method, equivalent node flow method, complex three-dimensional seepage field seepage calculation method and groundwater seepage flow calculation based on equipotential surface. Among them, the implementation scheme of the cross section method is: after obtaining the finite element numerical solution of the seepage field head function H(x, y, z), for any eight-node hexahedron isoparametric unit, select the cross section (cross section) formed by the midpoints of the four edges between a pair of faces as the flow section; in two-dimensional problems, the midline is used as the flow section. However, in this method, when the material partition or geological conditions of the calculation area are complex, the unit shape is very irregular, and the cross section is also an extremely irregular distorted surface. When the cross section of the unit is used as the flow surface, the accuracy of the calculated seepage flow is greatly reduced, and sometimes it cannot meet the requirements of engineering applications.

[0092] Among them, the equivalent node flow method is: the seepage on any flow section is expressed as the algebraic sum of the product of the conduction coefficient of the relevant unit and the corresponding node head. However, the selection of the water-passing section is more difficult, and the cross-section seepage obtained is the algebraic sum of the equivalent node flows of the front and rear units. The contribution of the front and rear units needs to be calculated separately during calculation, otherwise the calculation result is 0. Separating the contribution of the front and rear units requires given conditions or manual intervention, which is labor-intensive and prone to errors. The number of nodes covered by the flow surface is difficult to guarantee. When the number of nodes is too small, the calculation accuracy is reduced.

[0093] Among them, the above-mentioned complex three-dimensional seepage field seepage calculation method takes the seepage flow through the two adjacent end faces of a flow tube (the distance tends to zero) as the seepage flow through any of the two sections as the theoretical basis, finds the centroid point of the square grid unit on the calculation section, obtains the local coordinates of the centroid point through iterative calculation, and obtains the node head of any plane discrete grid based on this, and then obtains the seepage flow through any section. However, this method requires assuming the distance between the centroids of the two planes, which theoretically affects the calculation accuracy of the flow; it is necessary to use an iterative method to solve the local coordinates of the unit coordinate system of the unit through which the centroid of the cross-sectional unit passes in the overall area, and the calculation time is relatively long; the calculated head value is the head at the centroid of the unit obtained by dividing any calculation section according to any rule, and some units will happen to be on the material stratification surface, so the division form will affect the calculation accuracy; the accuracy of calculating the seepage flow of the unit based only on the node head at the centroid of the cross-sectional unit is limited; it requires that the cross-sectional mesh be discretized into square grids with equal length and width and small enough, and the adaptability to the cross-sectional shape is weak.

[0094] Among them, the above-mentioned calculation of groundwater seepage flow based on equipotential surface uses conventional finite element method to solve the seepage field and find the head value of each unit node; determine the head value of the equipotential surface, and the equal head nodes constitute the equipotential surface; for any unit, judge the size relationship between the head value of each node and the head value of the equipotential surface, and construct the equipotential surface patch; calculate the seepage flow through each equipotential surface patch, and sum up each patch to get the flow through the specified head equipotential surface. However, this method calculates the seepage flow based on the equipotential surface. Although the flow direction is perpendicular to the calculation section and the calculation is simple, in a complex seepage field, the equipotential surface may have a complex shape and cannot match the section of concern to the project. For example, the project may be more concerned about the seepage flow through a certain section, or the total seepage flow of the reservoir. The seepage flow calculated by this method is difficult to be directly applied to the analysis and evaluation of the seepage degree of the project; for a certain equipotential surface, the number of nodes in the three-dimensional model whose node head is equal to the head of the equipotential surface cannot be guaranteed. When the number of nodes is small, the calculation accuracy decreases.

[0095] The method for calculating the reservoir seepage provided in the embodiment of the present invention can improve the calculation accuracy by converting the Gaussian point flow velocity to calculate the seepage of the cross-section to be calculated of the reservoir; and by establishing a surface model, it is not restricted by the grid of the original three-dimensional model of the reservoir, and the density of the interpolation points and the interpolated points can be improved, thereby improving the interpolation accuracy and further ensuring the accuracy of the seepage calculation; the seepage calculation can be performed on any cross-section, thereby improving the flexibility of the calculation; through the combination of any cross-sections, the total seepage of the reservoir through any combination of cross-sections can be calculated, thereby improving the practicality of the project.

[0096] Figure 4 is a structural diagram of a reservoir seepage calculation device provided by an embodiment of the present invention, such as Figure 4 As shown, the device 20 may include:

[0097] Discretization module 201, used for performing three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model;

[0098] A first acquisition module 202, configured to use the Gaussian point in the discretized three-dimensional model as a first Gaussian point, acquire the flow velocity of the first Gaussian point, and acquire the coordinates of the first Gaussian point;

[0099] A second acquisition module 203 is used to use a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and acquire the flow velocity of the second Gauss point based on the coordinates and flow velocity of the first Gauss point and the coordinates of the second Gauss point;

[0100] The third acquisition module 204 is used to acquire the seepage volume of the cross section to be calculated based on the surface model and the flow velocity of the second Gaussian point.

[0101] Optionally, the second acquisition module is specifically used to:

[0102] An interpolation calculation is performed based on the flow velocity of the first Gaussian point, the coordinates of the first Gaussian point in the three-dimensional model, and the coordinates of the second Gaussian point in the three-dimensional model to obtain the flow velocity of the second Gaussian point.

[0103] Optionally, the device further comprises:

[0104] An establishing module, used for establishing a surface model of the section to be calculated in the discretized three-dimensional model;

[0105] A partitioning module is used to perform finite element partitioning on the surface model, and based on the partitioned surface model, obtain the coordinates of the second Gaussian points in the three-dimensional model, and obtain the Jacobian determinant value of each of the second Gaussian points.

[0106] Optionally, the third acquisition module includes:

[0107] A first determination submodule, configured to determine a unit normal vector of a unit surface in the surface model based on the node coordinates of each unit in the surface model;

[0108] a projection submodule, configured to project the flow velocity of any second Gaussian point based on the unit normal vector to obtain the flow velocity of the second Gaussian point in the direction of the unit normal vector as the normal flow velocity;

[0109] A second determination submodule is used to determine the seepage amount of each unit in the surface model based on the normal flow velocity of each second Gaussian point contained in each unit in the surface model, the Jacobian determinant value of the second Gaussian point and the weight of each second Gaussian point;

[0110] The fourth acquisition submodule is used to acquire the seepage amount of the cross section to be calculated based on the seepage amount of each unit in the surface model.

[0111] Optionally, the first acquisition module includes:

[0112] A fifth acquisition submodule is used to obtain the water head value of each node in the discretized three-dimensional model;

[0113] The sixth acquisition submodule is used to acquire the flow velocity of the first Gaussian point based on the head value of each node, the permeability matrix and the unit shape function.

[0114] Optionally, the reservoir to be calculated has a plurality of sections to be calculated, and the device further comprises:

[0115] The adding module is used to add the seepage volume of the multiple sections to be calculated to obtain the total seepage volume of the reservoir to be calculated through the multiple sections to be calculated.

[0116] In summary, the reservoir seepage calculation device provided by the embodiment of the present invention performs a three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model; the Gaussian point in the discretized three-dimensional model is used as the first Gaussian point, and the coordinates and flow velocity of the first Gaussian point are obtained; the Gaussian point in the surface model of the section to be calculated is used as the second Gaussian point, and the flow velocity of the second Gaussian point is obtained based on the coordinates and flow velocity of the first Gaussian point and the coordinates of the second Gaussian point; the seepage of the section to be calculated is obtained based on the surface model and the flow velocity of the second Gaussian point. In this way, the seepage of the section to be calculated of the reservoir is calculated by converting the Gaussian point flow velocity, which can improve the calculation accuracy; and by establishing a surface model, it is not restricted by the grid of the original three-dimensional model of the reservoir, and the density of the interpolation points and the interpolated points can be improved, thereby improving the interpolation accuracy and further ensuring the accuracy of the seepage calculation; the seepage calculation can be performed on any section to improve the flexibility of the calculation; through the combination of any sections, the total seepage of the reservoir through any combination of sections can be calculated, which can improve the engineering practicality.

[0117] The present invention also provides an electronic device, see Figure 5 , including: a processor 301, a memory 302, and a computer program 3021 stored in the memory and executable on the processor, wherein when the processor executes the program, the method for calculating the reservoir seepage volume of the aforementioned embodiment is implemented.

[0118] The present invention also provides a readable storage medium. When the instructions in the storage medium are executed by a processor of an electronic device, the electronic device can execute the method for calculating the reservoir seepage volume of the aforementioned embodiment.

[0119] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0120] The algorithm and display provided herein are not inherently related to any particular computer, virtual system or other device. Various general purpose systems can also be used together with the teachings based on this. According to the above description, it is obvious that the structure required for constructing such systems. In addition, the present invention is not directed to any specific programming language either. It should be understood that various programming languages ​​can be utilized to realize the content of the present invention described herein, and the description of the above specific languages ​​is for disclosing the best mode of the present invention.

[0121] In the description provided herein, a large number of specific details are described. However, it is understood that embodiments of the present invention can be practiced without these specific details. In some instances, well-known methods, structures and techniques are not shown in detail so as not to obscure the understanding of this description.

[0122] Similarly, it should be understood that in order to streamline the present invention and aid in understanding one or more of the various inventive aspects, in the above description of exemplary embodiments of the present invention, the various features of the present invention are sometimes grouped together into a single embodiment, figure, or description thereof. However, this disclosed method should not be interpreted as reflecting the following intention: that the claimed invention requires more features than the features explicitly recited in each claim. More specifically, as reflected in the claims below, inventive aspects lie in less than all the features of the individual embodiments disclosed above. Therefore, the claims that follow the specific embodiment are hereby expressly incorporated into the specific embodiment, with each claim itself serving as a separate embodiment of the present invention.

[0123] Those skilled in the art will appreciate that the modules in the devices in the embodiments may be adaptively changed and arranged in one or more devices different from the embodiments. The modules or units or components in the embodiments may be combined into one module or unit or component, and further may be divided into a plurality of submodules or subunits or subcomponents. All features disclosed in this specification (including the accompanying claims, abstracts and drawings) and all processes or units of any method or device so disclosed may be combined in any combination, except that at least some of such features and / or processes or units are mutually exclusive. Unless otherwise expressly stated, each feature disclosed in this specification (including the accompanying claims, abstracts and drawings) may be replaced by an alternative feature that provides the same, equivalent or similar purpose.

[0124] The various component embodiments of the present invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. It should be understood by those skilled in the art that a microprocessor or a digital signal processor (DSP) may be used in practice to implement some or all of the functions of some or all of the components in the sorting device according to the present invention. The present invention may also be implemented as a device or apparatus program for executing part or all of the methods described herein. Such a program for implementing the present invention may be stored on a computer-readable medium, or may be in the form of one or more signals. Such a signal may be downloaded from an Internet website, or provided on a carrier signal, or provided in any other form.

[0125] It should be noted that the above embodiments illustrate the present invention rather than limit it, and that those skilled in the art may devise alternative embodiments without departing from the scope of the appended claims. In the claims, any reference symbol between brackets shall not be construed as a limitation on the claims. The word "comprising" does not exclude the presence of elements or steps not listed in the claims. The word "one" or "an" preceding an element does not exclude the presence of a plurality of such elements. The present invention may be implemented by means of hardware comprising a number of different elements and by means of a suitably programmed computer. In a unit claim enumerating a number of devices, several of these devices may be embodied by the same hardware item. The use of the words first, second, and third, etc., does not indicate any order. These words may be interpreted as names.

[0126] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

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

[0128] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed by the present invention, which should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for calculating reservoir seepage, characterized in that: The method comprises: Perform three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model; Taking the Gaussian point in the discretized three-dimensional model as a first Gaussian point, obtaining the flow velocity of the first Gaussian point, and obtaining the coordinates of the first Gaussian point; Taking a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and performing interpolation calculation based on the flow velocity of the first Gauss point, the coordinates of the first Gauss point in the three-dimensional model, and the coordinates of the second Gauss point in the three-dimensional model to obtain the flow velocity of the second Gauss point; Determining a unit normal vector of a unit surface in the surface model based on the node coordinates of each unit in the surface model; For any second Gaussian point, project the flow velocity of the second Gaussian point based on the unit normal vector to obtain the flow velocity of the second Gaussian point in the direction of the unit normal vector as the normal flow velocity; Determine the seepage amount of each unit in the surface model based on the normal flow velocity of each second Gaussian point contained in each unit in the surface model, the Jacobian determinant value of the second Gaussian point and the weight of each second Gaussian point; Based on the seepage amount of each unit in the surface model, the seepage amount of the cross section to be calculated is obtained.

2. The method according to claim 1, characterized in that Before performing interpolation calculation based on the flow velocity of the first Gaussian point, the coordinates of the first Gaussian point in the three-dimensional model, and the coordinates of the second Gaussian point in the three-dimensional model to obtain the flow velocity of the second Gaussian point, the method further includes: Establishing a surface model of the cross section to be calculated in the discretized three-dimensional model; Finite element division is performed on the surface model, and based on the divided surface model, coordinates of second Gauss points in the three-dimensional model are obtained, as well as Jacobian determinant values ​​of each of the second Gauss points.

3. The method according to claim 1, characterized in that: The step of taking the Gaussian point in the discretized three-dimensional model as a first Gaussian point and obtaining the flow velocity of the first Gaussian point includes: Obtaining the water head value of each node in the discretized three-dimensional model; Based on the hydraulic head value of each node, the permeability matrix and the unit shape function, the flow velocity of the first Gaussian point is obtained.

4. The method according to any one of claims 1 to 3, characterized in that: The reservoir to be calculated has a plurality of sections to be calculated, and the method further includes: The seepage volumes of the multiple sections to be calculated are added together to obtain the total seepage volume of the reservoir to be calculated through the multiple sections to be calculated.

5. A reservoir seepage calculation device, characterized in that: The device comprises: A discrete module is used to perform a three-dimensional finite element discretization operation on the three-dimensional model of the reservoir to be calculated to obtain a discretized three-dimensional model; A first acquisition module, configured to use the Gaussian point in the discretized three-dimensional model as a first Gaussian point, acquire the flow velocity of the first Gaussian point, and acquire the coordinates of the first Gaussian point; A second acquisition module, configured to use a Gauss point in the surface model of the cross section to be calculated as a second Gauss point, and acquire the flow velocity of the second Gauss point based on the coordinates and flow velocity of the first Gauss point and the coordinates of the second Gauss point; A third acquisition module, configured to acquire the seepage volume of the cross section to be calculated based on the surface model and the flow velocity of the second Gaussian point; The second acquisition module is specifically used for: performing interpolation calculation based on the flow velocity of the first Gaussian point, the coordinates of the first Gaussian point in the three-dimensional model, and the coordinates of the second Gaussian point in the three-dimensional model to obtain the flow velocity of the second Gaussian point; The third acquisition module includes: A first determination submodule, configured to determine a unit normal vector of a unit surface in the surface model based on the node coordinates of each unit in the surface model; a projection submodule, configured to project the flow velocity of any second Gaussian point based on the unit normal vector to obtain the flow velocity of the second Gaussian point in the direction of the unit normal vector as the normal flow velocity; A second determination submodule is used to determine the seepage amount of each unit in the surface model based on the normal flow velocity of each second Gaussian point contained in each unit in the surface model, the Jacobian determinant value of the second Gaussian point and the weight of each second Gaussian point; The fourth acquisition submodule is used to acquire the seepage amount of the cross section to be calculated based on the seepage amount of each unit in the surface model.

6. The device according to claim 5, characterized in that The device also includes: An establishing module, used for establishing a surface model of the section to be calculated in the discretized three-dimensional model; A partitioning module is used to perform finite element partitioning on the surface model, and based on the partitioned surface model, obtain the coordinates of the second Gaussian points in the three-dimensional model, and obtain the Jacobian determinant value of each of the second Gaussian points.

7. The device according to claim 5, characterized in that The first acquisition module includes: A fifth acquisition submodule is used to obtain the water head value of each node in the discretized three-dimensional model; The sixth acquisition submodule is used to acquire the flow velocity of the first Gaussian point based on the head value of each node, the permeability matrix and the unit shape function.

8. The device according to any one of claims 5 to 7, characterized in that: The reservoir to be calculated has a plurality of sections to be calculated, and the device further comprises: The adding module is used to add the seepage volume of the multiple sections to be calculated to obtain the total seepage volume of the reservoir to be calculated through the multiple sections to be calculated.

9. An electronic device, characterized in that: include: A processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method according to any one of claims 1 to 4 when executing the program.

10. A readable storage medium, characterized in that: When the instructions in the storage medium are executed by a processor of an electronic device, the electronic device is enabled to execute the method according to any one of claims 1 to 4.

Images