CN117724167B - Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body - Google Patents
Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body Download PDFInfo
- Publication number
- CN117724167B CN117724167B CN202410172603.3A CN202410172603A CN117724167B CN 117724167 B CN117724167 B CN 117724167B CN 202410172603 A CN202410172603 A CN 202410172603A CN 117724167 B CN117724167 B CN 117724167B
- Authority
- CN
- China
- Prior art keywords
- velocity
- longitudinal wave
- physical
- dimensional
- wave velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000002689 soil Substances 0.000 title claims abstract description 39
- 239000011435 rock Substances 0.000 title claims description 15
- 239000006185 dispersion Substances 0.000 claims description 43
- 238000003384 imaging method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000013473 artificial intelligence Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 229910052704 radon Inorganic materials 0.000 claims description 3
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 3
- 238000011065 in-situ storage Methods 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000005553 drilling Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000003325 tomography Methods 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for estimating the physical and mechanical parameters of a near-surface rock-soil body in a three-dimensional evolution rule, which belongs to the technical field of earth detection and information technology. The invention provides a method for estimating the physical and mechanical parameter three-dimensional evolution law of a near-surface rock-soil body, which can rapidly and nondestructively acquire the physical and mechanical parameter of the near-surface three-dimensional regional rock-soil body in situ.
Description
Technical Field
The invention belongs to the technical field of earth detection and information, and particularly relates to a method for estimating a three-dimensional evolution rule of a physical and mechanical parameter of a near-surface rock-soil body.
Background
Physical and mechanical parameters are a central concern in near-surface geotechnical engineering construction. Laboratory core measurements and geophysical log interpretation are common methods in the acquisition of the physical mechanical parameters of the rock and earth mass, but they generally provide one-dimensional rock mechanical parameter distribution at the borehole, and are relatively costly, long-term and destructive. How to rapidly and accurately represent the physical and mechanical parameters of the near-surface three-dimensional rock-soil body without damage is a practical problem faced by the current rock-soil body engineering construction research. There is a close relationship between the velocity of the longitudinal and transverse waves of the earthquake and the physical and mechanical parameters of various rock and soil bodies.
The current method for estimating the evolution law of the physical and mechanical parameters of the rock and soil mass mainly comprises a laboratory core measurement method and a geophysical logging data interpretation method. The two methods can accurately provide one-dimensional rock physical and mechanical parameters at the drilling position, but have high cost and long period, and only can reflect static physical and mechanical information at the drilling position, but can not describe regional three-dimensional rock and soil body physical and mechanical change information. How to rapidly and nondestructively acquire the physical and mechanical properties and the change rule of the near-surface three-dimensional rock-soil body is a difficult problem of current engineering geology and hydrogeology research. Aiming at the problem, the invention provides a method for carrying out joint inversion on different period data by adopting time-lapse seismic exploration data to obtain longitudinal wave speed and transverse wave speed so as to estimate the physical and mechanical parameter change of a rock-soil body, and the method can obtain the physical and mechanical parameter change information of the rock-soil body with three-dimensional near-surface high density in a lossless, rapid and accurate manner.
Disclosure of Invention
The invention aims to provide a method for estimating the three-dimensional evolution law of the physical and mechanical parameters of a near-surface rock-soil body, which solves the problems that the traditional calculation method in the prior art only can provide one-dimensional rock-physical and mechanical parameters at a drilling position, has high cost and long period, only can reflect static physical and mechanical information at the drilling position, and can not describe regional three-dimensional rock-soil body physical and mechanical change information.
In order to achieve the purpose, the invention provides a method for estimating the three-dimensional evolution law of the physical and mechanical parameters of a near-surface rock-soil body, which comprises the following steps:
firstly, carrying out multi-period three-dimensional seismic exploration, and then carrying out observation system definition, bad cannon and bad channel rejection pretreatment on the acquired three-dimensional seismic exploration data;
step 2, performing surface wave separation on all the preprocessed period three-dimensional seismic exploration data, and then performing dispersion energy imaging on the separated surface waves and extracting a dispersion curve;
step 3, carrying out first arrival pickup on all the preprocessed secondary three-dimensional seismic exploration data;
step 4, performing surface wave joint inversion according to the dispersion curve of the first-period three-dimensional seismic exploration data extracted in the step 2 to obtain transverse wave speed;
step 5, converting the transverse wave velocity obtained in the step 4 into a longitudinal wave velocity according to the Poisson ratio, and taking the longitudinal wave velocity as an initial longitudinal wave velocity model of the first-stage tomographic inversion, and carrying out the first-stage tomographic inversion to obtain the longitudinal wave velocity by using the first-stage information picked up in the step 3 on the basis of the initial longitudinal wave velocity model;
step 6, evaluating whether the longitudinal wave speed, the transverse wave speed and the existing geological data accord with geological rules or not; if the longitudinal wave velocity is not consistent with the first period, converting the longitudinal wave velocity into the transverse wave velocity to serve as an initial model of the face wave joint inversion in the step 4, and repeating the steps 4-5 until the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule are obtained, and taking the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule as final longitudinal wave velocity and transverse wave velocity in the first period;
step 7, using the longitudinal wave speed and the transverse wave speed obtained in the step 6 as initial models of initial to chromatographic inversion and surface wave inversion in the later period respectively, calculating the longitudinal wave speed and the transverse wave speed in the next period according to the steps 4-6, and repeating the steps to obtain the longitudinal wave speed and the transverse wave speed in all periods;
and 8, estimating Poisson's ratio, compression modulus, elastic modulus and internal friction angle rock-soil physical mechanical parameters of all stages according to the longitudinal wave speed and the transverse wave speed calculated in each stage, and solving the change of adjacent stages.
Preferably, in step 2, the surface wave separation is performed on all the preprocessed three-dimensional seismic exploration data, and then the specific process of performing the dispersive energy imaging on the separated surface waves and extracting the dispersive curve is as follows:
s21, performing surface wave separation by adopting 3DFKK or Lato transformation according to the frequency and apparent speed range of the surface waves, and then cutting off the area beyond the surface waves so as to maximally improve the surface wave signal-to-noise ratio in the three-dimensional seismic data;
s22, performing surface wave phase velocity and group velocity dispersion energy imaging on the separated seismic surface waves by adopting a high-resolution radon transformation, multiple signal classification or time-frequency analysis method;
step S23, picking up a dispersion curve by adopting an artificial intelligence or manual mode on the basis of the dispersion energy imaging data.
Preferably, step 3 comprises the steps of:
step S31, filtering and linear dynamic correction are carried out on three-dimensional seismic exploration data so as to improve the quality of the seismic signals in the first arrival area;
and S32, carrying out earthquake first arrival picking by adopting an artificial picking or artificial intelligent automatic first arrival picking method.
Preferably, in step 4, surface wave joint inversion is performed by using one or more dispersion curves of group velocity and phase velocity, and base and higher order dispersion curves.
Preferably, in step 5, the poisson ratio obtaining mode for converting the transverse wave velocity into the longitudinal wave velocity is a mode of using logging curves, indoor physical and mechanical tests and empirical estimation.
Preferably, in the step 7, the initial model of the initial to tomographic inversion in the later period is the final longitudinal wave speed of the tomographic inversion in the earlier period, and the initial model of the face wave inversion in the later period is the final transverse wave speed of the face wave joint inversion in the earlier period.
Preferably, the physical and mechanical parameters of the rock-soil body in the step 8 are specifically calculated as follows:
step S81, according to the final longitudinal wave speed and transverse wave speed of each period, the ratio of the longitudinal wave speed to the transverse wave speed and the Poisson ratio of each point in the three-dimensional space are obtained:
；
in the middle ofPoisson's ratio->、/>Longitudinal wave velocity and transverse wave velocity respectivelyA degree;
step S82, the density of each point is obtained through a spatial interpolation mode according to a logging curve, an indoor physical mechanical test or an empirical formula, wherein the empirical formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is the density;
step S83, calculating the compression modulus, the dynamic elastic modulus and the dynamic shear modulus according to the longitudinal wave speed, the transverse wave speed and the density obtained in step S82:
(1) Dynamic compression modulusThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>For density (I)>Is transverse wave velocity;
(2) Modulus of dynamic elasticityThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is transverse wave velocity;
(3) Dynamic shear modulusThe calculation formula is as follows:
；
in the middle ofFor transverse wave velocity, < >>Is the density;
step S84, obtaining an internal friction angle according to the Poisson ratio calculated in the step S81:
；
in the middle ofIs poisson's ratio.
Therefore, the method for estimating the physical and mechanical parameters of the near-surface rock and soil body by adopting the three-dimensional evolution rule is adopted, and longitudinal wave speed and transverse wave speed are synchronously obtained by adopting three-dimensional geological exploration data for joint inversion so as to estimate the physical and mechanical parameters of the rock and soil body. The method can rapidly and nondestructively acquire the in-situ physical mechanical parameters of the near-surface three-dimensional regional rock-soil body, provides reliable basis for subsequent geotechnical engineering construction and disaster evaluation, and has great practical value and innovation.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
FIG. 1 is an overall flow chart of a method for estimating the three-dimensional evolution law of the physical and mechanical parameters of a near-surface rock-soil body;
FIG. 2 is a separated face view of an embodiment of the present invention;
FIG. 3 is a graph of a scattered energy imaging of an embodiment of the present invention;
FIG. 4 is a graph of the picked-up dispersion according to an embodiment of the present invention;
FIG. 5 is a first arrival view of a picked up seismic in accordance with an embodiment of the invention;
FIG. 6 is a three-dimensional plot of the final estimated elastic modulus of an embodiment of the present invention;
FIG. 7 is a three-dimensional plot of the dynamic elastic modulus of an embodiment of the present invention, and horizontal slices at different depths.
Detailed Description
The following detailed description of the embodiments of the invention, provided in the accompanying drawings, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1-6, a method for estimating the three-dimensional evolution law of the physical and mechanical parameters of a near-surface rock-soil body comprises the following steps:
firstly, carrying out multi-period three-dimensional seismic exploration, and then carrying out observation system definition, bad cannon and bad channel rejection pretreatment on the acquired three-dimensional seismic exploration data;
step 2, performing surface wave separation on all the preprocessed period three-dimensional seismic exploration data, and then performing dispersion energy imaging on the separated surface waves and extracting a dispersion curve, wherein the specific process is as follows:
s21, performing surface wave separation by adopting 3DFKK or Lato transformation according to the frequency and apparent speed range of the surface waves, and then cutting off the area beyond the surface waves so as to maximally improve the surface wave signal-to-noise ratio in the three-dimensional seismic data;
s22, performing surface wave phase velocity and group velocity dispersion energy imaging on the separated seismic surface waves by adopting a high-resolution radon transformation, multiple signal classification or time-frequency analysis method;
s23, picking up a dispersion curve by adopting an artificial intelligence or manual mode on the basis of the dispersion energy imaging data;
step 3, carrying out first arrival pickup on all the preprocessed secondary three-dimensional seismic exploration data, wherein the method comprises the following steps of:
step S31, filtering and linear dynamic correction are carried out on three-dimensional seismic exploration data so as to improve the quality of the seismic signals in the first arrival area;
step S32, carrying out earthquake first arrival picking by adopting an artificial picking or artificial intelligent automatic first arrival picking method;
step 4, performing surface wave joint inversion according to the dispersion curve of the first-period three-dimensional seismic exploration data extracted in the step 2 to obtain transverse wave speed; the surface wave joint inversion is carried out by adopting one or more dispersion curves in group velocity and phase velocity, base order and high-order dispersion curves;
step 5, converting the transverse wave velocity obtained in the step 4 into a longitudinal wave velocity according to the Poisson ratio, and taking the longitudinal wave velocity as an initial longitudinal wave velocity model of the first-stage tomographic inversion, and carrying out the first-stage tomographic inversion to obtain the longitudinal wave velocity by using the first-stage information picked up in the step 3 on the basis of the initial longitudinal wave velocity model; the Poisson ratio acquisition mode for converting the transverse wave speed into the longitudinal wave speed is a mode of adopting a logging curve, an indoor physical mechanical test and an empirical estimation;
step 6, evaluating whether the longitudinal wave speed, the transverse wave speed and the existing geological data accord with geological rules or not; if the longitudinal wave velocity is not consistent with the first period, converting the longitudinal wave velocity into the transverse wave velocity to serve as an initial model of the face wave joint inversion in the step 4, and repeating the steps 4-5 until the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule are obtained, and taking the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule as final longitudinal wave velocity and transverse wave velocity in the first period;
step 7, using the longitudinal wave speed and the transverse wave speed obtained in the step 6 as initial models of initial to chromatographic inversion and surface wave inversion in the later period respectively, calculating the longitudinal wave speed and the transverse wave speed in the next period according to the steps 4-6, and repeating the steps to obtain the longitudinal wave speed and the transverse wave speed in all periods; the initial model of the first-arrival tomography inversion in the later period is the final longitudinal wave speed of the tomography inversion in the former period, and the initial model of the face wave inversion in the later period is the final transverse wave speed of the face wave joint inversion in the former period;
step 8, estimating Poisson's ratio, compression modulus, elastic modulus and internal friction angle rock-soil physical mechanical parameters of all stages according to the longitudinal wave speed and the transverse wave speed calculated in each stage, and solving the change of adjacent stages; the physical and mechanical parameters of the rock-soil body are calculated as follows:
step S81, according to the final longitudinal wave speed and transverse wave speed of each period, the ratio of the longitudinal wave speed to the transverse wave speed and the Poisson ratio of each point in the three-dimensional space are obtained:
；
in the middle ofPoisson's ratio->、/>Longitudinal wave velocity and transverse wave velocity respectively;
step S82, the density of each point is obtained through a spatial interpolation mode according to a logging curve, an indoor physical mechanical test or an empirical formula, wherein the empirical formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is the density;
step S83, calculating the compression modulus, the dynamic elastic modulus and the dynamic shear modulus according to the longitudinal wave speed, the transverse wave speed and the density obtained in step S82:
(1) Dynamic compression modulusThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>For density (I)>Is transverse wave velocity;
(2) Modulus of dynamic elasticityThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is transverse wave velocity;
(3) Dynamic shear modulusThe calculation formula is as follows:
；
in the middle ofFor transverse wave velocity, < >>Is the density;
step S84, obtaining an internal friction angle according to the Poisson ratio calculated in the step S81:
；
in the middle ofIs poisson's ratio.
Examples
The flow of the method is described by taking four-dimensional seismic exploration data of a certain coal field as an example.
Step 1, four-period three-dimensional seismic exploration is carried out in the coal field research area (the four-period three-dimensional seismic data adopts the same position and the same observation system), and the adjacent interval time is about half a year. Performing pretreatment such as definition of an observation system, removal of bad cannons and bad tracks on the four-stage three-dimensional seismic data;
and 2, setting the frequency range to be 8-50 hz and the apparent velocity range to be 200-1500 m/s, and cutting off the area beyond the surface wave after performing surface wave separation on the preprocessed seismic data according to the set surface wave frequency and apparent velocity range by adopting 3DFKK so as to maximally improve the surface wave signal-to-noise ratio in the three-dimensional seismic data. On the basis, a multiple signal classification method and a time-frequency analysis method are adopted to respectively image phase velocity and group velocity dispersion energy, and finally, the phase velocity and group velocity dispersion curve of the surface wave is manually picked up. The separated surface wave is shown in fig. 2. The dispersive energy imaging diagram is shown in fig. 3. The picked-up surface dispersion curve is shown in fig. 4. Fig. 2 shows a face wave obtained by performing face wave separation and ablation on seismic data by using a crisscross domain 3 DFKK. From this figure, the rayleigh wave signal with strong reflected amplitude and obvious linear characteristic can be clearly seen. Fig. 3 is a graph of phase velocity dispersion energy imaging using a multiple signal classification method on the data of fig. 2. From this figure, two distinct dispersive energy clusters can be seen, of which the left is the fundamental mode rayleigh wave dispersive energy; on the right is the first Gao Jierui rayleigh dispersed energy. Fig. 4 is a graph of the dispersion energy picked up according to the maximum value of the dispersion energy of fig. 3, wherein picking up the maximum value of the dispersion energy of the rayleigh wave of the left side mode of fig. 3 can obtain a fundamental mode dispersion curve, and picking up the maximum value of the dispersion energy of the rayleigh wave of the first higher order mode of fig. 3 can obtain a first higher order mode dispersion curve. From the dispersion curve picked up in FIG. 3, it can be seen that the phase velocity of the point base mode corresponding to the shallow high frequency is about 280m/s, and the phase velocity of the low frequency end corresponding to the deep is about 900m/s; the lowest phase velocity in the first higher order mode is 550m/s, and the highest phase velocity is greater than 1000m/s. The high-quality and multi-mode dispersion curve provides a good basis for subsequent face wave inversion.
And step 3, manually picking up the first arrival of the preprocessed seismic data by adopting a manual picking method. The first arrival of the picked earthquake is shown in fig. 5, and the white point in fig. 5 is the first arrival time of the earthquake wave picked by adopting the manual picking method, and the subsequent refraction wave tomographic inversion can be performed by using the first arrival time.
And 4, dividing the whole research area into 3 typical areas according to 35 geological drilling data in the area, the logging curve and the dispersion curve change rule picked up in the step 2.
And 5, selecting 1 characteristic point (a drilling hole is arranged near the characteristic point and the quality of the phase velocity and group velocity dispersion curve is preferably used as a characteristic point selection standard) in each typical area, and performing joint inversion on the first-period surface dispersion curve by adopting a Rayleigh wave phase velocity and group velocity joint inversion method. Taking the transverse wave velocity obtained by the joint inversion as a plane wave inversion initial model of all points in the region;
step 6, inverting the phase velocity of the dispersion curves of all points in the area by using the initial model obtained in the step 5 to obtain the transverse wave velocity;
step 7, smoothing the transverse wave speed and converting the transverse wave speed into a longitudinal wave speed according to the Poisson ratio of the well logging curve;
step 8, performing tomographic inversion by taking the longitudinal wave velocity as an initial model and utilizing the seismic first arrival information picked up in the step 3 to obtain the longitudinal wave velocity of the near surface;
and 9, extracting the longitudinal wave velocity and transverse wave velocity horizontal slice and vertical slice obtained in the step 6 and the step 8, and finding out that the detail description of the local area and the existing geological information is not clear. The longitudinal wave velocity of step 8 is converted into the transverse wave velocity according to the poisson ratio of step 7 and is used as an initial model of the face wave inversion again. And repeating the steps 6-8 on the basis to obtain more reasonable longitudinal wave speed and transverse wave speed.
And 10, respectively taking the longitudinal wave speed and the transverse wave speed obtained in the step 9 as initial models of second-stage surface wave inversion and first-arrival tomography inversion, and calculating the longitudinal wave speed and the transverse wave speed in the second stage according to the steps 6-8. And the like to obtain longitudinal wave speeds and transverse wave speeds of all four periods.
And step 11, calculating physical and mechanical parameters of the rock-soil body such as poisson ratio, compression modulus, elastic modulus, internal friction angle and the like according to the four-stage longitudinal wave speed and the transverse wave speed in the step 10. The three-dimensional plot of the final estimated modulus of elasticity is shown in FIG. 6. A three-dimensional plot of the final estimated elastic modulus and horizontal slices at different depths are shown in fig. 7. Fig. 7 is a three-dimensional graph of elastic modulus calculated by using longitudinal wave velocity and transverse wave velocity, from which three-dimensional spatial distribution characteristics of elastic modulus of the whole region can be seen, horizontal slices with a height of 1300-1250 m and an interval of 10 m are further extracted, from the horizontal slices, the elastic modulus of the research region can be divided into three regions in total, the elastic modulus of the north-east side and the south-west side is relatively small, the elastic modulus of the middle region is relatively large, and the whole is in north-south east trend. Further comparing horizontal slices with different depths, the elastic modulus area with larger middle part becomes larger along with the deeper depth, which is consistent with the general trend of structural fluctuation of the bedrock surface of the area, and the estimation of physical and mechanical parameters, particularly the elastic modulus estimation, of the seismic data adopted by the scheme is reasonable and reliable from the side surface.
Therefore, the method for estimating the physical and mechanical parameters of the near-surface rock and soil body by adopting the three-dimensional evolution rule is adopted, and longitudinal wave speed and transverse wave speed are synchronously obtained by adopting three-dimensional geological exploration data for joint inversion so as to estimate the physical and mechanical parameters of the rock and soil body. The method can rapidly and nondestructively acquire the in-situ physical mechanical parameters of the near-surface three-dimensional regional rock-soil body, provides reliable basis for subsequent geotechnical engineering construction and disaster evaluation, and has great practical value and innovation.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention and not for limiting it, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that: the technical scheme of the invention can be modified or replaced by the same, and the modified technical scheme cannot deviate from the spirit and scope of the technical scheme of the invention.
Claims (7)
1. The method for estimating the three-dimensional evolution rule of the physical and mechanical parameters of the near-surface rock and soil mass is characterized by comprising the following steps of:
firstly, carrying out multi-period three-dimensional seismic exploration, and then carrying out observation system definition, bad cannon and bad channel rejection pretreatment on the acquired three-dimensional seismic exploration data;
step 2, performing surface wave separation on all the preprocessed period three-dimensional seismic exploration data, and then performing dispersion energy imaging on the separated surface waves and extracting a dispersion curve;
step 3, carrying out first arrival pickup on all the preprocessed secondary three-dimensional seismic exploration data;
step 4, performing surface wave joint inversion according to the dispersion curve of the first-period three-dimensional seismic exploration data extracted in the step 2 to obtain transverse wave speed;
step 5, converting the transverse wave velocity obtained in the step 4 into a longitudinal wave velocity according to the Poisson ratio, and taking the longitudinal wave velocity as an initial longitudinal wave velocity model of the first-stage tomographic inversion, and carrying out the first-stage tomographic inversion to obtain the longitudinal wave velocity by using the first-stage information picked up in the step 3 on the basis of the initial longitudinal wave velocity model;
step 6, evaluating whether the longitudinal wave speed, the transverse wave speed and the existing geological data accord with geological rules or not; if the longitudinal wave velocity is not consistent with the first period, converting the longitudinal wave velocity into the transverse wave velocity to serve as an initial model of the face wave joint inversion in the step 4, and repeating the steps 4-5 until the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule are obtained, and taking the longitudinal wave velocity and the transverse wave velocity consistent with the geological rule as final longitudinal wave velocity and transverse wave velocity in the first period;
step 7, using the longitudinal wave speed and the transverse wave speed obtained in the step 6 as initial models of initial to chromatographic inversion and surface wave inversion in the later period respectively, calculating the longitudinal wave speed and the transverse wave speed in the next period according to the steps 4-6, and repeating the steps to obtain the longitudinal wave speed and the transverse wave speed in all periods;
and 8, estimating Poisson's ratio, compression modulus, elastic modulus and internal friction angle rock-soil physical mechanical parameters of all stages according to the longitudinal wave speed and the transverse wave speed calculated in each stage, and solving the change of adjacent stages.
2. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock-soil body according to claim 1, wherein the specific process of performing surface wave separation on all the preprocessed three-dimensional seismic exploration data, then performing frequency dispersion energy imaging on the separated surface waves and extracting a frequency dispersion curve in the step 2 is as follows:
s21, performing surface wave separation by adopting 3DFKK or Lato transformation according to the frequency and apparent speed range of the surface wave, and then cutting off the area beyond the surface wave;
s22, performing surface wave phase velocity and group velocity dispersion energy imaging on the separated seismic surface waves by adopting a high-resolution radon transformation, multiple signal classification or time-frequency analysis method;
step S23, picking up a dispersion curve by adopting an artificial intelligence or manual mode on the basis of the dispersion energy imaging data.
3. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock-soil body according to claim 1, wherein the step 3 comprises the following steps:
step S31, filtering and linear dynamic correction are carried out on three-dimensional seismic exploration data;
and S32, carrying out earthquake first arrival picking by adopting an artificial picking or artificial intelligent automatic first arrival picking method.
4. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock-soil body according to claim 1, which is characterized by comprising the following steps of: and 4, performing surface wave joint inversion by adopting one or more dispersion curves of group velocity and phase velocity, base order and high-order dispersion curves.
5. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock-soil body according to claim 1, which is characterized by comprising the following steps of: in step 5, the poisson ratio obtaining mode for converting the transverse wave speed into the longitudinal wave speed is a mode of adopting a logging curve, an indoor physical and mechanical test and an empirical estimation.
6. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock-soil body according to claim 1, which is characterized by comprising the following steps of: in the step 7, the initial model of the initial to tomographic inversion in the later period is the final longitudinal wave speed of the tomographic inversion in the earlier period, and the initial model of the face wave inversion in the later period is the final transverse wave speed of the face wave joint inversion in the earlier period.
7. The method for estimating the three-dimensional evolution law of the physical and mechanical parameters of the near-surface rock and soil mass according to claim 1, wherein the concrete calculation of the physical and mechanical parameters of the rock and soil mass in the step 8 is as follows:
step S81, according to the final longitudinal wave speed and transverse wave speed of each period, the ratio of the longitudinal wave speed to the transverse wave speed and the Poisson ratio of each point in the three-dimensional space are obtained:
；
in the middle ofPoisson's ratio->、/>Longitudinal wave velocity and transverse wave velocity respectively;
step S82, the density of each point is obtained through a spatial interpolation mode according to a logging curve, an indoor physical mechanical test or an empirical formula, wherein the empirical formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is the density;
step S83, calculating the compression modulus, the dynamic elastic modulus and the dynamic shear modulus according to the longitudinal wave speed, the transverse wave speed and the density obtained in step S82:
(1) Dynamic compression modulusThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>For density (I)>Is transverse wave velocity;
(2) Modulus of dynamic elasticityThe calculation formula is as follows:
；
in the middle ofFor longitudinal wave velocity +.>Is transverse wave velocity;
(3) Dynamic shear modulusThe calculation formula is as follows:
；
in the middle ofFor transverse wave velocity, < >>Is the density;
step S84, obtaining an internal friction angle according to the Poisson ratio calculated in the step S81:
；
in the middle ofIs poisson's ratio.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410172603.3A CN117724167B (en) | 2024-02-07 | 2024-02-07 | Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410172603.3A CN117724167B (en) | 2024-02-07 | 2024-02-07 | Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117724167A CN117724167A (en) | 2024-03-19 |
| CN117724167B true CN117724167B (en) | 2024-04-12 |
Family
ID=90207284
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410172603.3A Active CN117724167B (en) | 2024-02-07 | 2024-02-07 | Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117724167B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016187252A1 (en) * | 2015-05-20 | 2016-11-24 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
| CN111538075A (en) * | 2020-05-11 | 2020-08-14 | 中国地质调查局水文地质环境地质调查中心 | Hot dry rock exploration method and device, electronic equipment and storage medium |
| CN117452491A (en) * | 2023-08-21 | 2024-01-26 | 四川省自然资源投资集团物探勘查院有限公司 | Combined exploration method for identifying characteristics of gas reservoirs of coal series under complicated mountain land surface conditions |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11947063B2 (en) * | 2022-04-28 | 2024-04-02 | Saudi Arabian Oil Company | Method of conditioning seismic data for first-break picking using nonlinear beamforming |
-
2024
- 2024-02-07 CN CN202410172603.3A patent/CN117724167B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016187252A1 (en) * | 2015-05-20 | 2016-11-24 | Conocophillips Company | Surface wave tomography using sparse data acquisition |
| CN111538075A (en) * | 2020-05-11 | 2020-08-14 | 中国地质调查局水文地质环境地质调查中心 | Hot dry rock exploration method and device, electronic equipment and storage medium |
| CN117452491A (en) * | 2023-08-21 | 2024-01-26 | 四川省自然资源投资集团物探勘查院有限公司 | Combined exploration method for identifying characteristics of gas reservoirs of coal series under complicated mountain land surface conditions |
Non-Patent Citations (2)
| Title |
|---|
| Estimating picking errors in near‐surface seismic data to enable their time‐lapse interpretation of hydrosystems;M. Dangeard 等;Near Surface Geophysics;20181231;第16卷(第6期);613-625 * |
| 面波信息约束的初至波走时层析反演方法;张利振 等;物探与化探;20231031;第47卷(第5期);1198-1205 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117724167A (en) | 2024-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108931814B (en) | Multi-attribute fusion based bedrock fracture prediction method | |
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| WO2016041189A1 (en) | Method for evaluating shale gas reservoir and seeking desert area | |
| CN106646601B (en) | The shallow mid-deep strata three-dimensional Q body method for building up of multi information joint constraint | |
| CN105277982A (en) | Shale total organic carbon content earthquake prediction method | |
| CN111175815B (en) | Method and system for solving micro-seismic monitoring crack seismic source mechanism in oil reservoir transformation | |
| CN105629308A (en) | Phase control heterogeneous mechanical parameter crustal stress prediction method | |
| CN107678064B (en) | Real-time extraction method for sound wave time difference | |
| CN112946752B (en) | Method for predicting fracture probability body based on fracture density model | |
| CN109541689B (en) | Method for evaluating compactness of medium based on reflected wave energy characteristics | |
| CN112578435A (en) | Rock ultrasonic testing first arrival picking method and system | |
| CN111308558B (en) | Shale gas horizontal well longitudinal wave time difference correction method | |
| CN109188528A (en) | Elastic wave chromatographic imaging system and method between well | |
| CN117724167B (en) | Estimation method for physical and mechanical parameter three-dimensional evolution rule of near-surface rock and soil body | |
| CN110244383B (en) | Geological lithology comprehensive model establishing method based on near-surface data | |
| CN106291675A (en) | A kind of geological data reconstructing method based on base tracer technique | |
| CN107679614B (en) | Particle swarm optimization-based real-time sound wave time difference extraction method | |
| CN110989034B (en) | Method for inverting logging transverse wave time difference by regression-fractal interpolation method | |
| CN113050168B (en) | Crack effectiveness evaluation method based on array acoustic logging and acoustic remote detection logging data | |
| CN106291676A (en) | A kind of geological data reconstructing method based on matching pursuit algorithm | |
| CN108363739A (en) | A kind of seismic data low-and high-frequency expanding method based on sparse acquisition | |
| CN110632660B (en) | Thin sand body characterization method and device based on seismic data body | |
| CN113296149A (en) | Landslide mass stratum distribution condition detection method based on micro-motion detection technology | |
| CN110297264B (en) | Low-permeability gas reservoir thin reservoir dessert earthquake prediction method | |
| CN112906465A (en) | Coal measure stratum acoustic curve reconstruction method and system based on stratum factors |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |