Underground water solid tide response fracture aquifer parameter inversion method and system

Abstract

The invention relates to an underground water solid tide response fissure aquifer parameter inversion method and system, and the method comprises the steps: building a fissure aquifer underground water solid tide response mathematical model considering pore-fissure dual media, and solving the mathematical model to obtain an analytical solution; hydrological borehole water level data of the to-be-estimated fissure aquifer are obtained, and the amplitude ratio and phase difference characteristic parameters of the water level responding to the change of the earth tide stress are obtained; fitting wiring is carried out on the analytic solution, the amplitude ratio and the phase difference characteristic parameters, and a parameter value corresponding to a fitting curve with the maximum determination coefficient R2 is selected as an inversion result of the fractured aquifer parameters; the dual-medium model can truly describe the movement law of underground water in a pore-fracture system in a fractured aquifer under the earth tide response, and can visually understand the influence of hydrogeological parameters of dual media in the fractured aquifer on the earth tide response of underground water. And a more reliable theoretical basis and method are provided for accurate inversion of fractured aquifer parameters.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of groundwater parameter inversion, more particularly to a groundwater solid tide response fractured aquifer parameter inversion method and system. BACKGROUND

[0002] As a key freshwater resource, groundwater plays a vital role in global water supply and ecosystem maintenance. However, with the increasing uneven distribution of precipitation caused by climate change and the intensifying water resource shortage, the rational development and sustainable utilization of groundwater resources have become a key issue that needs to be addressed by countries. In this context, the determination of aquifer parameters has become a core factor in addressing these challenges. As a natural hydrodynamic phenomenon, the solid tide provides an important basis for aquifer parameter inversion. Accurate aquifer parameter inversion results not only provide reliable basic data for numerical simulation and water flow prediction of aquifers, but also are the key to effective management and protection of water resources.

[0003] Fractured aquifers, as geological media for groundwater storage and movement, have complex internal structure and properties that significantly affect groundwater flow patterns. In real geological environments, they exhibit obvious dual-medium characteristics, consisting of both pore and fracture media. Traditional equivalent porous media have obvious shortcomings in describing this heterogeneity, typically assuming that aquifers are homogeneous and isotropic continuous media, ignoring the coupling between pore and fracture media, resulting in low inversion accuracy. Building a groundwater solid tide response fractured aquifer parameter inversion method and corresponding system based on a dual-medium model is an economical and effective means. Combined with long-term water level observation data, it can more accurately invert the true physical properties of aquifers, providing more accurate and efficient data support for related work. SUMMARY

[0004] The technical problem to be solved by the present application is to provide a groundwater solid tide response fractured aquifer parameter inversion method and a groundwater solid tide response fractured aquifer parameter inversion system to address the above-mentioned defects of the prior art.

[0005] The technical solution adopted by the present application to solve the technical problem is:

[0006] A groundwater solid tide response fractured aquifer parameter inversion method is constructed, wherein the method comprises the following steps:

[0007] A dual-medium mathematical model considering the groundwater response to solid tide in the pore-fracture system of fractured aquifers is established, and the mathematical model based on dual-medium theory is solved to obtain the analytical solution of the water head amplitude ratio and phase difference of fractured aquifers under tidal stress;

[0008] Obtain the hydrological borehole water level data of the aquifer to be estimated, and obtain the amplitude ratio and phase difference characteristic parameters of the water level responding to the solid tide stress change;

[0009] The analytical solution is fitted with the obtained amplitude ratio and phase difference characteristic parameters of the water level responding to the solid tide stress change, and the determination coefficient R 2 The parameter value corresponding to the maximum fitting curve is taken as the inversion result of the fractured aquifer parameter.

[0010] The underground water solid tide response fractured aquifer parameter inversion method provided by the application, wherein the double medium mathematical model considering the underground water responding to the solid tide in the pore-fracture system of the fractured aquifer is established based on the conditions that

[0011] The confined aquifer is considered to be horizontally infinite, uniform in thickness, and the vertical leakage recharge is ignored;

[0012] The aquifer pores and fractures are uniformly distributed in space;

[0013] Each fracture and matrix continuum is assumed to be homogeneous and isotropic;

[0014] At any point, there are double water heads of pores and fractures, and there is water exchange between the two, and the exchange amount depends on the water head difference between the media and obeys Darcy's law;

[0015] The pumping well is located in the fractured medium, the well wall uniformly intakes water, and the well is complete in pumping water.

[0016] The underground water solid tide response fractured aquifer parameter inversion method provided by the application, wherein the double medium mathematical model considering the underground water responding to the solid tide in the pore-fracture system of the fractured aquifer is established, the mathematical model based on the double medium theory is solved, and the analytical solution of the amplitude ratio and phase difference of the water head of the fractured aquifer under the action of the tidal stress is obtained, and the analytical solution includes:

[0017] The water flow control equation in the double medium model can be described by the following equation:

[0018]

[0019] In the formula, T is the aquifer transmissivity; h f , h m Water head in fracture, matrix respectively; s f , s m Reservoir rate in fracture, matrix respectively; B is the Skempton coefficient of the aquifer; K u It is the undrained bulk modulus of the aquifer; ρ is the density of water; g is the acceleration of gravity; ε is the pore-elastic volume strain generated by the earth tide; r is the radial distance from the wellbore.

[0020] Fluid exchange between the pore and the fracture is given by the equation:

[0021]

[0022] where μ is the water exchange parameter between the pore and the fracture;

[0023] The boundary conditions are set as:

[0024] h f (r,t) = h ∞ (t), r = 0 (2)

[0025] h m (r,t) = h ∞ (t), r = ∞ (3)

[0026] h f (r,t) = h w (t), r = r w (4)

[0027] h m (r,t) = h w (t), r = r w (5)

[0028]

[0029] where r w and r c are the outer and inner diameters of the wellbore, respectively;

[0030] The seepage continuity equation of formula (1) (2) at infinity is:

[0031]

[0032] Since infinity h ∞ (m) = h ∞,0 e iωt , h ∞ (m) = h ∞,0 e iωt , ε0 is the amplitude of ε, and by substituting the above equation, we get:

[0033]

[0034] where h w (m) = h w,0 e iωt is the periodic water level in the well, with complex amplitude h w,0 ; ω = 2π / τ is the angular frequency, and τ is the tidal oscillation period. It is assumed that:

[0035] h f (r,t) = Δhf (r, t) + h ∞ (t) (11)

[0036] h m (r, t) = Ah m (r, t) + h ∞ (t) (12)

[0037] Substituting equation (11) into equation (1) gives:

[0038]

[0039] According to equation (8), we have:

[0040]

[0041] Substituting equation (12) into equation (2) gives:

[0042]

[0043] According to equation (9), we have:

[0044]

[0045] Let Ah f = Ah f,0 (r) e iωt , Ah m = Ah m,0 (r) e iωt , equations (15) and (17) are expressed as:

[0046]

[0047] iωs m Ah m,0 = μ(Ah f,0 - Ah m,0 ) (19)

[0048] Then,

[0049]

[0050] Substituting equation (20) into equation (18) gives:

[0051]

[0052] At this time, the boundary condition is:

[0053] Ah f,0 (r→∞) = 0 (22)

[0054]

[0055] The general solution of formula (21) is:

[0056] Δh f,0 =C I I0(βr)+C K K0(βr) (25)

[0057] In the formula, I0 and K0 are the first-order and second-order modified Bessel functions of zero order respectively;

[0058] Then,

[0059]

[0060] When r approaches infinity, the first-order modified Bessel function of zero order I0 increases exponentially, and according to the boundary condition (18), C I = 0, then:

[0061] Δh f,0 =C K K0(βr) (27)

[0062] According to formula (24), C K Expression:

[0063]

[0064] Bring formula (29) into formula (23):

[0065]

[0066] Wherein,

[0067]

[0068] Define the amplitude ratio:

[0069]

[0070] The phase difference is:

[0071]

[0072] The groundwater solid tide response fissure aquifer parameter inversion method, wherein the hydrological borehole water level data of the fissure aquifer to be estimated is obtained. Through the baytap solid tide analysis program, the amplitude ratio and the phase difference of the water head of the fissure aquifer under the action of the solid tide stress are obtained.

[0073] A groundwater solid tide response fissure aquifer parameter inversion system is applied to the groundwater solid tide response fissure aquifer parameter inversion method, wherein the system comprises a model construction unit, a data acquisition and processing unit and a fitting wiring unit.

[0074] The model construction unit is configured to establish a double medium mathematical model of groundwater response to earth tide in a fissure aquifer pore-fracture system, solve the mathematical model based on the double medium theory, and obtain an analytical solution of a water head amplitude ratio and a phase difference of the fissure aquifer under the action of tidal stress;

[0075] The data acquisition and processing unit is configured to obtain hydrological borehole water level data of the fissure aquifer to be estimated, and obtain amplitude ratio and phase difference characteristic parameters of the water level response to the earth tide stress change;

[0076] The fitting and wiring unit is configured to fit and wire the analytical solution with the obtained amplitude ratio and phase difference characteristic parameters of the water level response to the earth tide stress change, select a determination coefficient R 2 The parameter value corresponding to the maximum fitting curve is taken as the inversion result of the fissure aquifer parameters.

[0077] The application has the beneficial effects that by using the double medium model, the application can more realistically depict the groundwater movement law in the pore-fracture system of the fissure aquifer under the response of the earth tide, can intuitively understand the influence of the hydrogeological parameters of the double medium in the fissure aquifer on the groundwater earth tide response, and can provide more reliable theoretical basis and method for the accurate inversion of the fissure aquifer parameters. BRIEF DESCRIPTION OF DRAWINGS

[0078] In order to more clearly illustrate the technical solutions in the embodiments of the application or the prior art, the application will be further described below with reference to the drawings and embodiments. The drawings in the following description are only some embodiments of the application, and those skilled in the art can obtain other drawings according to these drawings without creative labor:

[0079] Figure 1 is a groundwater earth tide response fissure aquifer parameter inversion method flow chart of a preferred embodiment of the application;

[0080] Figure 2 is a groundwater earth tide response fissure aquifer parameter inversion method earth tide schematic diagram of a preferred embodiment of the application;

[0081] Figure 3 is a groundwater earth tide response fissure aquifer parameter inversion method phase difference fitting diagram of a preferred embodiment of the application;

[0082] Figure 4 is a principle block diagram of a groundwater earth tide response fissure aquifer parameter inversion system of a preferred embodiment of the application. DETAILED DESCRIPTION

[0083] In order to make the purposes, technical solutions and advantages of the embodiments of the present application clearer, the following will make a clear and complete description of the technical solutions in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort belong to the protection scope of the present application.

[0084] The groundwater solid tide response fractured aquifer parameter inversion method of the preferred embodiment of the present application, as shown in Figure 1 , and referring to Figure 2 and Figure 3 , comprises the steps of:

[0085] S01: establishing a double medium mathematical model considering the groundwater response to the solid tide in the pore-fracture system of the fractured aquifer, solving the mathematical model based on the double medium theory, and obtaining the analytical solution of the amplitude ratio and phase difference of the water head of the fractured aquifer under the tidal stress;

[0086] S03: fitting and wiring the analytical solution and the amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress change, selecting the fitting curve corresponding to the maximum determination coefficient R 2 , and taking the parameter value as the inversion result of the aquifer parameter;

[0087] By using the double medium model, the groundwater movement law in the pore-fracture system of the fractured aquifer under the solid tide response can be truly depicted, and the influence of the hydrogeological parameters of the double medium in the fractured aquifer on the groundwater solid tide response can be intuitively understood, thereby providing a more reliable theoretical basis and method for the accurate inversion of the fractured aquifer parameters.

[0088] In the application aspect, the inversion method based on the double medium model can significantly improve the accuracy of parameter inversion. By combining with the solid tide signal, more effective information can be extracted from the hydrological monitoring data, and after processing by the optimization algorithm, the aquifer parameter estimation value closer to the actual situation is obtained. This has important practical significance for realizing the fine numerical simulation and prediction of the aquifer, optimizing the development and utilization scheme of the groundwater resources, strengthening the prevention and control of the groundwater pollution, and ensuring the sustainable utilization of the groundwater resources.

[0089] The double medium mathematical model considering the groundwater response to the solid tide in the pore-fracture system of the fractured aquifer is established based on the condition that:

[0090] The confined aquifer is considered to be horizontally infinite, uniform in thickness, and the vertical leakage recharge is ignored;

[0091] The pore and fracture of the aquifer are evenly distributed in space;

[0092] Each fracture and matrix continuum is assumed to be homogeneous and isotropic;

[0093] At any point, there is a double water head of pore and fracture, and there is a water exchange between the two, and the exchange amount depends on the water head difference between the media and obeys Darcy's law;

[0094] The pumping well is located in the fractured medium, and the well wall is uniformly watered, and the water is pumped out of the complete well. Figure 2 , a double medium mathematical model considering the response of groundwater in the pore-fracture system of the fractured aquifer to the earth tide is established, the mathematical model based on the double medium theory is solved, and the analytical solution of the amplitude ratio and phase difference of the water head of the fractured aquifer under the action of tidal stress is obtained. More specific operations include:

[0095]

[0096] In the formula, T is the aquifer transmissibility; h f , h m are the water heads in the fracture and matrix respectively; s f , s m are the water storage rates in the fracture and matrix respectively; B is the Skempton coefficient of the aquifer; K u is the undrained bulk modulus of the aquifer; ρ is the density of water; g is the acceleration of gravity; ε is the pore-elastic volume strain generated by the earth tide; r is the radial distance from the wellbore.

[0097] The fluid exchange between the pore and the fracture is calculated by the equation:

[0098]

[0099] Where μ is the water exchange parameter between the pore and the fracture;

[0100] The boundary conditions are set as:

[0101] h f (r,t)=h ∞ (t),r=∞ (3)

[0102] h m (r,t)=h ∞ (t),r=∞ (4)

[0103] h f (r,t)=h w (t),r=r w (5)

[0104] h m (r,t)=hw (t), r = r w (6)

[0105]

[0106] where r w , r c are the outer and inner diameters of the wellbore, respectively.

[0107] The flow continuity equation of formula (1) (2) at infinity is:

[0108]

[0109] Since the infinity h ∞ (m) = h ∞,0 e iωt , h ∞ (m) = h ∞,0 e iωt , ε0 is the amplitude of ε, and the above equation can be obtained by substituting it into the equation:

[0110]

[0111] where h w (m) = h w,0 e iωt is the periodic water level in the well, with complex amplitude h w,0 ; ω = 2π / τ is the angular frequency, τ is the tidal oscillation period, and it is assumed that:

[0112] h f (r, t) = Δh f (r, t) + h ∞ (t) (11)

[0113] h m (r, t) = Δh m (r, t) + h ∞ (t) (12)

[0114] Substituting formula (11) into formula (1) can be obtained:

[0115]

[0116] According to formula (8), we can get:

[0117]

[0118] Substituting formula (12) into formula (2) can be obtained:

[0119]

[0120] According to formula (9), we can get:

[0121]

[0122] Let Δh f = Δh f,0 (r) e iωt , Δh m = Δh m,0 (r) e iωt , formula (15) (17) is expressed as:

[0123]

[0124] iωs m Δh m,0 = μ(Δh f,0 - Δh m,0 ) (19)

[0125] Then

[0126]

[0127] Substitute formula (20) into formula (18) to obtain:

[0128]

[0129] At this time, the boundary condition is:

[0130] Δh f,0 (r→∞) = 0 (22)

[0131]

[0132] The general solution of formula (21) is:

[0133] Δh f,0 = C I I0(βr) + C K K0(βr) (25)

[0134] In the formula, I0, K0 are the first type and the second type zero-order modified Bessel functions respectively;

[0135] Then,

[0136] When r→∞, the first type zero-order modified Bessel function I0 exponentially increases, and according to the boundary condition (18), C I = 0, then:

[0137] Δh f,0 = C K K0(βr) (27)

[0138] According to formula (24), C K Expression:

[0139]

[0140]

[0141] Substitute equation (29) into equation (23):

[0142]

[0143] where,

[0144]

[0145] Define the amplitude ratio:

[0146]

[0147] The phase difference is:

[0148]

[0149] Adjust the upper and lower limits of the parameters, compare the amplitude ratio and phase difference expressions of the analytical model with the monitoring data, and perform fitting ( Figure 3 ), find the best parameter combination, and determine the estimated values of each permeability coefficient: T = 1 × 10 -8 m 2 / d, S m = 0.016, S f = 0.023, μ = 0.1.

[0150] As a natural and stable periodic stress load acting on the shallow crust for a long time, the solid tide causes the wellbore water level to produce corresponding periodic microdynamic changes. As shown in Figure 2 , the combined effect of the lunar gravity and the inertial centrifugal force generated by the earth's rotation and revolution, the tidal stress causes the elastic deformation of the solid skeleton of the aquifer, causing the groundwater to flow into and out of the wellbore with the compression and expansion of the solid skeleton, resulting in periodic microdynamic changes in the wellbore water level.

[0151] After obtaining the periodically changing water level data, the microdynamic components affected by the solid tide can be extracted from the macro water level dynamics by filtering method according to the known frequency range of the solid tide, and the rest can be removed. Then, the harmonic analysis program baytap for solid tide analysis can be used to analyze the microdynamic components to obtain the corresponding amplitude ratio and phase difference data.

[0152] After the amplitude ratio and phase difference analytical solution code is compiled in matlab, the obtained data is fitted and compared, and finally the fitting result graph shown in Figure 3 is drawn.

[0153] The hydrological borehole water level data of the fractured aquifer to be estimated is preferably acquired. The amplitude ratio and phase difference of the water head of the fractured aquifer under the action of the solid tide stress are obtained through a baytap solid tide analysis program.

[0154] A system for inverting parameters of a groundwater solid tide response fractured aquifer, applied to the method for inverting parameters of a groundwater solid tide response fractured aquifer as described above, as shown in the accompanying drawings, the system comprises a model construction unit 10, a data acquisition and processing unit 11 and a fitting wiring unit 12. Figure 4

[0155] The model construction unit 10 is configured to establish a double medium mathematical model for considering the groundwater response to the solid tide in the pore-fracture system of the fractured aquifer, to solve the mathematical model based on the double medium theory, and to obtain an analytical solution of the amplitude ratio and phase difference of the water head of the fractured aquifer under the action of the tidal stress.

[0156] The data acquisition and processing unit 11 is configured to acquire the hydrological borehole water level data of the fractured aquifer to be estimated, and to obtain the amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress variation.

[0157] The fitting wiring unit 12 is configured to fit and wire the analytical solution with the obtained amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress variation, to select the parameter value corresponding to the fitting curve with the maximum determination coefficient R 2 as the inversion result of the parameters of the fractured aquifer.

[0158] By using the system of the present application, the groundwater movement law in the pore-fracture system of the fractured aquifer under the response of the solid tide can be truly depicted by using the double medium model, and the influence of the hydrogeological parameters of the double medium in the fractured aquifer on the groundwater solid tide response can be intuitively understood, thereby providing a more reliable theoretical basis and method for the accurate inversion of the parameters of the fractured aquifer. It should be understood that, for those skilled in the art, improvements or changes can be made according to the above description, and all such improvements and changes shall fall within the protection scope of the appended claims of the present application.​

Claims

1. A method for inverting parameters of fractured aquifers in response to groundwater solid tide, characterized in that: The method comprises the steps of: A dual-medium mathematical model is established to consider the groundwater response to solid tides in the pore-fracture system of a fractured aquifer. The dual-medium mathematical model is solved to obtain analytical solutions for the hydraulic head amplitude ratio and phase difference of the fractured aquifer under tidal stress. Obtain the hydrological borehole water level data of the fractured aquifer to be estimated, and obtain the amplitude ratio and phase difference characteristic parameters of the water level response to the solid tidal stress change; The analytical solution is fitted with the obtained amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress change, and the determination coefficient R is selected. 2 The parameter value corresponding to the largest fitting curve is taken as the inversion result of the fracture aquifer parameters.

2. The method for inverting parameters of fractured aquifers in response to groundwater solid tide according to claim 1, characterized in that: The dual-medium mathematical model for considering groundwater response to solid tide in the pore-fracture system of a fractured aquifer is established based on the following conditions: The confined aquifer is considered to be infinitely extended horizontally with uniform thickness, and vertical overflow recharge is neglected; The pores and cracks in the aquifer are evenly distributed in space; Each fracture and matrix continuum is assumed to be homogeneous and isotropic; At any point, there is a dual hydraulic head of pores and cracks, and there is water exchange between the two. The exchange amount depends on the hydraulic head difference between the media and obeys Darcy's law. The pumping well is located in a fractured medium, water is evenly introduced into the well wall, and the pumping is done from a complete well.

3. The method for inverting parameters of fractured aquifers in response to groundwater solid tide according to claim 1 or 2, characterized in that: The establishment of a dual-medium mathematical model considering the groundwater response to solid tide in the pore-fracture system of the fractured aquifer, solving the dual-medium mathematical model, and obtaining an analytical solution of the hydraulic head amplitude ratio and phase difference of the fractured aquifer under tidal stress include: The water flow control equation in the dual medium model can be described by the following equation: Where T is the aquifer conductivity; h f 、h m are the hydraulic heads in the cracks and matrix, respectively; s f 、s m are the water storage rates in fractures and matrix respectively; B is the Skempton coefficient of the aquifer; K u is the undrained bulk modulus of the aquifer; ρ is the density of water; g is the acceleration of gravity; ε is the poroelastic volume strain due to the Earth's tides; r is the radial distance from the wellbore; The fluid exchange between pores and fractures is calculated using the equation: Wherein, μ is the water exchange parameter between pores and fractures; Set the boundary conditions: h f (r,t)=h ∞ (t),r=∞ (3) h m (r,t)=h ∞ (t),r=∞ (4) h f (r,t)=h w (t),r=r w (5) h m (r,t)=h w (t),r=r w (6) Where r w 、r c are the outer and inner diameters of the wellbore, respectively; The seepage continuity equation of formula (1) (2) at infinity is: Since the infinite h ∞ (m) = h ∞,0 e iωt , h ∞ (m) = h ∞,0 e iωt , ε0 is the amplitude of ε, and substituting it into the above equation, we can get: Where h w (m) = h w,0 e iωt is the periodic water level in the well, with complex amplitude h w,0 ; ω = 2π / τ is the angular frequency, τ is the tidal oscillation period, assuming: h f (r,t)=Δh f (r,t)+h ∞ (t) (11) h m (r,t)=Δh m (r,t)+h ∞ (t) (12) Substituting formula (11) into formula (1) yields: According to formula (8), we can get: Substituting formula (12) into formula (2) yields: According to formula (9), we can get: Let Δh f =Δh f,0 (r)e iωt , Δh m =Δh m,0 (r)e iωt , formula (15) (17) is expressed as: yes m Dh m,0 =μ(Δh f,0 -Dh m,0 ) (19) but, Substituting formula (20) into formula (18) yields: The boundary conditions at this time are: Δh f,0 (r→∞)=0 (22) The general solution of formula (21) is: Δh f,0 =C I I0(βr)+C K K0(βr) (25) Where I0 and K0 are the first and second kind zero-order modified Bessel functions, respectively; but, When r→∞, the first-order zero-order modified Bessel function I0 grows exponentially. According to the boundary condition (18), C I =0, then: Δh f,0 =C K K0(βr) (27) According to formula (24), we can know that C K expression: Substitute equation (29) into equation (23): in, Define the amplitude ratio: The phase difference is:

4. The method for inverting parameters of fractured aquifers in response to groundwater solid tide according to claim 1, characterized in that: The hydrological borehole water level data of the fractured aquifer to be estimated is obtained, and the amplitude ratio and phase difference of the water head of the fractured aquifer under the action of solid tidal stress are obtained through the baytap solid tidal analysis program.

5. A groundwater solid tide response fractured aquifer parameter inversion system, applied to the groundwater solid tide response fractured aquifer parameter inversion method according to any one of claims 1 to 4, characterized in that: The system includes: a model building unit, a data acquisition and processing unit and a fitting wiring unit; The model building unit is used to establish a dual-medium mathematical model that considers the groundwater response to solid tides in the pore-fracture system of the fractured aquifer, solve the dual-medium mathematical model, and obtain analytical solutions for the head amplitude ratio and phase difference of the fractured aquifer under tidal stress; The data acquisition and processing unit is used to obtain the hydrological borehole water level data of the fractured aquifer to be estimated, and obtain the amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress change; The fitting wiring unit is used to fit the analytical solution with the obtained amplitude ratio and phase difference characteristic parameters of the water level response to the solid tide stress change, and select the determination coefficient R 2 The parameter value corresponding to the largest fitting curve is taken as the inversion result of the fracture aquifer parameters.

Images