CN115542393A - Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction - Google Patents
Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction Download PDFInfo
- Publication number
- CN115542393A CN115542393A CN202211254264.0A CN202211254264A CN115542393A CN 115542393 A CN115542393 A CN 115542393A CN 202211254264 A CN202211254264 A CN 202211254264A CN 115542393 A CN115542393 A CN 115542393A
- Authority
- CN
- China
- Prior art keywords
- model
- correction
- inversion
- updating
- waveform inversion
- 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.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001514 detection method Methods 0.000 claims abstract description 28
- 238000003825 pressing Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 9
- 230000002159 abnormal effect Effects 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 10
- 238000009412 basement excavation Methods 0.000 description 9
- 230000005641 tunneling Effects 0.000 description 7
- 239000003245 coal Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a roadway driving full-waveform inversion method based on multi-parameter constraint and structural correction, which comprises the following steps of: step 1, constructing a full-waveform inversion initial model, and performing single-scale inversion on the initial model by adopting a multi-scale elastic wave full-waveform inversion method; step 2, carrying out multi-parameter weighted constraint structure correction on the single-scale inversion result to obtain a primary correction result; step 3, based on preset limiting conditions, carrying out one-dimensional wave velocity profile space structure correction and smooth constraint on the primary correction result to obtain a secondary correction result; step 4, taking the secondary correction result as an initial model of the next scale, and continuing to perform full waveform inversion; and (4) repeating the step (2) to the step (4) until all scales of inversion are finished, and obtaining a full waveform inversion result. The method can effectively solve the problems of strong limitation of an observation system, small detection data amount, small offset distance and strong waveform inversion multi-solution in the tunneling-following roadway advanced detection, and improves the inversion effect.
Description
Technical Field
The invention relates to the technical field of mine exploration imaging, in particular to an along-excavation roadway full-waveform inversion method based on multi-parameter constraint and structural correction, and more particularly relates to an along-excavation roadway detection elastic-waveform inversion method based on multi-parameter weighting constraint and wave velocity profile space structural correction.
Background
As one of two core links of coal mine production, the development requirement of intelligent tunneling is very urgent, but in the process of tunneling, disasters such as coal and gas outburst, water burst and the like seriously threaten the safety of tunneling production and the personal safety of miners. The geological guarantee technology is the basis of coal intelligent production safety guarantee, is a basic data source for realizing geological prediction, disturbance perception and risk assessment before, during and after tunnel excavation construction, and is a precondition guarantee for implementing all intelligent excavation key technologies.
However, most of the current mine earthquake advanced detection imaging interpretation results are 'arc drawing', many false abnormal interfaces are provided, the imaging precision is low (as shown in figure 1), and the requirement of intelligent roadway tunneling geological guarantee is difficult to meet.
The full waveform inversion method can fully utilize the kinematics and the dynamics characteristics of seismic waves to obtain the parameter information of the underground model, has the advantages of high imaging precision of complex structures, good physical parameter inversion effect and the like, has good application effect in surface seismic exploration, is the optimal selection of advanced detection imaging of mine seismic in the future, and can meet the requirements of intelligent tunneling geological guarantee of a roadway. However, the earthquake advanced detection mode is different from the earth surface detection, the limitation of an observation system is strong, a series of linearly arranged shot points and wave detection points are only distributed on a central line behind a tunneling head, the abnormal body in front of the tunneling is similar to single offset detection, and the quantity of the distributed shot points and wave detection points is limited by a detection space, so that the detection mode has small data quantity and small offset, the multi-solution of full waveform inversion is enhanced, and the calculation error of physical property parameters is increased.
Therefore, how to provide a high-precision full-waveform inversion method suitable for a mine roadway detection mode along with excavation is a problem to be solved urgently by the technical personnel in the field.
Disclosure of Invention
In view of the above, the invention provides a roadway driving full-waveform inversion method based on multi-parameter constraint and structure correction, which combines multi-parameter weighting constraint and wave velocity profile space structure correction, effectively solves the problems of strong limitation of an observation system, small detection data amount, small offset distance and strong waveform inversion multi-solution in roadway advance detection, improves the inversion effect, and realizes high-precision imaging of an abnormal geological structure in the front of driving.
In order to achieve the purpose, the invention adopts the following technical scheme:
a roadway driving full-waveform inversion method based on multi-parameter constraint and structure correction comprises the following steps:
step 2, carrying out multi-parameter weighted constraint structure correction on the single-scale inversion result obtained in the step 1 to obtain a primary correction result;
step 3, based on preset limiting conditions, carrying out one-dimensional wave velocity profile space structure correction and smooth constraint on the primary correction result to obtain a secondary correction result;
step 4, taking the secondary correction result obtained in the step 3 as an initial model of the next scale, and continuing to perform full waveform inversion according to the mode of the step 1;
and 5, repeatedly executing the step 2 to the step 4 until all scales of inversion are finished, and obtaining a full waveform inversion result of the roadway advanced detection elastic wave.
Further, in step 2, correcting the multi-parameter weighting constraint structure according to the following formula:
wherein, Δ m is the updating amount of the single iteration model; Δ m' corrected for multi-parameter weighting structureThe amount of model updates; m is an initial model; i and j respectively represent the positions of the grid nodes in the z and x directions; m is i,j p The parameter values of the positions of the nodes of the i, j grids of the initial model are single parameters of longitudinal wave velocity, transverse wave velocity or density; nz is the number of model vertical grid points, nx is the number of model horizontal grid points; v p Is the longitudinal wave velocity; v s Is the transverse wave velocity; den is the density.
Further, in step 3, the preset limiting conditions are:
setting a first limiting condition according to actual needs, wherein the first limiting condition takes the spatial correction lower limit value as a judgment basis, and updating the model with the value lower than the spatial correction lower limit value
Pressing to 0;
and setting a second limiting condition according to actual needs, wherein the second limiting condition takes the range of the structural correction area and the distance between the areas as a judgment basis, when the model updating amount is a negative value area, both sides of the area have an updating amount which is a positive value area, or when the model updating amount is a positive value area, both sides of the area have an updating amount which is a negative value area, and when the area range and the distance between the areas meet the second limiting condition, the model updating amount in the positive value area existing on both sides of the negative value area or the negative value area existing on both sides of the positive value area is suppressed to be 1/5 of the original model updating value.
Further, step 3 comprises:
301, selecting a grid coordinate y of a longitudinal axis i Extracting one-dimensional wave velocity profile along the axis of the tunnel
Wherein x is a horizontal axis coordinate;
step 302, calculating a one-dimensional wave velocity profile
Corresponding model update amount
And slope of
Step 303, updating the model below the value of the lower value of the spatial correction limit
Pressing to 0 to obtain corrected model update amount
Step 304, updating the quantity according to the corrected model
Judging and correcting each grid point on the two-dimensional wave velocity profile one by one according to the limiting conditions to obtain the corrected model updating amount
Step 305, updating quantity according to corrected model
And slope
Judging each grid point on the one-dimensional wave velocity profile one by one according to the change condition of the slope sign to obtain the corrected model update quantity
Step 306, updating the quantity according to the corrected model
Calculating a corrected one-dimensional wave velocity profile
Step 307, take down a grid coordinate y of the vertical axis i+1 And repeating the step 302 to the step 306 until the correction of all the grid one-dimensional wave velocity profiles is completed, and obtaining a secondary correction result after the space structure correction and the smoothing constraint of the one-dimensional wave velocity profiles.
Further, step 305 includes:
when two points with zero model updating quantity exist between three adjacent points with changed slope symbols, the model updating quantity symbols of grid points before and after the two points are changed, and no point with 0 slope exists between the two points, the model updating quantity of the grid points in the area between the three adjacent points with changed slope symbols is corrected to be the average value of the model updating quantities of the two points of the first slope symbol changing point and the third slope symbol changing point, and the corrected model updating quantity is obtained
Then according to the updated quantity of the corrected model
Continuously correcting, when the sign of the model updating quantity of all the grid points in the range between two nonadjacent points with the changed slope sign is not changed, correcting the model updating quantity of the grid points in the region between the two changed slope sign points, changing the model updating quantity into the average value of the model updating quantity of all the grid points in the region, and obtaining the corrected model updating quantity
Further, the inversion result in step 5 includes: longitudinal wave inversion results, shear wave inversion results and density inversion results.
According to the technical scheme, compared with the prior art, the invention discloses a method for inverting the full waveform of the roadway along with the excavation based on multi-parameter constraint and structural correction, which has the following beneficial effects:
(1) The invention adopts an elastic wave full waveform inversion method, utilizes longitudinal wave velocity, transverse wave velocity and density parameters in observation records, and can obtain more accurate underground medium information than acoustic wave full waveform inversion.
(2) According to the invention, through carrying out multi-parameter weighted constraint structure correction and one-dimensional wave velocity profile space structure correction and smooth constraint on the inversion result under a single scale, twice correction is carried out, high-precision imaging of the abnormal geological structure in front of roadway excavation is realized, the position and the occurrence of a geological abnormal body can be accurately judged by the imaging result, and the blank of full waveform inversion technology in the field of mine advanced detection is filled.
(3) The inversion strategy of the invention not only provides an initial model with higher precision for the next-scale inversion, but also avoids the problem of excessive aggravation of multiple iteration false anomalies, and the false anomalies are eliminated by structural correction, so that the direction of the full waveform inversion is restrained.
(4) The method basically solves the problems of small detection data volume, small offset distance and strong waveform inversion multi-solution in the advance detection of the roadway, improves the inversion effect, and improves the full waveform inversion accuracy by about 20 percent.
(5) The invention can improve the inversion accuracy without increasing the calculation amount of full waveform inversion.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a conventional mine seismic advanced survey imaging result provided by the present invention;
FIG. 2 is a flow chart of a method for full waveform inversion of detection elastic waves along with a roadway driving provided by the invention;
FIG. 3 is a schematic diagram illustrating the determination of the update amount and slope change of the one-dimensional wave velocity profile model according to the present invention;
FIG. 4 is a schematic structural diagram of an initial model provided by the present invention;
FIG. 5 is a complex advanced geological theoretical model provided by the present invention;
FIG. 6 is a diagram illustrating inversion results obtained by the method of the present invention;
FIG. 7 is a full waveform inversion result of a conventional time domain multi-scale elastic wave.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
As shown in fig. 2, the embodiment of the invention discloses a method for performing full waveform inversion along with a roadway based on multi-parameter constraint and structural correction, which comprises the following steps:
Step 2, correcting the single-scale inversion result obtained in the step 1 by a multi-parameter weighting constraint structure to obtain a primary correction result;
step 3, based on preset limiting conditions, carrying out one-dimensional wave velocity profile space structure correction and smooth constraint on the primary correction result to obtain a secondary correction result;
step 4, taking the secondary correction result obtained in the step 3 as an initial model of the next scale, and continuing to perform full waveform inversion according to the mode of the step 1;
and 5, repeatedly executing the step 2 to the step 4 until all scales of inversion are finished, and obtaining a full waveform inversion result of the roadway advanced detection elastic wave. The inversion results include: longitudinal wave inversion results, shear wave inversion results and density inversion results.
In the invention, the seismic sources are distributed in the middle of the coal seam, a plurality of detectors with horizontal components and vertical components are linearly arranged with the seismic sources and are also distributed in the middle of the coal seam to receive seismic signals, and an inversion initial model is set according to the condition of a uniform layered coal system stratum medium.
Due to the fact that the tunnel is used for detecting special observation system conditions in advance, the multi-solution problem of the conventional time domain multi-scale elastic wave full waveform inversion method is enhanced, and inversion accuracy is difficult to further improve from a data angle. According to the invention, from the model structure, the special structural characteristics implied by geological abnormal bodies in the roadway advanced detection full waveform inversion result are found, the correction term is constructed according to the structural characteristics, the direction of full waveform inversion is constrained, the inversion result meeting the preset structural characteristics is obtained, and the inversion effect is improved.
The discovered special structural characteristics of the geological anomalous body in the tunnel advanced detection elastic wave full waveform inversion result are as follows:
(1) For the recovery of abnormal geological structures, the inversion effect of density parameters is best, and the suppression effect of false abnormality and disturbance interference in longitudinal and transverse wave parameters is better;
(2) Longitudinal, transverse and density waves have the same geologic structure characteristics, although they have large differences in physical parameters;
(3) False abnormal regions opposite to the attribute characteristics of the false abnormal regions exist near the boundaries on the two sides of the abnormal geological structure;
(4) The horizontal speed-depth curve of the grid points in the vertical direction has W or M type abnormal characteristics, and the update sign of the middle range model of the W or M area is opposite to the update sign of the range models on the two sides;
(5) The transverse velocity-depth curve of the grid points in the vertical direction in the large-volume structure has the characteristic of 'wavy lines', and the 'wavy lines' relate to consistent model updating signs in a range.
Based on the found characteristics, the inversion result is corrected for multiple times by introducing multi-parameter weighting constraint and one-dimensional wave velocity profile space structure correction, the problems of small detection data amount, small offset distance and strong waveform inversion multi-solution in roadway advanced detection can be effectively solved, the inversion effect is improved, and high-precision imaging of the abnormal geological structure in the front of the excavation is realized.
Next, the correction method in step 2 and step 3 will be described in detail.
In a specific embodiment, in step 2, the multi-parameter weighted constraint structure is corrected according to the following formula:
wherein, Δ m is the updating amount of the single iteration model (i.e. the variation between the current iteration and the last iteration model); Δ m' is the model update amount after the multi-parameter weighting structure correction; m is an initial model; i and j respectively represent the positions of the grid nodes in the z and x directions; m is i,j p The parameter values of the positions of the nodes of the i, j grids of the initial model are single parameters of longitudinal wave velocity, transverse wave velocity or density; nz is the number of vertical grids of the model, and nx is the number of horizontal grids of the model; v p Is the velocity of the longitudinal wave; v s Is the transverse wave velocity; den is density
In one embodiment, in step 3, the preset limiting conditions are:
the first limiting condition takes the value of the lower limit of the spatial correction as a judgment basis, and updates the model with the value lower than the value of the lower limit of the spatial correction
Pressing to 0; in practical application, the autonomous setting is carried out according to the data processing effect, and is generally determined through experiments before formal inversion. For example, the value is set to less than 5% of the initial model parameters.
And the second limiting condition takes the structural correction area range and the distance between the areas as a judgment basis, when the model updating amount is a negative value area, an updating amount is a positive value area on both sides of the area, or when the model updating amount is a positive value area on both sides of the area, the updating amount is a negative value area on both sides of the area, and when the distance between the area range and the area meets the second limiting condition, the model updating amount in the positive value area on both sides of the negative value area or the negative value area on both sides of the positive value area is suppressed to be 1/5 of the original model updating value.
Wherein, the "distance between the regions" means: when both sides of the negative value area have an updating amount which is a positive value area, the distance between the positive value area at one side and the negative value area at the middle part is the distance between the positive value area at one side and the negative value area at the middle part; or the following steps: when both sides of the positive value area of the update amount have negative value areas of the update amount, the distance between the negative value area on one side and the positive value area in the middle is increased.
"area range" means: when there is a positive value area on both sides of the negative value area, the range of the single-side positive value area is as follows: when there is a region where the update amount is negative on both sides of the region where the update amount is positive, the range of the "negative region" on one side is large.
The extent of the regions and the distance between the regions are generally determined by experiments prior to formal inversion. For example: "area range: when the updating amount is a positive value area on both sides of the negative value area, the range size of the single-side positive value area cannot be smaller than 1/2 of the range of the middle positive value area.
"distance between regions": when the updating amount is a positive value area on both sides of the negative value area, the distance between the single positive value area and the middle negative value area cannot be larger than the size of the middle negative value area.
Specifically, step 3 includes:
301, selecting a grid coordinate y of a longitudinal axis i Extracting one-dimensional wave velocity profile along the axis of the tunnel
Wherein x is a horizontal axis coordinate;
step 302, calculate a one-dimensional wave velocity profile
Corresponding model update amount
And slope
Wherein, the model updating amount is as follows:
is a one-dimensional wave velocity profile at the e-th iteration,
is the one-dimensional wave velocity profile at the e-1 st iteration.
Slope:
is a one-dimensional wave velocity profile
The value of the jth grid point in the lateral direction,
is a one-dimensional wave velocity profile
The value of the j-1 horizontal grid point, dx is the model horizontal grid distance.
Step 303, updating the model below the value of the lower value of the spatial correction limit
Performing suppression to become 0 to obtain the corrected model update amount
Step 304, updating the quantity according to the corrected model
Judging and correcting each grid point on the two-dimensional wave velocity profile one by one according to the limiting conditions to obtain the corrected model updating amount
Step 305, updating amount according to corrected model
And slope
Judging each grid point on the one-dimensional wave velocity profile one by one according to the change condition of the slope sign to obtain the corrected model update quantity
When two points (as shown in fig. 3) with zero model update quantity exist between three adjacent points with changed slope signs, the model update quantity signs of the grid points before and after the two points change, and no point with 0 slope exists between the two points, the model update quantity of the grid point in the region between the three adjacent points with changed slope signs is corrected to become the average value of the model update quantities of the two points at the first slope sign change point and the third slope sign change point, and the corrected model update quantity is obtained
Then according to the updated quantity of the corrected model
The correction is continued and the correction is continued,when the model updating quantity symbols of all grid points in the range between two nonadjacent points with changed slope symbols are not changed, correcting the model updating quantity of the grid points in the region between the two changed slope symbols to obtain the average value of the model updating quantity of all grid points in the region, and obtaining the corrected model updating quantity
In order to control the correction range, search determination is performed according to a set limit value of the area range. That is, the maximum correction range is controlled, for example, to be set to 10 grid points, and the "search determination is performed according to the set area range limit" is to sequentially search for 10 grid points in one turn.
Step 306, updating the quantity according to the corrected model
Calculating a corrected one-dimensional wave velocity profile
Step 307, take the grid coordinate y of the next vertical axis i+1 And repeating the step 302 to the step 306 until the correction of all the grid one-dimensional wave velocity profiles is completed, and obtaining a secondary correction result after the spatial structure correction and the smoothing constraint of the one-dimensional wave velocity profiles.
In order to further verify the effectiveness and the high efficiency of the tunnel-along-excavation exploration elastic wave full waveform inversion method based on multi-parameter weighting constraint and wave velocity profile space structure correction, the scheme provided by the invention is used in a complex advanced geological theoretical model (figure 5), and the inversion result is obtained as shown in figure 6.
Comparing fig. 6 and 7, it can be seen that the full waveform inversion result obtained by the present invention is closer to the real model, taking the conventional time domain multi-scale elastic wave full waveform inversion result (as shown in fig. 7) as the comparison study object. The boundaries of two sides of the structure of the collapse column and the fault fracture zone in the inversion result are clear, the recovery of the internal parameters of the structure is good, and the lithologic boundary of the small-size fault layer is recovered and displayed; most obviously, the false abnormity and disturbance interference in the result are basically completely suppressed, and the inversion accuracy is greatly improved. And the numerical difference degree of the inversion result and the real model is compared, so that the full waveform inversion precision is improved by about 20%.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (6)
1. A roadway driving full-waveform inversion method based on multi-parameter constraint and structural correction is characterized by comprising the following steps:
step 1, inputting theoretical simulation data or actual measurement conventional roadway seismic advanced detection data, constructing a full-waveform inversion initial model, and performing single-scale inversion on the initial model by adopting a multi-scale elastic wave full-waveform inversion method;
step 2, carrying out multi-parameter weighted constraint structure correction on the single-scale inversion result obtained in the step 1 to obtain a primary correction result;
step 3, based on preset limiting conditions, carrying out one-dimensional wave velocity profile space structure correction and smooth constraint on the primary correction result to obtain a secondary correction result;
step 4, taking the secondary correction result obtained in the step 3 as an initial model of the next scale, and continuing to perform full waveform inversion according to the mode of the step 1;
and 5, repeatedly executing the step 2 to the step 4 until all scales are inverted, and obtaining a full waveform inversion result of the roadway advanced detection elastic wave.
2. The method as claimed in claim 1, wherein in step 2, the multiparameter weighted constraint structure correction is performed according to the following formula:
wherein, Δ m is the updating amount of the single iteration model; Δ m' is the model update amount after the multi-parameter weighting structure correction; m is an initial model; i and j respectively represent the positions of the grid nodes in the z and x directions; m is a unit of i,j p Parameter values of the positions of nodes of an i, j grid of a single parameter initial model of longitudinal wave velocity, transverse wave velocity or density are obtained; nz is the number of model vertical grid points, nx is the number of model horizontal grid points; v p Is the velocity of the longitudinal wave; v s Is the transverse wave velocity; den is the density.
3. The method for performing full waveform inversion along with the roadway based on multi-parameter constraint and structural correction as recited in claim 1, wherein in step 3, the preset limiting conditions are as follows:
setting a first limiting condition according to actual needs, wherein the first limiting condition takes the spatial correction lower limit value as a judgment basis, and the model update amount lower than the spatial correction lower limit value is set
Pressing to become 0;
and setting a second limiting condition according to actual needs, wherein the second limiting condition takes the range of the structural correction area and the distance between the areas as a judgment basis, when the model updating amount is a negative value area, both sides of the area have an updating amount which is a positive value area, or when the model updating amount is a positive value area, both sides of the area have an updating amount which is a negative value area, and when the area range and the distance between the areas meet the second limiting condition, the model updating amount in the positive value area existing on both sides of the negative value area or the negative value area existing on both sides of the positive value area is suppressed to be 1/5 of the original model updating value.
4. The method for full waveform inversion along with the roadway based on multi-parameter constraint and structure correction as recited in claim 2, wherein the step 3 comprises:
301, selecting a grid coordinate y of a longitudinal axis i Extracting one-dimensional wave velocity profile along the axis of the tunnel
Wherein x is a horizontal axis coordinate;
step 302, calculate a one-dimensional wave velocity profile
Corresponding model update amount
And slope
Step 303, updating the model below the value of the lower value of the spatial correction limit
Pressing to 0 to obtain corrected model update amount
Step 304, updating the quantity according to the corrected model
Judging and correcting each grid point on the two-dimensional wave velocity profile one by one according to the limiting conditions to obtain the corrected model updating amount
Step 305, updating quantity according to corrected model
And slope of
Judging each grid point on the one-dimensional wave velocity profile one by one according to the change condition of the slope symbol to obtain the corrected model updating amount
Step 306, updating the quantity according to the corrected model
Calculating the corrected one-dimensional wave velocity profile
Step 307, take down a grid coordinate y of the vertical axis i+1 And repeating the step 302 to the step 306 until the correction of all the grid one-dimensional wave velocity profiles is completed, and obtaining a secondary correction result after the spatial structure correction and the smoothing constraint of the one-dimensional wave velocity profiles.
5. The method of claim 4, wherein the step 305 comprises:
when two points with zero model updating quantity exist between three adjacent points with changed slope symbols, the symbols of the model updating quantity of grid points before and after the two points occurAnd if the two points are not provided with the point with the slope of 0, correcting the model updating amount of the grid point in the area between the three adjacent points with the changed slope sign to obtain the average value of the model updating amounts of the two points of the first slope sign change point and the third slope sign change point, and obtaining the corrected model updating amount
Then according to the updated quantity of the corrected model
Continuing to correct, when the signs of the model updating quantities of all the grid points in the range between two nonadjacent points with changed slope signs are not changed, correcting the model updating quantities of the grid points in the region between the two slope sign change points, changing the model updating quantities into the average value of the model updating quantities of all the grid points in the region, and obtaining the corrected model updating quantities
6. The method for performing full waveform inversion along with the roadway based on multi-parameter constraint and structural correction as claimed in claim 1, wherein the inversion result in step 5 comprises: longitudinal wave inversion results, shear wave inversion results and density inversion results.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211254264.0A CN115542393A (en) | 2022-10-13 | 2022-10-13 | Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction |
| PCT/CN2023/113748 WO2024078134A1 (en) | 2022-10-13 | 2023-08-18 | Excavation tunnel full-waveform inversion method based on multi-parameter constraint and structure correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211254264.0A CN115542393A (en) | 2022-10-13 | 2022-10-13 | Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115542393A true CN115542393A (en) | 2022-12-30 |
Family
ID=84733508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211254264.0A Pending CN115542393A (en) | 2022-10-13 | 2022-10-13 | Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN115542393A (en) |
| WO (1) | WO2024078134A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024078134A1 (en) * | 2022-10-13 | 2024-04-18 | 安徽理工大学 | Excavation tunnel full-waveform inversion method based on multi-parameter constraint and structure correction |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101318994B1 (en) * | 2012-02-24 | 2013-10-16 | 서울대학교산학협력단 | Method and apparatus of estimating underground structure using a plurality of weighting value |
| CN112698390B (en) * | 2020-11-11 | 2022-12-02 | 中国石油天然气股份有限公司 | Pre-stack seismic inversion method and device |
| CN114624779A (en) * | 2020-12-14 | 2022-06-14 | 中国石油化工股份有限公司 | Pre-stack multi-parameter inversion method for balanced model constraint |
| CN113552625B (en) * | 2021-06-21 | 2022-05-13 | 中国地质科学院地球物理地球化学勘查研究所 | Multi-scale full waveform inversion method for conventional land-domain seismic data |
| CN115542393A (en) * | 2022-10-13 | 2022-12-30 | 安徽理工大学 | Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction |
| CN115453624A (en) * | 2022-10-13 | 2022-12-09 | 安徽理工大学 | Mine following excavation seismic signal imaging method |
| CN115524744A (en) * | 2022-10-13 | 2022-12-27 | 安徽理工大学 | Roadway full-waveform inversion method based on cut-off frequency optimization and regularization constraints |
-
2022
- 2022-10-13 CN CN202211254264.0A patent/CN115542393A/en active Pending
-
2023
- 2023-08-18 WO PCT/CN2023/113748 patent/WO2024078134A1/en unknown
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024078134A1 (en) * | 2022-10-13 | 2024-04-18 | 安徽理工大学 | Excavation tunnel full-waveform inversion method based on multi-parameter constraint and structure correction |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2024078134A1 (en) | 2024-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10598819B2 (en) | Systems and methods to build sedimentary attributes | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| CN103913774B (en) | Reservoir geology mechanics parameter inversion method based on micro-seismic event | |
| CN115524744A (en) | Roadway full-waveform inversion method based on cut-off frequency optimization and regularization constraints | |
| CN112394404B (en) | Progressive reservoir fine characterization method | |
| CN103412333B (en) | A kind of static correction base level defining method | |
| EP2877880B1 (en) | System and method for migration velocity modeling | |
| CN105386756A (en) | Method for applying dependent variables to calculate porosity of brittle formation | |
| US20230384470A1 (en) | Method for three-dimensional velocity geological modeling with structures and velocities randomly arranged | |
| CN109632016A (en) | Rock And Soil adit digging and surrouding rock stress, strain monitoring experimental rig and its method | |
| CN103018778A (en) | Method and equipment for reservoir prediction by correcting acoustic logging curve | |
| CN112231974B (en) | Deep learning-based method and system for recovering seismic wave field characteristics of rock breaking seismic source of TBM (Tunnel boring machine) | |
| CN114035228A (en) | Tunnel seismic wave velocity inversion method and system based on deep learning | |
| CN115542393A (en) | Tunneling-following roadway full-waveform inversion method based on multi-parameter constraint and structural correction | |
| CN112364422A (en) | Shield construction earth surface deformation dynamic prediction method based on MIC-LSTM | |
| CN109459787A (en) | Underground coal mine structure imaging method and system based on SEISMIC CHANNEL WAVE METHOD full waveform inversion | |
| CN102590857A (en) | True surface relief prestack depth domain two-way wave imaging method | |
| CN111983683A (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| CN103217715B (en) | Multiple dimensioned regular grid Static Correction of Tomographic Inversion method | |
| CN111208568B (en) | Time domain multi-scale full waveform inversion method and system | |
| CN110907990B (en) | Quantitative prediction method and system for post-stack seismic cracks | |
| CN110687591B (en) | Method for determining physical property parameters of coal bed and surrounding rock based on waveform matching of prior data | |
| CN115453624A (en) | Mine following excavation seismic signal imaging method | |
| CN117270032A (en) | Sandstone uranium deposit uranium feature parameter inversion method and device based on three-dimensional earthquake | |
| CN116660974B (en) | Three-dimensional joint inversion method for bulk wave and surface wave based on structural coupling constraint |
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 |