CN114578414A - Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium - Google Patents
Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN114578414A CN114578414A CN202210065037.7A CN202210065037A CN114578414A CN 114578414 A CN114578414 A CN 114578414A CN 202210065037 A CN202210065037 A CN 202210065037A CN 114578414 A CN114578414 A CN 114578414A
- Authority
- CN
- China
- Prior art keywords
- candidate
- azimuth angle
- component
- azimuth
- seismic source
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000012806 monitoring device Methods 0.000 title description 5
- 230000001133 acceleration Effects 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000001615 p wave Methods 0.000 claims description 12
- 238000013519 translation Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 description 14
- 230000000875 corresponding effect Effects 0.000 description 13
- 230000010287 polarization Effects 0.000 description 12
- 238000001514 detection method Methods 0.000 description 11
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000000835 fiber Substances 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013589 supplement Substances 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
- G01V1/242—Seismographs
-
- 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
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a seismic source monitoring method, a device, computer equipment and a storage medium, wherein the method comprises the following steps: acquiring seismic data acquired by monitoring points in a preset historical time period; according to the first translational component, the second translational component and each candidate azimuth angle, obtaining a lateral acceleration corresponding to each candidate azimuth angle; determining a zero lag cross-correlation coefficient between each lateral acceleration and the third rotational component; and determining the final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross-correlation coefficient. The method can accurately determine the seismic source position of the earthquake.
Description
Technical Field
The invention relates to the technical field of seismic monitoring, in particular to a seismic source monitoring method, a seismic source monitoring device, computer equipment and a storage medium.
Background
When an earthquake with a large magnitude occurs, the fracture path and the velocity thereof are important factors for determining the destructive power thereof. And determining the kinematic characteristics of the fracture by using a finite fault inversion method based on the observed values of a large number of free parameters which are closely fitted. However, the creation of a kinematic model containing a priori information about the fault structure is characterized by intrinsic non-uniqueness and does not ensure mechanical consistency in terms of seismic dynamics. The ever-increasing computing resources allow observation-constrained dynamic fracturing models to be developed to supplement data-driven analysis. Such fracture scenarios provide a physically self-consistent description of complex fault slips, and their complexity limits the total number of possible numerical experiments. The rise of large-scale, dense seismic array instrumentation has made possible complementary techniques for tracking the seismic in space and time. This method requires very limited a priori knowledge to image the coherent high frequency energy radiation in a simple and fast manner.
Rotation can be derived from the spatial gradient of translational motion, and array observation is a measurement method that calculates the rotational component from the translational component. According to the method, a large number of translation sensors form a sensor array, and on a local uniform structure, a finite difference method is used for acting on translation components measured by the well-calibrated sensors, so that corresponding seismic rotation signals are obtained through calculation. This method requires the seismometer to be itself sensitive to rotational movement, and although the minimum number of stations required is three, it usually requires a few more stations with reasonably good azimuthal coverage to make stable observations. In addition, quality control of rotational component measurement based on array observation is very difficult. The relative spatial gradient at the reference station can be estimated by a certain number of dense arrays (OliveiraandBrolt, 1989; Spudichetal, 1995). FrohlichandPullam (1999) states that the classical single station method has some ambiguity, such as 180 ° back azimuth error, compared to the travel-time based method. Joint analysis of translational and rotational motion data may help overcome these disadvantages.
The most common techniques for imaging seismic properties using array data can be divided into two categories, both of which are based on analyzing the phase information of P-waves. The first type of method is based on conventional array measurements. Known as the backprojection method, the seismic energy radiation is imaged by applying an array beamforming technique. Backprojection was successfully demonstrated for the first time in 2004 by the sumatra-anderman seismic source (epicenter). Directional effects are used to characterize fault mechanisms (Kruger and Ohrnberger, 2005; Ishiietal, 2005). The second type of seismic fracture tracking method is to estimate the back azimuth (BAz) of the incident wave by a single station. In polarization analysis, the three translational components of a standard seismometer can be used to estimate the rear and incident angles of an incident wave. Coclardal (2006) proposes: on a free surface, the three rotational components in the x, y, z directions can be expressed as:
wherein
Representing the time inverse of the displacement wavefield. The rotation can therefore be derived from the spatial gradient of the transverse motion.
Bayer et al (2012) developed a single station method for tracking moving sources by polarization analysis of local and regional P-wave arrivals.
Disclosure of Invention
The application provides a seismic source monitoring method, a seismic source monitoring device, computer equipment and a storage medium.
A first aspect provides a seismic source monitoring method, the method comprising:
acquiring seismic data acquired by monitoring points in a preset historical time period, wherein the seismic data comprises: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis;
determining a candidate azimuth reversal angle from the seismic source to the monitoring point;
obtaining a lateral acceleration corresponding to each candidate anti-azimuth angle according to the first translational component, the second translational component and each candidate anti-azimuth angle;
determining a zero lag cross-correlation coefficient between each of the lateral acceleration and the third rotational component;
and determining the final azimuth reversal angle from the seismic source to the monitoring point according to the candidate azimuth reversal angle and the zero lag cross correlation coefficient.
In some embodiments, the determining a candidate azimuth back angle from the source to the monitor point comprises:
dividing the counter azimuth angles from 0 degree to 360 degrees into a plurality of parts at preset degree intervals to obtain a plurality of candidate counter azimuth angles.
In some embodiments, the determining a final azimuthing angle of the source to the monitor point based on the candidate azimuthing angle and the zero lag cross-correlation coefficient comprises:
selecting the candidate azimuths corresponding to the largest zero lag cross correlation coefficient as the final azimuths.
In some embodiments, the determining a candidate azimuth back angle for the source to the monitor point comprises:
obtaining an estimated negative azimuth angle according to the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude;
judging whether the estimated azimuth is less than 0, and if the estimated azimuth is not less than 0, taking the estimated azimuth as the candidate azimuth; and if the estimated azimuth angle is less than 0, adding 180 degrees to the estimated azimuth angle to obtain the candidate azimuth angle.
In some embodiments, the determining a candidate azimuth back angle from the source to the monitor point comprises:
judging whether the zero lag cross correlation coefficient is larger than 0, if so, adding 180 degrees to the candidate counter-azimuth angle to obtain the final counter-azimuth angle; and if the candidate azimuth angle is less than or equal to 0, the candidate azimuth angle is the final azimuth angle.
In some embodiments, the obtaining an estimated azimuth angle according to the first rotation component with the largest amplitude and the second rotation component with the largest amplitude includes:
substituting the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude into an estimated azimuth angle calculation formula to obtain the estimated azimuth angle, wherein the estimated azimuth angle calculation formula is as follows:
in the formula, thetaBAzTo estimate the azimuth-negation angle; reThe first rotation component with the maximum amplitude; rnThe second rotation component having the largest amplitude.
In some embodiments, the seismic data further comprises: the method comprises the following steps that the wave velocity of p waves and the wave velocity of s waves in seismic waves and the propagation time difference of the p waves and the s waves are obtained;
the method further comprises the following steps: and determining the distance from the seismic source of the seismic waves to a monitoring point according to the wave velocity of the p waves, the wave velocity of the s waves and the propagation time difference.
A second aspect provides a seismic source monitoring apparatus, the apparatus comprising:
the data input unit is used for acquiring seismic data acquired by monitoring points in a preset historical time period, and the seismic data comprises: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis;
the candidate determining unit is used for determining a candidate counter-azimuth angle from the seismic source to the monitoring point;
a lateral acceleration calculation unit, configured to obtain, according to the first translational component, the second translational component, and each candidate azimuths, a lateral acceleration corresponding to each candidate azimuths;
a cross-correlation coefficient calculation unit for determining a zero-lag cross-correlation coefficient between each of the lateral accelerations and the third rotation component;
and the result output unit is used for determining the final counter-azimuth angle from the seismic source to the monitoring point according to the candidate counter-azimuth angle and the zero lag cross-correlation coefficient.
A third aspect provides a computer apparatus comprising a memory and a processor, the memory having stored therein computer readable instructions which, when executed by the processor, cause the processor to perform the steps of the seismic source monitoring method described above.
A fourth aspect provides a storage medium having computer-readable instructions stored thereon which, when executed by one or more processors, cause the one or more processors to perform the steps of the seismic source monitoring method described above.
According to the seismic source monitoring method and device, the computer equipment and the storage medium, the translational components and the rotational motion components in three directions of the three-dimensional space acquired by the monitoring points are acquired; determining a reverse azimuth angle from a seismic source to a monitoring point according to the translation components and the rotation components in three directions of the three-dimensional space; and determining the distance from the seismic source to the monitoring point according to the wave velocity and the propagation time difference of the p wave and the s wave in the vibration wave. Therefore, the single method has no site limitation and the cost of manual arrangement is greatly reduced, and the counter-azimuth angle information can be acquired by combining the rotating component and the moving component of the single method. The cross-correlation method focuses on torsional motion (rotation of the vertical axis), and the polarization analysis method focuses on rocking motion (rotation of the horizontal axis). The estimation of the anti-azimuth angle based on the polarization analysis of the horizontal rotation component is not affected by a slight but significant asynchronization, since only one data logger is involved. Polarization analysis methods are more flexible to use and can be used in P, SV, SH, Rayleigh and Love waves.
Drawings
FIG. 1 is a flow chart of a seismic source monitoring method in one embodiment of the present application;
FIG. 2 is a flow chart of a seismic source monitoring method using a cross-correlation method to achieve source tracking according to an embodiment of the present application;
FIG. 3 is a waveform of a lateral acceleration and a third rotation component of a seismic source monitoring method according to an embodiment of the present application at a negative azimuth of 280; the abscissa corresponds to time (unit, second), and the ordinate corresponds to normalized amplitude; the solid line represents the third rotational component and the broken line represents the converted lateral acceleration;
FIG. 4 is a cross-correlation coefficient obtained by the cross-correlation method for each azimuth angle within a time window of 270s to 300s for the seismic source monitoring method according to an embodiment of the present application, where the abscissa corresponds to each azimuth angle and the ordinate is the corresponding cross-correlation coefficient;
FIG. 5 is a flow chart of a seismic source tracking implementation using a polarization analysis method of a seismic source monitoring method according to an embodiment of the present application;
FIG. 6 is a diagram illustrating the calculated azimuth and theoretical azimuth angles predicted by a polarization analysis method over a time window for a seismic source monitoring method according to an embodiment of the present application; wherein the black asterisk is the predicted azimuth of each time window, and the horizontal solid line is the theoretical numerical value of the azimuth; the abscissa is the sampling time (unit of second), and the ordinate is the angle of the negative azimuth;
FIG. 7 is a block diagram of a seismic source monitoring device according to an embodiment of the present application;
FIG. 8 is a block diagram showing an internal configuration of a computer device according to one embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
It will be understood that the terms "first," second, "third," and the like as used herein may be used herein to describe various elements, but these elements are not limited by these terms. These terms are only used to distinguish one element from another.
An embodiment of a seismic source monitoring method is provided, in which a computer device and a seismic monitor may be included.
The computer device has an interface, which may be, for example, an API (application programming interface). The seismic monitor is a device for collecting seismic data arranged at a detection point, and the computer equipment carries out seismic source monitoring (such as determining the counter-azimuth angle of a seismic source) according to the data collected by the seismic monitor.
It should be noted that the terminal and the computer device may be a smart phone, a tablet computer, a notebook computer, a desktop computer, and the like, but are not limited thereto. The computer device and the terminal may be connected through bluetooth, USB (universal serial bus), or other communication connection methods, which is not limited herein.
As shown in fig. 1 to 6, in one embodiment, a seismic source monitoring method is proposed, which may be applied in the computer device described above, and which may be implemented based on a python program. As shown in fig. 1, the method may specifically include the following steps:
it can be understood that, in this step, obtaining the translation components and rotation components of three coordinate axes of the stereo space acquired by the monitoring point includes:
adopting a six-component seismograph to acquire translation components and rotation components in three directions of a three-dimensional space; the six-component seismograph may include: the system comprises three fiber-optic gyroscopes, three accelerometers and a six-component signal processing module, wherein the sensitive axes of the three fiber-optic gyroscopes are mutually orthogonal, and the sensitive axes of the three accelerometers are mutually orthogonal; the sensitive axis of each fiber optic gyroscope is parallel to or coincided with the three sensitive axes of one accelerometer one by one; the input end of the six-component signal processing module is connected with the detection signal output end of each fiber-optic gyroscope and the detection signal output end of each accelerometer, the six-component signal processing module is used for generating modulation signals required by the fiber-optic gyroscopes and obtaining detection angular velocity according to detection signals output by the detection signal output ends of the fiber-optic gyroscopes, and the six-component signal processing module is used for performing error compensation on the detection angular velocity and obtaining detection translational acceleration according to detection signals output by the detection signal output ends of the three-axis accelerometers and performing error compensation on the detection translational acceleration. The method comprises the steps of establishing a space rectangular coordinate system, wherein the space rectangular coordinate system comprises a first coordinate axis (X axis), a second coordinate axis (Y axis) and a third coordinate axis (Z axis) which are mutually orthogonal, and sensitive axes of the fiber-optic gyroscope and the accelerometer are respectively arranged on the first coordinate axis (X axis), the second coordinate axis (Y axis) and the third coordinate axis (Z axis).
Furthermore, the six-component seismograph comprises three orthogonal fiber optic gyroscopes and three same orthogonal accelerometers, the fiber optic gyroscopes are used for measuring angular velocity information of the object, the three-axis accelerometers are used for measuring line motion information of the carrying body, and the six-component signal processing module can provide complete seismic wave field information by processing the outputs of the fiber optic gyroscopes and the accelerometers.
it will be appreciated that the candidate azimuths may be determined empirically or from seismic data.
In some embodiments, the determination of the final azimuth opposition angle of the seismic source is implemented by a cross-correlation method, and the empirically setting the candidate azimuth opposition angle in step 102 of the cross-correlation method may include:
step 1021a, dividing the counter azimuth angles from 0 degree to 360 degrees into a plurality of parts at preset degree intervals to obtain a plurality of candidate counter azimuth angles.
In some embodiments, the determining the final azimuth reversal angle of the seismic source is implemented by a polarization analysis method, and the determining the candidate azimuth reversal angle according to the seismic data in step 102 in the polarization analysis method may include:
step 1022a, obtaining an estimated negative azimuth angle according to the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude;
in this step, substituting the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude into an estimated azimuth angle calculation formula to obtain an estimated azimuth angle, wherein the estimated azimuth angle calculation formula is as follows:
in the formula, thetaBAzTo estimate the azimuth-negation angle; reThe first rotation component with the maximum amplitude; rnThe second rotation component having the largest amplitude.
Step 1022b, determining whether the estimated azimuth angle is less than 0, and if the estimated azimuth angle is not less than 0, taking the estimated azimuth angle as a candidate azimuth angle; and if the estimated azimuth angle is less than 0, adding 180 degrees to the estimated azimuth angle to obtain a candidate azimuth angle.
103, obtaining a lateral acceleration corresponding to each candidate azimuth according to the first translational component, the second translational component and each candidate azimuth;
in the step, the first translational component, the second translational component and each candidate azimuth angle are substituted into a lateral acceleration formula to obtain the lateral acceleration of each candidate azimuth angle.
Wherein, the lateral acceleration formula is:
in the formula, AeIs a first translational component; a. thenIs the second translational component; thetaBazIs a candidate anti-azimuth angle; a. thetIs the lateral acceleration.
104, determining a zero lag cross-correlation coefficient between each transverse acceleration and the third rotation component;
in the step, each transverse acceleration and the third rotation component are substituted into a correlation calculation formula to obtain a zero-lag cross correlation coefficient between the transverse acceleration and the third rotation component. Wherein the correlation calculation formula is
In the formula, AtThe transverse acceleration corresponding to each acquisition moment; r iszVertical rotation components collected for each collection time;
is the average of the lateral acceleration;
is the average of the vertical rotational components; xcorr (A)t,Rz) Is the zero lag cross-correlation coefficient between the lateral acceleration and the vertical rotational component.
It will be appreciated that the six-component seismometer acquires seismic data at predetermined time intervals (e.g., once every 1 second), with a corresponding lateral acceleration (a) at each acquisition timet) The method is obtained by substituting a first translational component, a second translational component and a counter-azimuth angle (candidate counter-azimuth angle or final counter-azimuth angle) acquired at each acquisition moment into a transverse acceleration calculation formula. Average value of lateral acceleration
Is the average of the lateral acceleration over a period of time (a window of time, e.g., 3 seconds).
And 105, determining a final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross-correlation coefficient.
In this step, there is a certain relation between the zero lag cross-correlation coefficient between the lateral acceleration and the third rotational component and the anti-square angle of the seismic source, so that the final anti-azimuth angle from the seismic source to the monitoring point can be determined according to the candidate anti-azimuth angle and the zero lag cross-correlation coefficient.
In some embodiments, the determining the final azimuth-reversal angle of the seismic source is implemented by using a cross-correlation method, and then in the step 105 of the cross-correlation method, determining the final azimuth-reversal angle from the seismic source to the monitoring point according to the candidate azimuth-reversal angle and the zero-lag cross-correlation coefficient may include:
step 1051a, selecting the candidate azimuth angle corresponding to the largest zero lag cross-correlation coefficient as the final azimuth angle.
In some embodiments, the determining the final azimuth reversal angle of the seismic source is implemented by using a polarization analysis method, and then in step 105 of the polarization analysis method, the determining the final azimuth reversal angle from the seismic source to the monitoring point according to the candidate azimuth reversal angle and the zero-lag cross-correlation coefficient may include:
step 1052a, judging whether the zero lag cross correlation coefficient is greater than 0, and if so, adding 180 degrees to the candidate azimuth angle to obtain a final azimuth angle; if the azimuth angle is less than or equal to 0, the candidate azimuth angle is the final azimuth angle.
Further, it may be understood that the seismic data acquired by the monitoring point within the preset historical time period may be acquired within 12s of the past time of the current time, the seismic data is acquired at a frequency of once every 1 second by the monitoring point, a final azimuth reversal angle from the seismic source to the monitoring point is determined according to the seismic data, the seismic data may be divided into 3 time windows (each time window has a duration of 4s), one time window is taken as a minimum unit, a candidate azimuth reversal angle in each time window is determined according to the seismic data in each time window in the above embodiment, a final azimuth reversal angle in each time window is further determined, and finally, a azimuth reversal angle from the seismic source to the monitoring point is determined according to the final azimuth reversal angles of the time windows. That is, the predicted final azimuth angles are mostly distributed near the theoretical azimuth angle, the result error of individual time windows is large, and may be influenced by other external factors, and the overall predicted result is matched with the theoretical value, so that the final azimuth angles of a plurality of time windows are fitted in a linear manner to obtain the azimuth angles from the seismic source to the monitoring point. Of course, after dividing the time window, "maximum amplitude", "maximum zero lag cross-correlation coefficient", and the like in the above embodiments are maximum within the unit time window.
As shown in fig. 2 to 4, in an application scenario, the feasibility of the cross-correlation method was tested by using 6.4-level seismic data of yunnan university collected by wuhan instruments in 5/21/2021. The theoretical azimuth angle is 280 degrees, and the distance between a seismic source (epicenter) and a test station is about 1583 km.
As shown in fig. 2, the method may include:
and step A1, acquiring three translational components and rotational component seismic data from a rotational seismograph, wherein the length of the intercepted and recorded data is 390s, and the method starts from the occurrence time of the earthquake. Then, preprocessing such as instrument response removal is carried out on the data;
step a2, dividing the counter-azimuth angles from 0 ° to 360 ° into 72 parts at 5 ° intervals, and for each given counter-azimuth angle, obtaining the lateral acceleration by using the first horizontal component and the second horizontal component according to a lateral acceleration formula. FIG. 3 is a waveform of lateral acceleration and a third rotation component at a negative azimuth angle of 280 °, and it can be seen that the waveform can be better fitted;
step A3, dividing recorded 390s seismic data into 13 parts uniformly according to 30s as a time window. Within each time window, all candidate azimuths and each time window are cycled. The obtained transverse acceleration AtWith a third rotational component RzSubstituting into the correlation calculation formula, a zero lag cross correlation coefficient between the lateral acceleration and the third rotational component. FIG. 4 is a cross-correlation coefficient for each of the azimuths within a 270s to 300s time window;
and step A4, selecting the azimuth angle corresponding to the maximum value position of the cross-correlation coefficient as the final result of the time window.
And after all the time window cycles are finished, the counter azimuth angle prediction result of the whole time channel can be obtained.
In an application scenario, the feasibility of the polarization analysis method was tested using a methane explosion experiment conducted in Guangzhou at 12.1.1.2021, as shown in FIGS. 5 and 6. The theoretical azimuth reversal angle is 190 degrees, and the distance between the methane source and the test station is about 305 m. As shown in fig. 5, the method may include:
and step B1, downloading from the seismograph to obtain three translational component and rotational component seismic data, wherein the length of the intercepted and recorded data is 15s, and the method starts from the explosion of the methane source. Then, preprocessing such as instrument response removal and the like is carried out on the downloaded data;
and step B2, uniformly dividing the recorded 15s seismic data into 15 parts according to 1s as a time window. Within each time window, the first and second rotational components of maximum magnitude are selected to compute candidate azimuths.
And step B3, converting the horizontal component into the lateral acceleration according to the calculated candidate azimuth angle in each 1s time window.
And step B4, substituting the first translational component, the second translational component and the candidate azimuth angles into a transverse acceleration formula to obtain the transverse acceleration of each candidate azimuth angle.
And step B5, if the transverse acceleration is positively correlated with the third rotation component, adding 180 degrees to the obtained counter-azimuth angle, and otherwise, keeping the counter-azimuth angle unchanged.
Step B6, looping each time window to obtain the prediction result of the azimuth reversal angle of the whole time channel, as shown in fig. 6. It can be seen that the predicted azimuth reversals are relatively uniformly distributed near the theoretical azimuth, but some time windows are influenced by other external factors, so that the result error is large.
In some embodiments, the wave velocity of p-waves and the wave velocity of s-waves and the propagation time difference of p-waves and s-waves in the seismic waves;
the method further comprises the following steps: and determining the distance from the seismic source of the seismic waves to the monitoring point according to the wave velocity of the p waves, the wave velocity of the s waves and the propagation time difference. Specifically, the method comprises the following steps of calculating according to the following formula:
wherein S is the distance from the seismic source to the monitoring point; v. ofpIs the wave velocity of the p-wave; and vsRespectively the wave velocity of the s-wave; t is t1And t2Respectively, the arrival time of the p-wave and the s-wave obtained from the waveform.
In one embodiment, there is provided an apparatus for determining seismic wave phase velocity, which may be integrated in the computer device described above, as shown in fig. 7, and may specifically include:
the data input unit 711 is configured to acquire seismic data acquired by a monitoring point in a preset historical time period, where the seismic data includes: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis;
a candidate determination unit 712 for determining candidate azimuths from the source to the monitoring points;
the transverse acceleration calculation unit is used for obtaining the transverse acceleration corresponding to each candidate anti-azimuth angle according to the first translational component, the second translational component and each candidate anti-azimuth angle;
a cross-correlation coefficient calculation unit 713 for determining a zero-lag cross-correlation coefficient between each lateral acceleration and the third rotation component;
and a result output unit 714, configured to determine a final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross-correlation coefficient.
As shown in fig. 8, in one embodiment, a computer device is presented that may include a processor, a storage medium, a memory, and a network API interface connected by a system bus. The computer device comprises a storage medium, an operating system, a database and computer readable instructions, wherein the database can store control information sequences, and the computer readable instructions can enable a processor to realize the seismic source monitoring method when being executed by the processor. The processor of the computer device is used for providing calculation and control capability and supporting the operation of the whole computer device. The memory of the computer device may have stored therein computer readable instructions that, when executed by the processor, may cause the processor to perform a seismic source monitoring method. The network API interface of the computer device is used for connecting and communicating with the terminal. Those skilled in the art will appreciate that the architecture shown in fig. 8 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
The processor, when executing the computer program, implements the steps of: acquiring seismic data acquired by monitoring points in a preset historical time period, wherein the seismic data comprises: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis; determining candidate counter-azimuth angles from the seismic source to the monitoring points; according to the first translational component, the second translational component and each candidate azimuth angle, obtaining a lateral acceleration corresponding to each candidate azimuth angle; determining a zero lag cross-correlation coefficient between each lateral acceleration and the third rotational component; and determining the final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross-correlation coefficient.
In one embodiment, a storage medium is presented having computer-readable instructions stored thereon which, when executed by one or more processors, cause the one or more processors to perform the steps of: acquiring seismic data acquired by monitoring points in a preset historical time period, wherein the seismic data comprises: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis; determining candidate counter-azimuth angles from the seismic source to the monitoring points; according to the first translational component, the second translational component and each candidate azimuth angle, obtaining a lateral acceleration corresponding to each candidate azimuth angle; determining a zero lag cross-correlation coefficient between each lateral acceleration and the third rotational component; and determining the final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross-correlation coefficient.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and can include the processes of the embodiments of the methods described above when the computer program is executed. The storage medium may be a non-volatile storage medium such as a magnetic disk, an optical disk, a Read-only memory (ROM), or a Random Access Memory (RAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above examples only show some embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent should be subject to the appended claims.
Claims (10)
1. A seismic source monitoring method, the method comprising:
acquiring seismic data acquired by monitoring points in a preset historical time period, wherein the seismic data comprises: a first translation component and a first rotation component on a first coordinate axis, a second translation component and a second rotation component on a second coordinate axis, and a third translation component and a third rotation component on a third coordinate axis, which are perpendicular to each other in a three-dimensional space;
determining candidate counter-azimuth angles from the seismic source to the monitoring points;
obtaining a lateral acceleration corresponding to each candidate anti-azimuth angle according to the first translational component, the second translational component and each candidate anti-azimuth angle;
determining a zero lag cross-correlation coefficient between each of the lateral acceleration and the third rotational component;
and determining the final azimuth angle from the seismic source to the monitoring point according to the candidate azimuth angle and the zero lag cross correlation coefficient.
2. The seismic source monitoring method of claim 1, wherein the determining candidate azimuths for sources to the monitoring points comprises:
dividing the counter azimuth angles from 0 degree to 360 degrees into a plurality of parts at preset degree intervals to obtain a plurality of candidate counter azimuth angles.
3. The seismic source monitoring method of claim 2, wherein the determining a final azimuthing angle of the source to the monitoring point based on the candidate azimuthing angle and the zero lag cross-correlation coefficient comprises:
selecting the candidate azimuths corresponding to the largest zero lag cross correlation coefficient as the final azimuths.
4. The seismic source monitoring method of claim 2, wherein the determining candidate azimuths for sources to the monitoring points comprises:
obtaining an estimated negative azimuth angle according to the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude;
judging whether the estimated azimuth is less than 0, and if the estimated azimuth is not less than 0, taking the estimated azimuth as the candidate azimuth; and if the estimated azimuth angle is less than 0, adding 180 degrees to the estimated azimuth angle to obtain the candidate azimuth angle.
5. The seismic source monitoring method of claim 4, wherein the determining a final azimuthing angle of the source to the monitoring point based on the candidate azimuthing angle and the zero lag cross-correlation coefficient comprises:
judging whether the zero lag cross correlation coefficient is larger than 0, if so, adding 180 degrees to the candidate counter-azimuth angle to obtain the final counter-azimuth angle; and if the candidate azimuth angle is less than or equal to 0, the candidate azimuth angle is the final azimuth angle.
6. The seismic source monitoring method of claim 4, wherein obtaining the estimated azimuths from the first rotating components with the largest amplitudes and the second rotating components with the largest amplitudes comprises:
substituting the first rotation component with the maximum amplitude and the second rotation component with the maximum amplitude into an estimated azimuth angle calculation formula to obtain the estimated azimuth angle, wherein the estimated azimuth angle calculation formula is as follows:
in the formula, thetaBAzTo estimate the azimuth of the reaction; r iseThe first rotation component with the maximum amplitude; r isnThe second rotation component with the largest amplitude.
7. The seismic source monitoring method of claim 1, wherein the seismic data further comprises: the method comprises the following steps that the wave velocity of p waves and the wave velocity of s waves in seismic waves and the propagation time difference of the p waves and the s waves are obtained;
the method further comprises the following steps: and determining the distance from the seismic source of the seismic waves to a monitoring point according to the wave velocity of the p waves, the wave velocity of the s waves and the propagation time difference.
8. A seismic source monitoring apparatus, the apparatus comprising:
the data input unit is used for acquiring seismic data acquired by monitoring points in a preset historical time period, and the seismic data comprises: a first translational component and a first rotational component on a first coordinate axis which are perpendicular to each other in a three-dimensional space, a second translational component and a second rotational component on a second coordinate axis, and a third translational component and a third rotational component on a third coordinate axis;
the candidate determining unit is used for determining a candidate counter-azimuth angle from the seismic source to the monitoring point;
a lateral acceleration calculation unit, configured to obtain, according to the first translational component, the second translational component, and each candidate azimuths, a lateral acceleration corresponding to each candidate azimuths;
a cross-correlation coefficient calculation unit for determining a zero-lag cross-correlation coefficient between each of the lateral accelerations and the third rotation component;
and the result output unit is used for determining the final counter-azimuth angle from the seismic source to the monitoring point according to the candidate counter-azimuth angle and the zero lag cross-correlation coefficient.
9. A computer apparatus comprising a memory and a processor, the memory having stored therein computer readable instructions which, when executed by the processor, cause the processor to perform the steps of the seismic source monitoring method of any of claims 1 to 7.
10. A storage medium having computer readable instructions stored thereon which, when executed by one or more processors, cause the one or more processors to perform the steps of the seismic source monitoring method as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210065037.7A CN114578414A (en) | 2022-01-20 | 2022-01-20 | Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210065037.7A CN114578414A (en) | 2022-01-20 | 2022-01-20 | Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114578414A true CN114578414A (en) | 2022-06-03 |
Family
ID=81772422
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210065037.7A Pending CN114578414A (en) | 2022-01-20 | 2022-01-20 | Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114578414A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2750253A1 (en) * | 2009-01-20 | 2010-07-29 | Spectraseis Ag | Time reverse imaging operators for source location |
| CN112051606A (en) * | 2020-09-10 | 2020-12-08 | 北京大学 | Six-component seismograph |
| CN113608262A (en) * | 2021-08-11 | 2021-11-05 | 中国地质大学（北京） | Seismic data processing method and device for calculating rotation component by using translation component |
-
2022
- 2022-01-20 CN CN202210065037.7A patent/CN114578414A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2750253A1 (en) * | 2009-01-20 | 2010-07-29 | Spectraseis Ag | Time reverse imaging operators for source location |
| CN112051606A (en) * | 2020-09-10 | 2020-12-08 | 北京大学 | Six-component seismograph |
| CN113608262A (en) * | 2021-08-11 | 2021-11-05 | 中国地质大学（北京） | Seismic data processing method and device for calculating rotation component by using translation component |
Non-Patent Citations (1)
| Title |
|---|
| 许川佩;梁光发;吴玉龙;: "一种地震预警监测传感器模块的设计及实现", 地震工程学报, no. 03, 30 September 2013 (2013-09-30) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cochard et al. | Rotational motions in seismology: theory, observation, simulation | |
| Yuan et al. | Six degree‐of‐freedom broadband ground‐motion observations with portable sensors: Validation, local earthquakes, and signal processing | |
| Geng et al. | A new seismogeodetic approach applied to GPS and accelerometer observations of the 2012 Brawley seismic swarm: Implications for earthquake early warning | |
| RU2494418C1 (en) | Method of measuring coordinates of microseismic sources and parameters of mechanisms of centres thereof in conditions of strong seismic interference (versions) | |
| Donner et al. | Comparing direct observation of strain, rotation, and displacement with array estimates at Piñon Flat Observatory, California | |
| CN103582826A (en) | Determining an indication of wavefield velocity | |
| RU2616650C2 (en) | Method and device of seismic data processing | |
| Barak et al. | Six-component seismic land data acquired with geophones and rotation sensors: Wave-mode selectivity by application of multicomponent polarization filtering | |
| US20130291637A1 (en) | System and Method For Monitoring Mechanically Coupled Structures | |
| EP2951618B1 (en) | Computing rotation data using a gradient of translational data | |
| CN107918156A (en) | Detect the method and device of subsea node collection earthquake data polarity | |
| TW201428297A (en) | Angular velocity estimation using a magnetometer and accelerometer | |
| US20180059278A1 (en) | Seismic sensor devices, systems, and methods including noise filtering | |
| Ross et al. | Low‐frequency tilt seismology with a precision ground‐rotation sensor | |
| Schneider et al. | Inversion of geo-magnetic full-tensor gradiometer data | |
| Chen et al. | Six‐Component Earthquake Synchronous Observations Across Taiwan Strait: Phase Velocity and Source Location | |
| Noe et al. | Anisotropic elastic parameter estimation from multicomponent ground-motion observations: a theoretical study | |
| Kovalyshen et al. | Inversion of ultrasonic data for transversely isotropic media | |
| Giammarinaro et al. | Seismic surface wave focal spot imaging: numerical resolution experiments | |
| US7388538B1 (en) | System and method for obtaining attitude from known sources of energy and angle measurements | |
| CN114578414A (en) | Earthquake focus monitoring method, earthquake focus monitoring device, computer equipment and storage medium | |
| CN112649889A (en) | Six-component seismic data and absolute gravity measuring instrument and measuring method | |
| US9354339B2 (en) | Systems and methods for combining triaxial geophone data for localizing nearby transient seismic sources | |
| Poppeliers et al. | Three‐dimensional wave gradiometry for polarized seismic waves | |
| Barak et al. | Recording active-seismic ground rotations using induction-coil magnetometers |
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 |