Burial depth design method and device based on upper limit theorem fault optimization model and medium

Abstract

The invention relates to the technical field of compressed air energy storage, in particular to a burial depth design method and device based on an upper limit theorem fault optimization model and a medium, and the burial depth design method comprises the following steps: S1, determining design parameters of a gas storage; s2, the design parameters are substituted into a single-hole upper limit theorem limit anti-lifting model to calculate the burial depth; s3, determining the horizontal projection length of the fracture surface corresponding to the burial depth in the S2 when the fracture surface extends to the earth surface; s4, whether the fault intersects with the overlying rock mass of the gas storage is judged, and if not, the burial depth in the S2 serves as the safe burial depth; and if yes, calculating the burial depth by adopting an upper limit theorem limit lift-resistant fault optimization model. The method can solve the problems that an existing model cannot consider the influence of fault weakening and the fracture surface at the fault intersection is not estimated accurately, meanwhile, the critical horizontal distance parameter can evaluate the position relation between the fault and the overlying rock mass and calculate the size of the fault influence area, and the critical horizontal distance parameter can also serve as the judgment basis of the fault transverse distance in gas storage site selection. And the burial depth calculation result is more reasonable and reliable.

Description

Technical Field

[0001] This invention relates to the field of compressed air energy storage technology, and in particular to a method, equipment and medium for burial depth design based on an upper bound theorem fault optimization model. Background Technology

[0002] In compressed air energy storage systems, the engineering design of artificial gas storage caverns must first determine their burial depth. Burial depth has a critical impact on the structural safety of the cavern: when the burial depth is shallow, the ground surface above the cavern may heave or deform under the action of high-pressure gas inside, thus threatening the overall structural stability. Furthermore, the burial depth directly relates to the feasibility of site selection; the geological structure and groundwater conditions of the area will further affect the safety and reliability during construction and operation. On the other hand, the selection of the burial depth of the gas storage cavern also relates to the engineering scale of construction access and auxiliary caverns. For the entire compressed air energy storage power station, the cost of these auxiliary projects will significantly affect the total project investment and expected returns, therefore a comprehensive evaluation is required during the design phase.

[0003] Existing burial depth theoretical models assume that the overlying rock mass is a complete, homogeneous, and continuous medium. However, when encountering special geological conditions such as faults, due to limitations such as site conditions and location selection in actual engineering projects, even after ensuring the safe distance between the gas storage tank and the fault, the surrounding rock within the failure surface may still intersect with the fault. The stress transmission path in the intersection area changes. If the calculation is still performed according to the original assumption of the failure surface, the burial depth calculation will be insufficient. Existing burial depth calculation theories cannot consider the adverse effects of faults weakening the model. Summary of the Invention

[0004] The purpose of this invention is to overcome the technical problem that existing calculation models cannot consider the influence of faults and the inaccurate prediction of the fracture surface at the fault intersection leads to insufficient burial depth calculation, which affects the safety of the burial depth of the cavern. The invention provides a burial depth design method, equipment and medium based on a fault optimization model with upper bound theorem.

[0005] In a first aspect, the present invention provides a method for designing burial depth based on an optimization model for limiting uplift faults using an upper bound theorem, comprising: S1: Determine the design parameters of the gas storage facility, including the inner radius R of the top sphere of the gas storage facility and the unit weight of the overlying rock mass of the gas storage facility. The design internal pressure P, initial burial depth H0, and uniaxial compressive strength of the intact rock in the gas storage facility are specified. ci Empirical constant A, empirical constant B, safety factor F s Fault dip angle The horizontal distance L between the fault initiation point and the center of the gas storage facility; the vertical distance D between the fault initiation point and the center of the gas storage facility dome; construct the failure function curve equation f(x, a0, b0, A, B, ... , ci ); S2: Substitute the design parameters into the single-tunnel upper limit theorem limit uplift resistance model to calculate the burial depth H1; S3: Determine the horizontal projection length L0 of the rupture surface corresponding to the burial depth H1 when it extends to the ground surface; S4: Determine whether the fault intersects with the overlying rock mass of the gas storage. If not, take the burial depth H1 in S2 as the safe burial depth H; if so, use the upper limit theorem limit uplift fault optimization model to calculate the burial depth H2. The calculation steps for the burial depth H2 include: S41: Take the burial depth H1 in S2 as the initial burial depth H0; S42: Solve for the volume V1 of the intersection region of the fault, calculate the volume V of the overlying rock mass in the fault optimization model based on V1, and then establish an equilibrium relationship based on V to solve for the ultimate internal pressure P. u : ; S43: Compare the ultimate internal pressure P u The accuracy of the design internal pressure P is judged to see if it meets the requirements. If it does, the initial burial depth H0 is used as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S42~S43 is repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H2.

[0006] This invention proposes a burial depth design method based on a fault optimization model with an upper bound theorem. This model, based on the upper bound theorem of limit analysis, derives a cavern uplift failure function under the assumption that the rock mass obeys the associated Hoek-Brown criterion. The failure curve is rotated along the y-axis to form a three-dimensional failure surface. When encountering adverse geological conditions at a fault, the original failure surface undergoes a sharp change at the fault, extending along the fault strike. By solving the balance relationship between soil gravity and internal pressure within the failure surface after fault cutting, the safe burial depth is obtained. This method solves the problems of existing calculation models failing to consider the impact of fault weakening and inaccurate prediction of the fracture surface at fault intersections. Simultaneously, the critical horizontal distance parameter can assess the positional relationship between the fault and the overlying rock mass of the gas storage facility, as well as the size of the fault-affected area. It can also serve as a basis for determining the lateral distance from the fault in gas storage facility site selection. By considering the impact of fault weakening, the burial depth calculation results are more reasonable and reliable.

[0007] Preferably, in S2, the step of calculating the burial depth H1 of the single-hole tunnel using the single-hole upper limit theorem limit uplift resistance model includes: S21: Determine the initial burial depth ; S22: Determine the destruction function constant b0, b0 = R + H0; S23: Set the boundary conditions f(x, a0, b0, A, B, ... , ci )=0, x=R, and the failure function constant b0, the unit weight of the rock The uniaxial compressive strength of the intact rock ci Substituting the empirical constants A and B into the equation f(x, a0, b0, A, B), , ci Solve for the failure function constant a0; S24: Substitute the destruction function constants a0 and b0 into the destruction function curve equation f(x, a0, b0, A, B, ... , ci The destruction function f(x) is obtained; the horizontal distance L0 from which the destruction function extends to the Earth's surface is calculated, where L0 = a0 / ; S25: Solve for the volume V0 of the damaged breadth soil: ; S26: Solving for the ultimate internal pressure P u : ; S27: Compare the ultimate internal pressure P u The accuracy of the design internal pressure P is judged to see if it meets the requirements. If it does, the initial burial depth H0 is used as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S22~S27 is repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H1.

[0008] Preferably, in S43, determining whether the accuracy meets the requirements includes: like Then the ultimate internal pressure P u The designed internal pressure P meets the accuracy requirements; if Then H0 = H0 + adjustment value; In S27, determining whether the accuracy meets the requirements includes: S271: If If the required accuracy is met; if Then proceed with S272; S272: Compare the ultimate internal pressure P u and the designed internal pressure P; S273: If P > P u If H0 = H0 + adjustment value; if P < P u Then H0 = H0 - adjustment value; Wherein, the allowable error = the design internal pressure P × error ratio, and the error ratio ≤ 5%; the adjustment value = the initial burial depth H0 / 5 × the error ratio.

[0009] Preferably, in step S42, the solution step for the volume V of the overlying rock mass in the fault optimization model includes: C421: The volume V0 of the damaged envelope soil is calculated by following the steps from S22 to S25 based on the initial burial depth H0 described in S41. C422: The volume V of the overlying rock mass in the fault optimization model = the volume V0 of the damaged envelope soil mass - the volume V1 of the intersection region of the fault, that is: .

[0010] Preferably, in step S42, the step of solving for the volume V1 of the intersection region of the fault includes: S421: Establish a coordinate system with the center of the gas storage dome as the origin. First, determine the coordinates of the fault starting point A1 as (x1, y1), the coordinates of the intersection of the fault and the ground B1 as (x2, y2), the coordinates of the starting point of the overlying rock fracture surface A2 as (x3, y3), and the coordinates of the intersection of the fracture surface and the ground B2 as (x4, y4); then calculate the coordinates of the intersection point c1 of the fault and the fracture surface as (x0, y0). S422: Calculate the height h1 of the intersection region, the width b1 of the cone base of the intersection region, and the chord length l1 of the cone base of the intersection region, respectively. h1 = y2 - y0; b1 = x4 - x2; ; S423: Calculate the cone-shaped base area S1 of the intersection region: ; S424: Calculate the volume V1 of the intersection region of the fault: .

[0011] Preferably, in S421, the calculation formulas for the coordinates (x1, y1) of the fault initiation point A1, the coordinates (x2, y2) of the fault-ground intersection point B1, the coordinates (x3, y3) of the initiation point A2 of the overlying rock mass fracture surface, and the coordinates (x4, y4) of the fracture surface-ground intersection point B2 include: x1=L, y1=B, x2=(H1+R-B) / tan +L, y2=y4=H1+R, x3=R, y3=0, x4=L0.

[0012] Preferably, in S421, the calculation steps for the coordinates (x0, y0) of the intersection point c1 of the fault and the rupture surface include: S4211: Constructing the fault line equation: ; S4212: Combining the fault line equation and the failure function curve equation f(x, a0, b0, A, B, ... , ci x0 and y0 are obtained by solving using software; where, ; ; .

[0013] Preferably, in step S4, the step of determining whether the fault intersects with the overlying rock mass of the gas storage tank includes: C41: Project the three-dimensional upper bound theorem limit uplift fault optimization model to form a two-dimensional fault optimization model, such that the two-dimensional fault optimization model is along the axis of the gas storage tank and perpendicular to the fault. C42: Based on the aforementioned two-dimensional fault optimization model, the equation of the failure function curve f(x, a0, b0, A, B, ...) is... , ci Differentiating f'(x, a0, b0, A, B) yields f'(x, a0, b0, A, B). , ci Substitute the interval of independent variables (0, L0) into f'(x, a0, b0, A, B, ...); , ci ), thus obtaining the maximum slope k of the failure surface. r k r =max(f'(x,a0,b0,A,B, , ci The maximum dip angle of the failure surface was calculated. ', =arctan(k) r ); C43: When the fault dip angle ≥ When there is no tangent point between the fault and the overlying rock mass, the critical intersection point M(L0, H+R) between the fault and the overlying rock mass is defined; when the fault dip angle is... < At that time, there is a tangent point between the fault and the failure surface of the overlying rock mass of the gas storage reservoir. The tangent point N(x) is defined as follows: c ,yc ), where x c =solve(f'(x c )=tan y c =f(x c ); C44: Determine the critical horizontal distance L r When the fault dip angle ≥ 'hour, When the fault dip angle < 'hour, ; C45: Compare with L r and L, if L r If L ≥ L, then the fault and the overlying rock mass of the gas storage are mutually opposed; if L r If the fault is less than L, then the fault intersects with the overlying rock mass of the gas storage facility.

[0014] In a second aspect, the present invention provides an electronic device comprising a memory and at least one processor, the memory storing a computer program, the processor executing the computer program to implement the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model as described above.

[0015] In a third aspect, the present invention provides a computer-readable storage medium storing a computer program that, when executed, implements the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model as described above.

[0016] Compared with the prior art, the beneficial effects of the present invention are as follows: This invention proposes a method, equipment, and medium for burial depth design based on an upper-boundary fault optimization model. This model is a three-dimensional rotating body model, which better conforms to the stress mode of large tank-type gas storage facilities. The spatial characteristics of this model make it more sensitive to adverse geological conditions such as faults, thus requiring consideration of the weakening effect of faults on the surrounding rock within the failure surface. This method determines the critical tangent point (or intersection point) by solving for the maximum slope of the failure surface. A critical horizontal distance is introduced to define the abscissa of the fault initiation point when the overlying rock mass and the fault are exactly tangent (intersecting). Based on this parameter, the positional relationship between the fault and the overlying rock mass of the gas storage facility can be quickly determined, thus serving as a parameter for avoidance conditions during the initial site selection of the gas storage facility. Simultaneously, by comparing the actual fault initiation coordinates with the critical horizontal distance, the size of the intersection area can be calculated, thereby further evaluating the adverse effects of the fault on the gas storage facility. When a fault intersects with the surrounding rock within the overlying failure surface of a gas storage facility, the fault optimization calculation model for large tank gas storage facilities recommended in this method should be adopted. By correcting the fracture angle of the surrounding rock in the fault intersection area to the fault dip angle, a more realistic stress and failure mode can be restored, making the calculation results more realistic and reliable. Taking the calculation example 1 in the embodiment as an example, after considering the fault conditions, the calculated burial depth increased from 73.5 to 83.79, an increase of 14%, and the correction result is significant. Attached Figure Description

[0017] Figure 1 This is a three-dimensional schematic diagram of the single-hole upper bound theorem limit resistance model of the present invention.

[0018] Figure 2 This is a planar schematic diagram of the single-hole upper bound theorem limit resistance model of the present invention.

[0019] Figure 3 This is a flowchart illustrating the establishment of the destruction function curve equation of the present invention.

[0020] Figure 4 This is a flowchart of the solution process for the single-hole upper bound theorem limit resistance model of the present invention.

[0021] Figure 5 This is a three-dimensional schematic diagram of the upper bound theorem limit uplift fault optimization model of the present invention.

[0022] Figure 6 This is a three-dimensional schematic diagram of the fault cutting the surrounding rock mass in the upper bound theorem limit uplift fault optimization model of the present invention.

[0023] Figure 7 This is a planar schematic diagram of the upper bound theorem limit uplift fault optimization model of the present invention.

[0024] Figure 8 This is a flowchart of the solution process for the upper bound theorem-based optimization model for uplift faults in this invention.

[0025] Marked in the image: 1. Gas storage facility; 2. Overlying rock mass; 3. Fault; 4. Fault-cutting body. Detailed Implementation

[0026] The present invention will now be described in further detail with reference to specific embodiments. However, this should not be construed as limiting the scope of the present invention to the following embodiments; all technologies implemented based on the content of the present invention fall within the scope of the present invention.

[0027] Example 1 This embodiment provides a method for designing burial depth based on an optimization model for limiting uplift faults using an upper bound theorem.

[0028] refer to Figure 8 The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model in this embodiment includes the following steps: S1: Determine the design parameters of the gas storage facility, including the inner radius R (in meters) of the top sphere of the gas storage facility and the unit weight of the overlying rock mass. (Unit: kN / m) 3 The design internal pressure P (in kPa) of the gas storage facility, the initial burial depth H0 (in m), and the uniaxial compressive strength of intact rock are also considered. ci (Unit: kPa) Rock tensile strength t (Unit: kPa), Empirical constant A, Empirical constant B, Safety factor F s Fault dip angle (Unit: degrees), horizontal distance L (unit: m) from the fault initiation point to the center of the gas storage facility, vertical distance D (unit: m) from the fault initiation point to the center of the gas storage facility dome; construct the failure function curve equation f(x, a0, b0, A, B, , ci ); S2: Substitute the design parameters into the single-tunnel upper limit theorem limit uplift resistance model to calculate the burial depth H1; S3: Determine the horizontal projection length L0 of the rupture surface corresponding to the burial depth H1 when it extends to the ground surface; S4: Determine whether the fault intersects with the overlying rock mass of the gas storage. If not, use the burial depth H1 in S2 as the safe burial depth H; if so, use the upper bound theorem limit uplift fault optimization model to calculate the burial depth H2. The calculation steps for burial depth H2 include: S41: Use the burial depth H1 in S2 as the initial burial depth H0; S42: Solve for the volume V1 of the fault intersection region, calculate the volume V of the overlying rock mass in the fault optimization model based on V1, and then establish an equilibrium relationship based on V to solve for the ultimate internal pressure P. u : ; S43: Compare the ultimate internal pressure P u The system determines whether the accuracy of the design internal pressure P meets the requirements. If it does, the initial burial depth H0 is set as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S42~S43 is repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H2.

[0029] This invention proposes a burial depth design method based on a fault optimization model with an upper bound theorem. This model, based on the upper bound theorem of limit analysis, derives a cavern uplift failure function under the assumption that the rock mass obeys the associated Hoek-Brown criterion. The failure curve is rotated along the y-axis to form a three-dimensional failure surface. When encountering adverse geological conditions at a fault, the original failure surface undergoes a sharp change at the fault, extending along the fault dip angle. By solving the balance relationship between soil gravity and internal pressure within the failure surface after fault cutting, the safe burial depth is obtained. This method solves the problems of existing calculation models failing to consider the impact of fault weakening and inaccurate prediction of the fracture surface at fault intersections. Simultaneously, the critical horizontal distance parameter can assess the positional relationship between the fault and the overlying rock mass of the gas storage facility, calculate the size of the fault-affected area, and can also serve as a basis for determining the lateral distance from the fault in gas storage facility site selection. By considering the impact of fault weakening, the burial depth calculation results are more reasonable and reliable.

[0030] In step S1 above, when determining the relevant design parameters of the gas storage facility, the following can be constructed: Figure 1 The single-hole upper bound theorem limit uplift resistance model shown and Figure 5 The upper bound theorem limit uplift fault optimization model shown is as follows: Figure 2 Then it is Figure 1 The two-dimensional planar model corresponding to the three-dimensional model shown is Figure 7 for Figure 5 The two-dimensional planar model corresponding to the three-dimensional model shown is Figure 6 From Figure 5 The three-dimensional model of the overlying rock mass cut off by the fault (i.e., fault cut body 4) is a model of the gas storage tank 1 and its overlying rock mass 2. The dimensions and other parameters are all in the model. Figure 1 , Figure 2 , Figure 5 and Figure 7 The calculation units for each design parameter have been marked in step S1 above. When performing calculations using the formulas in this embodiment, the calculation units in S1 must be followed. If you need to switch units, the units should be consistent; otherwise, calculation errors will occur.

[0031] The main idea of ​​the burial depth design method of this invention is to simulate the actual cavern structure and its overlying rock mass by constructing three-dimensional and two-dimensional model diagrams. By calculating the gravity of the overlying rock mass (i.e., the volume of the overlying rock mass multiplied by its unit weight) and the internal pressure of the cavern, the equilibrium relationship between the volume of the overlying rock mass and the internal pressure of the cavern is derived. The accuracy is then solved and verified to finally obtain a suitable cavern burial depth. For example... Figure 1 The single-cavity upper bound theorem limit uplift model structure shown is based on the assumption that the overlying rock mass 2 above the gas storage 1 is a three-dimensional structure with the fracture surface of the function f(x). Therefore, it is necessary to consider the corresponding design parameters (such as the initial burial depth H0, the inner radius R of the top sphere of the gas storage 1, and the fracture angle). The volume of the overlying rock mass 2 is calculated, and the balance relationship between the volume of the overlying rock mass 2 and the pressure inside the cavern is derived. Then, the accuracy is solved and verified. After iterative calculation, the final safe burial depth is obtained.

[0032] For example Figure 5 and Figure 7 The upper bound theorem-based optimization model structure for uplift faults, when calculating the volume of the overlying rock mass 2, requires subtracting the fault-cutting body 4 (i.e., the volume cut off by fault 3) from the volume of the three-dimensional frustum-shaped overlying rock mass 2. Figure 6 The volume of the overlying rock mass 2 (shown) is then subtracted from the volume of the fault-cutting body 4 to establish a new equilibrium relationship with the cavern pressure. This can be understood as, after considering the fault's influence, the equilibrium relationship is... Figure 1 The model shown has been optimized, and the safe burial depth calculated based on the optimized balance relationship is more accurate.

[0033] In step S1 above, the constructed destruction function curve equation f(x, a0, b0, A, B, , ci )for: .

[0034] refer to Figure 3 The equation of the function curve f(x, a0, b0, A, B) is broken. , ci The derivation steps of ) include: S11: Derive the total internal energy power W consumed by the rock mass on the failure function curve. in Derivation of the power W produced by gravity g ; S12: Power W based on internal pressure p The expression W p =W in -W g =p u ×u, the ultimate internal pressure p is derived. uThe expression for u, where u is the velocity field (m / s²). 2 ); S13: Based on the ultimate internal pressure p u The expression is derived through Lagrange transformation to obtain the equation of the destruction function curve f(x, a0, b0, A, B, γ, σ). ci ).

[0035] In step S11 above, the total internal energy power W consumed by the rock mass on the failure function curve is... in The derivation process is as follows: Assuming the rock mass satisfies the associated Hoek-Brown yield criterion, its plastic potential function is equal to its yield function:

[0036] Based on the flow law and the geometric relationship of the function curve, the normal strain rate and shear strain rate are obtained respectively:

[0037]

[0038]

[0039]

[0040] Combining the above formulas (2) to (5), we obtain the expression for normal stress:

[0041] The energy dissipation power per unit volume of rock mass within the failure function curve f(x) is:

[0042] Thus, the total internal energy power W consumed by the rock mass on the failure function curve is obtained. in :

[0043] In the above formulas (1) to (8), τ n The shear stress at the failure surface is expressed in kPa; σ n The stress is the normal stress on the failure surface, expressed in kPa; σ ci σ represents the uniaxial compressive strength of intact rock, expressed in kPa. tm Here, λ represents the tensile stress parameter of the rock mass, in kPa; A and B are empirical constants; λ is the flow law constant; and u represents the velocity field, in m / s. 2 w represents the thickness of the micro-element in meters; R represents the radius of the cavern in meters; and L0 represents the horizontal distance from the failure function to the surface in meters.

[0044] In step S11 above, the power W generated by gravity is derived. g The expression is:

[0045] Where V is the area of ​​the envelope within the destruction function curve.

[0046] In step S12 above, the derived limiting internal pressure p u The expression is:

[0047] In step S13 above, the equation of the broken function curve f(x, a0, b0, A, B, γ, σ) is broken. ci The specific derivation process is as follows: Define a generic function:

[0048] The Lagrange function can then be expressed as:

[0049] According to the Lagrange equation, we can obtain:

[0050] Integrating gives f'(x), and integrating again gives f(x): .

[0051] f(x, a0, b0, A, B, , ci This is equivalent to the above formula (14), i.e., f(x, a0, b0, A, B, ... , ci )=f(x).

[0052] refer to Figure 4 In this embodiment, step S2 above, the step of calculating the burial depth H1 of a single tunnel using the single-tunnel upper limit theorem limit uplift resistance model, includes: S21: Determine the initial burial depth ; S22: Determine the destruction function constant b0, b0 = R + H0; S23: Set the boundary conditions f(x, a0, b0, A, B, ... , ci )=0, x=R, and the failure function constant b0, rock unit weight Uniaxial compressive strength of intact rock ciSubstituting the empirical constants A and B into the equation of the failure function curve f(x, a0, b0, A, B), , ci Solve for the failure function constant a0; S24: Substitute the destruction function constants a0 and b0 into the destruction function curve equation f(x, a0, b0, A, B, ... , ci The destruction function f(x) is obtained; the horizontal distance L0 from which the destruction function extends to the Earth's surface is calculated, where L0 = a0 / ; S25: Solve for the volume V0 of the damaged breadth soil: ; S26: Solving for the ultimate internal pressure P u : ; S27: Compare the ultimate internal pressure P u The system determines whether the accuracy of the design internal pressure P meets the requirements. If it does, the initial burial depth H0 is set as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S22~S27 are repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H1.

[0053] Optionally, in step S43 above, determining whether the accuracy meets the requirements includes: like Then the ultimate internal pressure P u The internal pressure P of the design meets the accuracy requirements; if Then H0 = H0 + adjustment value.

[0054] In step S27 above, determining whether the accuracy meets the requirements includes: S271: If If the accuracy meets the requirements, then the accuracy is satisfactory; if Then proceed with S272; S272: Compare the ultimate internal pressure P u and the design internal pressure P; S273: If P > P u If H0 = H0 + adjustment value; if P < P u Then H0 = H0 - adjustment value; Wherein, the allowable error = design internal pressure P × error ratio, and the error ratio ≤ 5%; the adjustment value = initial burial depth H0 / 5 × error ratio.

[0055] In this embodiment, the step S42 above, which involves solving for the volume V of the overlying rock mass in the fault optimization model, includes: C421: The volume V0 of the damaged envelope soil is calculated by following the steps from S22 to S25 based on the initial burial depth H0 in S41. C422: The volume of the overlying rock mass in the fault optimization model is V = the volume of the damaged breadth soil V0 - the volume of the fault intersection region V1, i.e.: .

[0056] In this embodiment, the step of solving for the volume V1 of the fault intersection region in step S42 above includes: S421: Establish a coordinate system with the center of the gas storage dome as the origin. First, determine the coordinates of the fault initiation point A1 as (x1, y1), the coordinates of the intersection of the fault and the ground B1 as (x2, y2), the coordinates of the initiation point of the overlying rock fracture surface A2 as (x3, y3), and the coordinates of the intersection of the fracture surface and the ground B2 as (x4, y4). Then, calculate the coordinates of the intersection point c1 of the fault and the fracture surface as (x0, y0). S422: Calculate the height h1 of the intersection region, the width b1 of the cone base of the intersection region, and the chord length l1 of the cone base of the intersection region respectively: h1 = y2 - y0; b1 = x4 - x2; ; S423: Calculate the cone base area S1 of the intersection region: ; S424: Calculate the volume V1 of the intersection region of the fault. .

[0057] Alternatively, in step S421 above, the calculation formulas for the coordinates (x1, y1) of the fault initiation point A1, the coordinates (x2, y2) of the fault-ground intersection point B1, the coordinates (x3, y3) of the initiation point A2 of the overlying rock mass fracture surface, and the coordinates (x4, y4) of the fracture surface-ground intersection point B2 include: x1=L, y1=B, x2=(H1+R-B) / tan +L, y2=y4=H1+R, x3=R, y3=0, x4=L0.

[0058] In this embodiment, the calculation of the coordinates (x0, y0) of the intersection point c1 of the fault and the rupture surface in step S421 includes: S4211: Constructing the fault line equation: ; S4212: Combining the fault line equation and the failure function curve equation f(x, a0, b0, A, B, ... , ci x0 and y0 are obtained by solving using software; where, ; ; .

[0059] In this embodiment, step S4 above, the step of determining whether the fault intersects with the overlying rock mass of the gas storage tank, includes: C41: Projecting the three-dimensional upper bound theorem-limited uplift fault optimization model (e.g.) Figure 5 As shown) to form a two-dimensional fault optimization model (such as Figure 7 As shown), the two-dimensional fault optimization model is aligned with the axis of the gas storage tank 1 and perpendicular to the fault 3. C42: with Figure 7 Based on the two-dimensional fault optimization model shown, the equation of the failure function curve f(x, a0, b0, A, B, ...) is... , ci Differentiating f'(x, a0, b0, A, B) yields f'(x, a0, b0, A, B). , ci Substitute the interval of independent variables (0, L0) into f'(x, a0, b0, A, B, ...); , ci ), thus obtaining the maximum slope k of the failure surface. r k r =max(f'(x,a0,b0,A,B, , ci The maximum dip angle of the failure surface was calculated. ', =arctan(k) r ); C43: When the fault dip angle ≥ When there is no tangent point between the fault and the overlying rock mass, the critical intersection point M(L0, H+R) between the fault and the overlying rock mass is defined; when the fault dip angle is... < At that time, there is a tangent point between the fault and the failure surface of the overlying rock mass of the gas storage reservoir. The tangent point N(x) is defined as follows: c ,y c ), where x c =solve(f'(x c )=tan y c =f(x c ); C44: Determine the critical horizontal distance L r When the fault dip angle ≥ 'hour, When the fault dip angle < 'hour, ; C45: Compare with L r and L, if L r If L ≥ L, then the fault and the overlying rock mass of the gas storage are mutually opposed; if L r If the fault is less than L, then the fault intersects with the overlying rock mass of the gas storage facility.

[0060] The following is a calculation example using the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model of this embodiment: Calculation Example 1: The first phase of a domestic compressed air energy storage peak-shaving power station project has a total installed capacity of 700MW / 3500MWh. The project adopts advanced compressed air energy storage technology, and the gas storage facility is a newly excavated underground artificial cavern with a designed volume of 875,000 m³. 3 The design tunnel has a diameter of 43m and a maximum gas storage pressure of 14.7MPa. The layout of the underground gas storage facility in this project is controlled by the F1 fault on the north side. Geophysical data shows that the F1 fault is wide, with a fracture orientation of 350∠67 and an affected width of about 60m.

[0061] 1. Determine design parameters: Inner radius of the tunnel R = 21.5m, unit weight of rock =25kN / m 3 The initial burial depth H0 = 73.5m, the design internal pressure P = 14700kPa, and the uniaxial compressive strength of the intact rock are... ci =45000 kPa, rock tensile strength t =1100kPa, Hoek-Brown empirical constant A=0.39, Hoek-Brown empirical constant B=0.82, safety factor F s =1.5, Fault F1 dip angle =113°, the horizontal distance between the starting point of fault F1 and the center of the tank is L=65m, and the vertical distance between the starting point of fault F1 and the center of the tank dome is D=0m.

[0062] 2. Substitute the parameters into the single-tunnel upper limit theorem limit uplift resistance model to calculate the single-tunnel burial depth: the burial depth H1 is 73.5m.

[0063] 3. Determine the horizontal projection length of the rupture surface when it extends to the ground surface at the single burial depth H1: L0 = 84.414m.

[0064] 4. Substitute the parameters to determine the positional relationship between the fault and the gas storage tank, and select a suitable model: Determine the maximum dip angle of the failure surface ': The maximum slope on the failure surface is kr=max( f (x) = 1.84, maximum tilt angle =arctan(k) r =61.5°.

[0065] Solve for the coordinates of the critical tangent point (intersection point): fault dip angle ≥ At that time, there was no tangent point between the fault and the gas storage uplift failure surface, and the critical intersection point M (84.414, 95) between the fault and the overlying rock mass was defined.

[0066] Determine the critical horizontal distance L r : .

[0067] Determine the locational relationship between the fault and the overlying rock mass of the gas storage facility: due to Since the fault and the overlying rock mass of the gas storage are intersecting, the upper limit theorem limit uplift fault optimization model should be selected for calculation.

[0068] 5. Determine the initial burial depth H0 of the fault optimization model: The calculation result of the single-cavity upper limit theorem limit anti-uplift model of the large tank gas storage is used as the initial burial depth of the fault optimization model and iteratively set H0=73.5m.

[0069] 6. Determine the destruction function constant b0: .

[0070] 7. Solve for the failure function constant a0: Apply boundary conditions Substituting into equation (14), i.e. f(x, a0, b0, A, B, , ci From this, we can obtain: Solving for a, we get a0 = 2110.3.

[0071] 8. Determine the destructive function f(x): Substitute the constants a0 and b0 obtained above into equation (14), i.e., f(x, a0, b0, A, B, ... , ci The destruction function is obtained from ) .

[0072] 9. Solve for the horizontal distance L0 from the point where the failure function extends to the Earth's surface: .

[0073] 10. Solve for the volume V1 of the fault intersection region: Determine the coordinate parameters: x1=L=65, y1=B=0, x2=(H1+RB) / tan +L=24.67, y2=y4=H1+R=95, x3=R=21.5, y3=0, x4=L0=84.414:

[0074] The ordinate of the intersection point:

[0075] Determine the volume V1 of the overlapping region:

[0076] 11. Solve for the volume V of the overlying rock mass in the fault optimization model: .

[0077] 12. Solve for the ultimate internal pressure p u : .

[0078] 13. Accuracy check: Compare the solution with the ultimate internal pressure p obtained in the previous step. u Is the difference between the design internal pressure P and the design internal pressure within the required accuracy range? (In this case, the error ratio is required to be controlled at 5%, so the allowable error = 14700 × 5% = 735 kPa). The accuracy requirement is not met because P is greater than p. u Then H0 = H0 + adjustment value = H0 + 0.735 = 74.235m, and return to H0 = 74.235m in step 5 to recalculate.

[0079] ... (iterative loop) 14. Determine the initial burial depth H0 of the fault optimization model: H0=H1=83.79m.

[0080] 15. Determine the destruction function constant b0: .

[0081] 16. Solve for the failure function constant a0: Apply boundary conditions Substituting into equation (14), i.e. f(x, a0, b0, A, B, , ci From this, we can obtain: Solving for a, we get a0 = 2248.7.

[0082] 17. Determine the destructive function f(x): Substitute the constants a0 and b0 obtained above into equation (14), i.e., f(x, a0, b0, A, B, ... , ci The destruction function is obtained from ) .

[0083] 18. Solve for the horizontal distance L0 from the point where the failure function extends to the Earth's surface: .

[0084] 19. Solve for the volume V1 of the fault intersection region: Determine the coordinate parameters: x1=L=65, y1=B=0, x2=(H1+RD) / tanb+L=20.31, y2=y4=H1+R=105.29, x3=R=21.5, y3=0, x4=L0=89.95.

[0085] The ordinate of the intersection point:

[0086] Determine the volume V1 of the overlapping region:

[0087] 20. Solve for the volume V of the overlying rock mass in the fault optimization model: .

[0088] 21. Solve for the ultimate internal pressure p u : .

[0089] 22. Accuracy check: Compare the solution with the ultimate internal pressure p obtained in the previous step. u Is the difference between the design internal pressure P and the design internal pressure within the required accuracy range? (In this case, the error ratio is required to be controlled at 5%, so the allowable error = 14700 × 5% = 735 kPa). This meets the accuracy requirements.

[0090] 23. Output safe burial depth H=83.79m.

[0091] Calculation Example 2: The first phase of a domestic compressed air energy storage peak-shaving power station project has a total installed capacity of 700MW / 3500MWh. The project adopts advanced compressed air energy storage technology, and the gas storage facility is a newly excavated underground artificial cavern with a designed volume of 875,000 m³. 3 The design tunnel has a diameter of 43m and a maximum gas storage pressure of 14.7MPa. The layout of the underground gas storage facility in this project is controlled by the F1 fault on the north side. Geophysical data shows that the F1 fault is wide, with a fracture orientation of 350∠67 and an affected width of about 60m.

[0092] 1. Determine design parameters: Inner radius of the tunnel R = 21.5m, unit weight of rock =25kN / m3 The initial burial depth H0 = 73.5m, the design internal pressure P = 14700kPa, and the uniaxial compressive strength of the intact rock are... ci =45000 kPa, rock tensile strength t =1100kPa, Hoek-Brown empirical constant A=0.39, Hoek-Brown empirical constant B=0.82, safety factor F s =1.5, Fault F1 dip angle =113°, the horizontal distance L = 120m between the starting point of fault F1 and the center of the tank, and the vertical distance D = 20m between the starting point of fault F1 and the center of the tank dome.

[0093] 2. Substitute the parameters into the single-tunnel upper limit theorem limit uplift resistance model to calculate the single-tunnel burial depth: the burial depth H1 is 73.5m.

[0094] 3. Determine the horizontal projection length of the rupture surface when it extends to the ground surface at the single burial depth H1: L0 = 84.414m.

[0095] 4. Substitute the parameters to determine the positional relationship between the fault and the gas storage tank, and select a suitable model: Determine the maximum dip angle of the failure surface ': The maximum slope on the failure surface is kr=max( f (x) = 1.84, maximum tilt angle =arctan(k) r =61.5°.

[0096] Solve for the coordinates of the critical tangent point (intersection point): fault dip angle ≥ At that time, there was no tangent point between the fault and the gas storage uplift failure surface, and the critical intersection point M (84.414, 95) between the fault and the overlying rock mass was defined.

[0097] Determine the critical horizontal distance L r : .

[0098] Determine the locational relationship between the fault and the overlying rock mass of the gas storage facility: due to Since the fault and the overlying rock mass of the gas storage are in a divergent relationship (i.e., do not intersect), the upper limit theorem limit uplift resistance model of the large tank gas storage should be selected for calculation.

[0099] 5. Output safe burial depth H=H1=73.5m.

[0100] Example 2 This embodiment provides an electronic device.

[0101] The electronic device of this embodiment includes a memory and at least one processor. The memory stores a computer program, and the processor executes the computer program to implement the burial depth design method based on the upper bound theorem limit uplift fault optimization model as described in Embodiment 1.

[0102] Specifically, the processor may include a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application. The memory may include a mass storage device for data or instructions. For example, and not limitingly, the memory may include a hard disk drive (HDD), a floppy disk drive, a solid-state drive (SSD), flash memory, an optical disk drive, a magneto-optical disk drive, magnetic tape, or a Universal Serial Bus (USB) drive, or a combination of two or more of these. Where appropriate, the memory may include removable or non-removable (or fixed) media. Where appropriate, the memory may be internal or external to the data processing device. In a particular embodiment, the memory is non-volatile. Volatile memory. In a particular embodiment, the memory includes read-only memory. ROM (ROM-only memory) and RAM (Random Access Memory). Where appropriate, the ROM can be a mask-programmed ROM or a programmable ROM. Only Memory (PROM) and Erasable Programmable Read-Only Memory (EPRROM) Only Memory (FPROM) and Electrically Erasable Programmable Read-Only Memory (PROM) Only Memory (EFPROM) and Electrically Alterable Read ROM (EEPROM) Only Memory (EAROM) or Flash Memory, or a combination of two or more of these. Where appropriate, this RAM can be Static Random Access Memory (SRAM). Access Memory (SRAM) or Dynamic Random Access Memory (DRAM) can be either SRAM or DRAM. DRAM can be Fast Page Mode Dynamic Random Access Memory (FPMDRAM), Extended Data Out Dynamic Random Access Memory (EDODRAM), or Synchronous Dynamic Random Access Memory (SDRAM). Access Memory (SDRAM) and similar technologies.

[0103] Memory can be used to store or cache various data files that need to be processed and / or communicated, as well as possible computer program instructions executed by the processor.

[0104] The processor implements any of the methods described in the above embodiments by reading and executing computer program instructions stored in memory. In some embodiments, the electronic device may further include a communication interface and a bus. The processor, memory, and communication interface are connected via the bus and communicate with each other.

[0105] A bus, including hardware, software, or both, couples components of a computer device together. Buses include, but are not limited to, at least one of the following: Data Bus, Address Bus, Control Bus, Expansion Bus, and Local Bus. For example, and not as a limitation, a bus may include Accelerated Graphics Port (AGP) or other graphics buses, Extended Industry Standard Architecture (EISA) buses, Front Side Bus (FSB), Hyper Transport (HT) interconnects, Industry Standard Architecture (ISA) buses, InfiniBand interconnects, Low Pin Count (LPC) buses, memory buses, Micro Channel Architecture (MCA) buses, Peripheral Component Interconnect (PCI) buses, and PCI... Express (PCI The bus may be an X-bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local Bus (VLB) bus, or other suitable buses, or a combination of two or more of these. Where appropriate, the bus may include one or more buses. Although specific buses are described and illustrated in embodiments of this application, this application contemplates any suitable bus or interconnect.

[0106] Example 3 This embodiment provides a computer-readable storage medium.

[0107] The computer-readable storage medium of this embodiment stores a computer program, which, when executed, implements the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model as described in Embodiment 1.

[0108] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium, and when executed, it performs the steps of the above method embodiments. When the integrated unit of the present invention is implemented as a software functional unit and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiments of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the method of Embodiment 1 of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, ROMs, magnetic disks, or optical disks.

[0109] In summary, this invention proposes a burial depth design method, equipment, and medium based on an upper-boundary fault optimization model. This model is a three-dimensional rotating body model, which better conforms to the stress mode of large tank-type gas storage facilities. The spatial characteristics of this model make it more sensitive to adverse geological conditions such as faults, thus requiring consideration of the weakening effect of faults on the surrounding rock within the failure surface. This method determines the critical tangent point (or intersection point) by solving for the maximum slope of the failure surface, and introduces a critical horizontal distance to define the abscissa of the fault initiation point when the overlying rock mass and the fault are just tangent (intersecting). Based on this parameter, the positional relationship between the fault and the overlying rock mass of the gas storage facility can be quickly determined, thus serving as a parameter for avoidance conditions during the early site selection of the gas storage facility. At the same time, by comparing the actual fault initiation coordinates with the critical horizontal distance, the size of the intersection area can be calculated, thereby further evaluating the adverse effects of the fault on the gas storage facility. When a fault intersects with the surrounding rock within the overlying failure surface of a gas storage facility, the fault optimization calculation model for large tank gas storage facilities recommended in this method should be adopted. By correcting the fracture angle of the surrounding rock in the fault intersection area to the fault dip angle, a more realistic stress and failure mode can be restored, making the calculation results more realistic and reliable. Taking the calculation example 1 in the embodiment as an example, after considering the fault conditions, the calculated burial depth increased from 73.5 to 83.79, an increase of 14%, and the correction result is significant.

[0110] Compared to the single-cavity model, the fault optimization model takes into account the influence of faults. Therefore, it is necessary to determine the positional relationship between the fault and the gas storage uplift body, calculate the volume of the intersection zone, and use the fault optimization model calculation process for calculation.

[0111] The main differences between this invention and the three-dimensional rigid frustum fault optimization model are as follows: 1. The failure surfaces are different. The failure surface established in this invention considers more parameters than the three-dimensional rigid frustum model. The failure surface definition is more accurate and has a wider range of applications. The failure angle of 30-45° of the three-dimensional rigid frustum is suitable for harder rock masses. The failure surface definition of the three-dimensional rigid frustum model is relatively simple, and the model parameters are easy to obtain. Therefore, this method can help with calculation and evaluation when specific geotechnical parameters cannot be obtained in the early stage of the project.

[0112] 2. The two limit criteria are different: Due to the difference in model definition, the three-dimensional rigid frustum model uses the limit dip angle to define the critical state and determines the positional relationship by comparing the actual fault dip angle with the critical dip angle; the upper limit theorem limit uplift fault optimization model of this invention defines the critical state by the critical horizontal distance and determines the positional relationship by comparing the horizontal coordinate of the fault starting point with the critical horizontal distance, which makes the evaluation of fault influence more intuitive. At the same time, the critical horizontal distance Lr in this method can also be used as a basis for not being disturbed by faults in the early site selection.

[0113] 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 within the protection scope of the present invention.

Claims

1. A method for designing burial depth based on an optimization model for limit uplift faults using an upper bound theorem, characterized in that, include: S1: Determine the design parameters of the gas storage facility, including the inner radius R of the top sphere of the gas storage facility and the unit weight of the overlying rock mass of the gas storage facility. The design internal pressure P, initial burial depth H0, and uniaxial compressive strength of the intact rock in the gas storage facility are specified. ci Empirical constant A, empirical constant B, safety factor F s Fault dip angle The horizontal distance L between the fault initiation point and the center of the gas storage facility; the vertical distance D between the fault initiation point and the center of the gas storage facility dome; construct the failure function curve equation f(x, a0, b0, A, B, ... , ci ); S2: Substitute the design parameters into the single-tunnel upper limit theorem limit uplift resistance model to calculate the burial depth H1; S3: Determine the horizontal projection length L0 of the rupture surface corresponding to the burial depth H1 when it extends to the ground surface; S4: Determine whether the fault intersects with the overlying rock mass of the gas storage. If not, take the burial depth H1 in S2 as the safe burial depth H; if so, use the upper limit theorem limit uplift fault optimization model to calculate the burial depth H2. The calculation steps for calculating the burial depth H2 using the upper bound theorem-based limit uplift fault optimization model include: S41: Take the burial depth H1 in S2 as the initial burial depth H0; S42: Solve for the volume V1 of the intersection region of the fault, calculate the volume V of the overlying rock mass in the fault optimization model based on V1, and then establish an equilibrium relationship based on V to solve for the ultimate internal pressure P. u : ; S43: Compare the ultimate internal pressure P u The accuracy of the design internal pressure P is judged to see if it meets the requirements. If it does, the initial burial depth H0 is used as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S42~S43 is repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H2.

2. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 1, characterized in that, In S2, the steps for calculating the burial depth H1 of the single-hole tunnel using the single-hole upper limit theorem limit uplift resistance model include: S21: Determine the initial burial depth ; S22: Determine the destruction function constant b0, b0 = R + H0; S23: Set the boundary conditions f(x, a0, b0, A, B, ... , ci )=0, x=R, and the failure function constant b0, the unit weight of the rock The uniaxial compressive strength of the intact rock ci Substituting the empirical constants A and B into the equation f(x, a0, b0, A, B), , ci Solve for the failure function constant a0; S24: Substitute the destruction function constants a0 and b0 into the destruction function curve equation f(x, a0, b0, A, B, ... , ci The destruction function f(x) is obtained; the horizontal distance L0 from which the destruction function extends to the Earth's surface is calculated, where L0 = a0 / ; S25: Solve for the volume V0 of the damaged breadth soil: ; S26: Solving for the ultimate internal pressure P u : ; S27: Compare the ultimate internal pressure P u The accuracy of the design internal pressure P is judged to see if it meets the requirements. If it does, the initial burial depth H0 is used as the safe burial depth H and output. If it does not meet the requirements, the initial burial depth H0 is adjusted, and S22~S27 is repeated until the accuracy meets the requirements. The safe burial depth H is then output: H=H1.

3. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 2, characterized in that, In S43, determining whether the accuracy meets the requirements includes: like Then the ultimate internal pressure P u The designed internal pressure P meets the accuracy requirements; if Then H0 = H0 + adjustment value; In S27, determining whether the accuracy meets the requirements includes: S271: If If the required accuracy is met; if Then proceed with S272; S272: Compare the ultimate internal pressure P u and the designed internal pressure P; S273: If P > P u If H0 = H0 + adjustment value; if P < P u Then H0 = H0 - adjustment value; Wherein, the allowable error = the design internal pressure P × error ratio, and the error ratio ≤ 5%; the adjustment value = the initial burial depth H0 / 5 × the error ratio.

4. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 2, characterized in that, In S42, the steps for solving the volume V of the overlying rock mass in the fault optimization model include: C421: The volume V0 of the damaged envelope soil is calculated by following the steps from S22 to S25 based on the initial burial depth H0 described in S41. C422: The volume V of the overlying rock mass in the fault optimization model = the volume V0 of the damaged envelope soil mass - the volume V1 of the intersection region of the fault, that is: 。 5. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 1, characterized in that, In S42, the steps for solving the volume V1 of the intersection region of the fault include: S421: Establish a coordinate system with the center of the gas storage dome as the origin. First, determine the coordinates of the fault starting point A1 as (x1, y1), the coordinates of the intersection of the fault and the ground B1 as (x2, y2), the coordinates of the starting point of the overlying rock fracture surface A2 as (x3, y3), and the coordinates of the intersection of the fracture surface and the ground B2 as (x4, y4); then calculate the coordinates of the intersection point c1 of the fault and the fracture surface as (x0, y0). S422: Calculate the height h1 of the intersection region, the width b1 of the cone base of the intersection region, and the chord length l1 of the cone base of the intersection region, respectively. h1 = y2 - y0; b1 = x4 - x2; ； S423: Calculate the cone-shaped base area S1 of the intersection region: ； S424: Calculate the volume V1 of the intersection region of the fault: 。 6. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 5, characterized in that, In S421, the calculation formulas for the coordinates (x1, y1) of the fault initiation point A1, the coordinates (x2, y2) of the fault-ground intersection point B1, the coordinates (x3, y3) of the initiation point A2 of the overlying rock mass fracture surface, and the coordinates (x4, y4) of the fracture surface-ground intersection point B2 are as follows: include: x1=L,y1=B,x2=(H1+R-B) / tan +L,y2=y4=H1+R,x3=R,y3=0,x4=L0。 7. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 5, characterized in that, In S421, the calculation steps for the coordinates (x0, y0) of the intersection point c1 of the fault and the rupture surface include: S4211: Constructing the fault line equation: ; S4212: Combining the fault line equation and the failure function curve equation f(x, a0, b0, A, B, ... , ci x0 and y0 are obtained by solving using software; where, ； ； 。 8. The burial depth design method based on the upper bound theorem-based limit uplift fault optimization model according to claim 1, characterized in that, In S4, the step of determining whether the fault intersects with the overlying rock mass of the gas storage tank includes: C41: Project the three-dimensional upper bound theorem limit uplift fault optimization model to form a two-dimensional fault optimization model, such that the two-dimensional fault optimization model is along the axis of the gas storage tank and perpendicular to the fault. C42: Based on the aforementioned two-dimensional fault optimization model, the equation of the failure function curve f(x, a0, b0, A, B, ...) is... , ci Differentiating f'(x, a0, b0, A, B) yields f'(x, a0, b0, A, B). , ci Substitute the interval of independent variables (0, L0) into f'(x, a0, b0, A, B, ...); , ci ), thus obtaining the maximum slope k of the failure surface. r k r =max(f'(x,a0,b0,A,B, , ci The maximum dip angle of the failure surface was calculated. ', =arctan(k) r ); C43: When the fault dip angle ≥ When there is no tangent point between the fault and the overlying rock mass, the critical intersection point M(L0, H+R) between the fault and the overlying rock mass is defined; when the fault dip angle is... < At that time, there is a tangent point between the fault and the failure surface of the overlying rock mass of the gas storage reservoir. The tangent point N(x) is defined as follows: c ,y c ), where x c =solve(f'(x c )=tan y c =f(x c ); C44: Determine the critical horizontal distance L r When the fault dip angle ≥ 'hour, When the fault dip angle < 'hour, ; C45: Compare with L r and L, if L r If L ≥ L, then the fault and the overlying rock mass of the gas storage are mutually opposed; if L r If the fault is less than L, then the fault intersects with the overlying rock mass of the gas storage facility.

9. An electronic device, characterized in that, The electronic device includes a memory and at least one processor, the memory storing a computer program, and the processor executing the computer program to implement the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model as described in any one of claims 1 to 8.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, which, when executed, implements the burial depth design method based on the upper bound theorem-based limit uplift fault optimization model as described in any one of claims 1 to 8.

Images