US20230324577A1 - Real-Time Microseismic Magnitude Calculation Method and Device Based on Deep Learning - Google Patents
Real-Time Microseismic Magnitude Calculation Method and Device Based on Deep Learning Download PDFInfo
- Publication number
- US20230324577A1 US20230324577A1 US18/304,026 US202318304026A US2023324577A1 US 20230324577 A1 US20230324577 A1 US 20230324577A1 US 202318304026 A US202318304026 A US 202318304026A US 2023324577 A1 US2023324577 A1 US 2023324577A1
- Authority
- US
- United States
- Prior art keywords
- microseismic
- magnitude
- data
- real
- deep learning
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 82
- 238000013135 deep learning Methods 0.000 title claims abstract description 49
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000012544 monitoring process Methods 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 230000004927 fusion Effects 0.000 claims abstract description 5
- 239000002184 metal Substances 0.000 claims description 21
- 238000009826 distribution Methods 0.000 claims description 19
- 230000003287 optical effect Effects 0.000 claims description 18
- 239000013307 optical fiber Substances 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 12
- 239000004568 cement Substances 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 9
- 239000002002 slurry Substances 0.000 claims description 9
- 230000007246 mechanism Effects 0.000 claims description 8
- 238000011176 pooling Methods 0.000 claims description 8
- 230000004913 activation Effects 0.000 claims description 7
- 238000003062 neural network model Methods 0.000 claims description 7
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000005086 pumping Methods 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000013527 convolutional neural network Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013434 data augmentation Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- AKPLHCDWDRPJGD-UHFFFAOYSA-N nordazepam Chemical compound C12=CC(Cl)=CC=C2NC(=O)CN=C1C1=CC=CC=C1 AKPLHCDWDRPJGD-UHFFFAOYSA-N 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0464—Convolutional networks [CNN, ConvNet]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/047—Probabilistic or stochastic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/09—Supervised learning
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Definitions
- the present disclosure relates to the field of microseismic monitoring and deep learning.
- the disclosure relates to a real-time microseismic magnitude calculation method based on deep learning and a corresponding device, wherein the method includes using a training data set to train a magnitude calculation module, analyzing and processing field data, and inputting the processed data to the magnitude calculation module so as to obtain the calculated magnitude, etc. Therefore, the present disclosure belongs to the field of seismic data processing in seismic monitoring.
- Hydraulic fracturing technology relates to injection of high-pressure fluid into shale reservoirs to generate complex artificial fractures, which can increase the connectivity of reservoirs and increase well production. Monitoring and evaluating the effects of hydraulic fracturing on reservoir substructures in different stages is a prerequisite for efficient exploitation and safe production. Microseismic monitoring technology is an important means for real-time monitoring of hydraulic fracturing stimulation of unconventional resources. Microseismic monitoring technology evaluates fracturing effects by monitoring microseismic signals generated during hydraulic fracturing, and then improves and optimizes engineering parameters.
- Magnitude estimation results often lead to overestimation of small earthquakes and underestimation of larger earthquakes, to a greater or lesser extent. Further, we can try to estimate the magnitude by using two or more of the above three types of parameters in combination, so as to improve the accuracy of earthquake magnitude estimation. Magnitude estimation methods also include methods of characteristic frequency parameters based on period, such as the maximum predominant period Pd (peak displacement) method, the average predominant period ⁇ c method, and so on. In addition, if there are data from multiple stations, the magnitudes from each station can be estimated, and then the weighted average can be used to obtain a magnitude estimate.
- Pd peak displacement
- a recently proposed method uses strong motion data to establish a NN network for earthquake pre-warning magnitude estimation, using waveform spectrum information as model input. It was found that, within the 3 s time period, this method is superior to the traditional single parameter method, and it is also believed that the use of multi-characteristic combination input will lead to a better magnitude estimation model. Mousavi & Beroza (2020) constructed a magnitude estimation model composed of a convolutional neural network and a recurrent neural network based on a 30 seconds of single-station waveform input, which effectively improved the magnitude estimation results.
- a convolutional neural network model with multiple fully connected layers is established, and the frequency domain data of P wave at different periods of 3-9 s, the source information (epicentral distance and source depth) and site information of the corresponding earthquake event are selected as full data set, to train the proposed model and evaluate the estimation performance.
- Lomax (2019) used three-component acceleration time history data of a single station as input data of a first convolutional layer, normalized the acceleration peak value in the three-component band, and used the normalized value as auxiliary input parameters of the fully connected layer, then using the above two parts of data to iteratively train the model to achieve the purpose of automatic detection and rapid characterization of earthquake magnitude.
- the purpose of the present disclosure is to at least partly overcome the defects of the prior art, and provide a real-time microseismic magnitude calculation method and a corresponding device based on deep learning.
- the present disclosure provides a real-time microseismic magnitude calculation method based on deep learning, in which the strain data collected by a DAS system may be inputted directly without converting the DAS strain data to conventional data (generally displacement, velocity or acceleration data), thereby improving the ability to quickly estimate the microseismic magnitude.
- the present disclosure provides a real-time microseismic magnitude calculation method based on deep learning to improve the accuracy of the microseismic magnitude estimation.
- the present disclosure provides a real-time microseismic magnitude calculation device based on deep learning, which is used for implementing the aforementioned real-time microseismic magnitude calculation method based on deep learning.
- a real-time microseismic magnitude calculation method based on deep learning comprising:
- the training data set comprises Nk DAS microseismic strain data
- each DAS microseismic strain data comprises Nt channels
- each channel comprises Ns sampling points, forming a data body of Nk*Nt*Ns, where Nk, Nt and Ns are all 1000.
- the theoretical magnitude Gaussian distribution label is a Gaussian distribution centered on the theoretical magnitude, the magnitude ranges from -3 to 3, the resolution is 0.1, and it is a vector with 61 values.
- the constructing a training data set comprises:
- the single initial microseismic data is generated by:
- the magnitude distribution is a Gaussian distribution of the current microseismic theoretical magnitude, with a distribution calculation formula of
- z 0 is a theoretical magnitude
- ⁇ is a Gaussian distribution value
- the magnitude ranges from -3 to 3
- the resolution is 0.1
- the number of type labels is 61
- the displacement M(z) corresponds to the current magnitude.
- the data-augmenting the initial microseismic data comprises:
- the preprocessing all data comprises:
- the magnitude calculation module uses a 2D convolutional layer
- the convolutional layer group comprises a convolutional layer, a max pooling layer, and a Dropout layer
- the convolution operation uses a ReLu activation function to extract characteristic from the inputted data
- the max pooling layer is after the convolutional layer
- a rectified linear unit is used for activation in each hidden convolutional layer
- the kernel size of the convolution is set to 64 ⁇ 3 ⁇ 3
- the final output layer is a vector of 1*61, corresponding to magnitude -3 to magnitude 3 of seism magnitude, then a maximum value is taken as the current calculated magnitude.
- the 3-layer convolutional structure comprises a convolutional layer, a pooling layer and a fully connected layer.
- the training the magnitude calculation module by using the training data set comprises:
- the DAS-based horizontal well microseismic monitoring system comprises a metal sleeve and an armored optical cable fixedly arranged outside the metal sleeve, and a single-mode optical fiber is arranged in the armored optical cable;
- constructing a DAS-based horizontal well microseismic monitoring system comprises:
- the preprocessing the field data comprises:
- a real-time microseismic magnitude calculation device based on deep learning wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to any one of the foregoing embodiments.
- the present disclosure proposes a real-time microseismic magnitude calculation method and a corresponding device based on deep learning, which integrates waveform and frequency spectrum characteristics with a deep model network, and trains the model through continuous DAS strain records, without the need for converting the DAS strain data to conventional seismic data (generally displacement, velocity or acceleration data), improving the ability to quickly estimate the microseismic magnitude and improving the accuracy of the microseismic magnitude estimation.
- FIG. 1 is a flowchart of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure
- FIG. 2 shows a software data processing flow of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure
- FIG. 3 is a schematic diagram of a network structure of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure.
- DAS Distributed Acoustic Sensing
- DAS is an emerging data acquisition technology that has developed rapidly in recent years.
- DAS generally consists of a demodulator (including an optical system and a signal acquisition system) and a normal optical fiber or a special optical fiber for sensing.
- DAS realizes the measurement of strain signal along the axial direction of the optical fiber by detecting a phase change of backward Rayleigh scattering light generated by the laser pulse in the internal scatterer of the optical fiber.
- DAS optical fibers are generally arranged along the horizontal well. Compared with conventional geophone detection, DAS can obtain the measurement data of the whole well section, which greatly increases the azimuth of data collection; DAS has 1 m spatial sampling, which significantly improves the density and spatial continuity of the monitoring data.
- the deep learning method is a good tool to solve difficult problems in the field of seismology and geophysical research.
- the application of deep learning in the field of seismology continues to expand, for example, to seismic phase identification, first arrival picking, source location, etc.
- the present disclosure proposes a method for quickly estimating seismic magnitude based on deep learning and feature fusion, which integrates waveform and frequency spectrum characteristics with deep model network, and trains the model through continuous DAS strain records, improving the ability to quickly estimate the seismic magnitude.
- a real-time microseismic magnitude calculation method based on deep learning is proposed, and the calculation method is carried out in accordance with FIG. 1 ( FIG. 1 ):
- a DAS-based horizontal well microseismic monitoring system is constructed.
- the DAS-based horizontal well microseismic monitoring system comprises a metal sleeve and an armored optical cable fixedly arranged outside the metal sleeve, and a single-mode optical fiber is arranged in the armored optical cable; the DAS modem is placed near a wellhead, the DAS signal port of the DAS modem is connected to the special single-mode optical fiber outside the metal sleeve.
- a DAS-based horizontal well microseismic monitoring system includes:
- the DAS-based horizontal well microseismic monitoring system may use the armored optical cables permanently arranged outside the metal sleeve to perform hydraulic fracturing microseismic monitoring.
- the single-mode optical fibers outside the metal sleeve arranged in the well are used to collect data, the data is transmitted to the DAS modem near the wellhead for demodulation, then the microseismic events are obtained, wherein the microseismic events are generated from the stimulated fracture of the underground reservoir of the offset well or in-well when the continuously- hydraulic fracturing operation is performed.
- relevant information of the microseismic events can be obtained, and such information includes the time, three-dimensional spatial position and energy magnitude of the microseismic events.
- focal mechanism analysis and magnitude analysis can be performed according to the recorded longitudinal wave signal characteristics of the microseismic events, so as to obtain fracture mechanism of the most microseismic events, and the distribution characteristic and regularity of fractures after the hydraulic fracturing operation may be analyzed. Based on the above information, the effective and reliable, qualitative and quantitative evaluation of the hydraulic fracturing effect of reservoir can be carried out.
- the real-time microseismic magnitude calculation method based on deep learning further includes:
- the structure of the data set is like this: (2D microseismic data, 1D theoretical magnitude Gaussian distribution label).
- 2D microseismic data a total of Nt (1000) channels are included, each channel includes Ns (1000) sampling points, and the training data set includes Nk (1000) DAS microseismic strain data.
- the theoretical magnitude Gaussian distribution label is a 1D Gaussian distribution centered on the theoretical magnitude. The magnitude ranges from -3 to 3 Ml, and the resolution is 0.1, and it is a vector with 61 values. Therefore, the data body is Nk*Nt*Ns and the number of the type labels is Nk*61.
- the magnitude distribution is a Gaussian distribution of the current microseismic theoretical magnitude, with a distribution calculation formula of
- z 0 is a theoretical magnitude
- ⁇ is a Gaussian distribution value
- the magnitude ranges from -3 to 3
- the resolution is 0.1
- the network has a total of 61 labels
- the displacement M(z) corresponds to the current magnitude
- Data augmentation data-augmenting the initial microseismic data: magnifying or shrinking the synthetic DAS microseismic strain data according to different amplitude scales; constructing abnormal traces (bad traces, missing traces); and magnifying or shrinking the background noise data according to different amplitude scales.
- Data preprocessings performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- the magnitude calculation module consists of two parts, including two input branches of frequency spectrum and time waveform, the two input branches use a three-layer convolution structure (convolution layer, pooling layer, activation function, fully connected layer), to extract frequency characteristic and waveform characteristic of the microseismic event, and then a model fusion is performed, and then 2 fully connected layers are used, and finally a calculated magnitude is outputted.
- a three-layer convolution structure convolution layer, pooling layer, activation function, fully connected layer
- the magnitude calculation module uses a 2D convolutional layer, the convolutional layer group includes a convolutional layer, a max pooling layer, and a Dropout layer, the convolution operation uses a ReLu activation function to extract characteristic from the inputted data, the max pooling layer is right after the convolutional layer, its main functions are dimensionality reduction, feature compression, and network complexity simplification.
- the Dropout layer refers to discarding neurons from the network according to a certain probability, thereby reducing the risk of over-fitting of the model.
- a rectified linear unit is used for activation in each hidden convolutional layer, and the kernel size of the convolution is set to (64 ⁇ 3 ⁇ 3); padding and stride are 1 and 2, respectively.
- the final output layer is a vector of 1*61, corresponding to magnitude -3 to magnitude 3 of seism magnitude, then a maximum value is taken as the current calculated magnitude.
- Module IO The input of the computational neural network model is the single DAS microseismic strain data, and the output is the calculated magnitude. An error between the calculated magnitude and theoretical magnitude Gaussian distribution is used to update network parameters. The network uses a mean squared error (MSE) as a loss function to calculate the error of the neural network.
- MSE mean squared error
- FIG. 3 shows a network structure of the real-time microseismic magnitude calculation method based on deep learning according to the embodiment of the present disclosure.
- the network consists of two sub-networks, which are respectively inputted with the waveform and the corresponding frequency spectrum of 2D data, and the output is a 1D Gaussian function of the magnitude, and the peak value corresponds to the accurate value of the magnitude.
- a neural network model adopt a stochastic gradient descent optimization method; setting a dynamic learning rate, the initial value being set to 0.0001, reduced by half every 50 times; the batch size being set to 40; the number of iterations being 200.
- Preprocessing the field data performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- a real-time microseismic magnitude calculation device based on deep learning wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to any one of the foregoing embodiments.
- the present disclosure proposes a real-time microseismic magnitude calculation method and a corresponding device based on deep learning, which integrates waveform and frequency spectrum characteristics with deep model network, and trains the model through continuous DAS strain records, without the need for converting the DAS strain data, improving the ability to quickly estimate the microseismic magnitude and improving the accuracy of the microseismic magnitude estimation.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Geophysics (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Probability & Statistics with Applications (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Embodiments of the present disclosure provide a real-time microseismic magnitude calculation method based on deep learning and a corresponding device. The method includes: constructing a DAS-based horizontal well microseismic monitoring system; constructing a training data set; constructing a magnitude calculation module, wherein the magnitude calculation module comprises two input branches of frequency spectrum and time waveform, the two input branches use a 3-layer convolution structure to extract frequency characteristic and waveform characteristic of a microseismic event, and then a model fusion is performed, and then 2 fully connected layers are used, and finally a calculated magnitude is outputted; training the magnitude calculation module; and analyzing and processing field data. The microseismic magnitude calculation method in the present disclosure improves the ability to quickly estimate the microseismic magnitude, without the need for converting the strain data, and improves the accuracy of the microseismic magnitude estimation.
Description
- This application claims priority to Chinese Patent Application No. 202211206267.7, filed on Sep. 30, 2022 before the China National Intellectual Property Administration, the disclosure of which is incorporated herein by reference in entirety.
- The present disclosure relates to the field of microseismic monitoring and deep learning. In particular, the disclosure relates to a real-time microseismic magnitude calculation method based on deep learning and a corresponding device, wherein the method includes using a training data set to train a magnitude calculation module, analyzing and processing field data, and inputting the processed data to the magnitude calculation module so as to obtain the calculated magnitude, etc. Therefore, the present disclosure belongs to the field of seismic data processing in seismic monitoring.
- Hydraulic fracturing technology relates to injection of high-pressure fluid into shale reservoirs to generate complex artificial fractures, which can increase the connectivity of reservoirs and increase well production. Monitoring and evaluating the effects of hydraulic fracturing on reservoir substructures in different stages is a prerequisite for efficient exploitation and safe production. Microseismic monitoring technology is an important means for real-time monitoring of hydraulic fracturing stimulation of unconventional resources. Microseismic monitoring technology evaluates fracturing effects by monitoring microseismic signals generated during hydraulic fracturing, and then improves and optimizes engineering parameters.
- Quickly and accurately estimating the magnitude of a microseismic event is an important part of microseismic monitoring, and the magnitude of a microseismic event is usually around -3 to 3 Ml. In the earthquake study, key indicators of the traditional methods are extracted from initial earthquake motion information as characteristic parameters, and then an empirical formula is used to relate the characteristic parameters to the magnitude of the earthquake. There are many options for selecting characteristic parameters, which may include amplitude parameters, period (frequency) parameters, and energy parameters of the initial earthquake motion. Although it has been proved that there is a generally linear relationship between the above characteristic parameters and the magnitude, some characteristic parameters are only loosely related to the magnitude, leading to errors in magnitude estimation results. Magnitude estimation results often lead to overestimation of small earthquakes and underestimation of larger earthquakes, to a greater or lesser extent. Further, we can try to estimate the magnitude by using two or more of the above three types of parameters in combination, so as to improve the accuracy of earthquake magnitude estimation. Magnitude estimation methods also include methods of characteristic frequency parameters based on period, such as the maximum predominant period Pd (peak displacement) method, the average predominant period τc method, and so on. In addition, if there are data from multiple stations, the magnitudes from each station can be estimated, and then the weighted average can be used to obtain a magnitude estimate.
- Since earthquake magnitude determination is a very complicated process, it is necessary to comprehensively consider the influence of multiple factors such as source process, propagation medium, site condition, and instrument response. Although a single amplitude or periodic parameter reflects the scale of an earthquake, it is difficult to capture all the characteristics of the earthquake, resulting in insufficient accuracy of earthquake magnitude determination. In recent years, with the rapid development of deep learning technology, more and more researchers have begun to apply machine learning and deep learning to various elements of earthquake waveform data research, and artificial intelligence has been introduced into the research of earthquake magnitude estimation. Deep learning-based techniques can generally be divided into two categories. One is to combine multiple earthquake pre-warning parameters with machine learning methods to estimate the final magnitude. Researchers have integrated various types of 12 initial characteristic parameters of P wave (amplitude parameters, period parameters, energy parameters, and derivative parameters) as an input of a convolutional neural network, and constructed a DCNN-M model in 3 s time period after the arrival of P wave. This has been used for earthquake pre-warning magnitude estimation, with better results than the traditional single parameter method. Ochoa (2018) used 25 parameters (parameters related to azimuth, parameters related to magnitude, parameters related to epicentral distance) obtained after 5 seconds of P wave triggering from a single station to estimate the Richter magnitude scale, based on the Support Vector Regression algorithm. For earthquakes below 4.5, this magnitude estimation has a small error. Another method is to directly combine earthquake waveform data with machine learning methods to improve magnitude estimation. For example, a recently proposed method uses strong motion data to establish a NN network for earthquake pre-warning magnitude estimation, using waveform spectrum information as model input. It was found that, within the 3 s time period, this method is superior to the traditional single parameter method, and it is also believed that the use of multi-characteristic combination input will lead to a better magnitude estimation model. Mousavi & Beroza (2020) constructed a magnitude estimation model composed of a convolutional neural network and a recurrent neural network based on a 30 seconds of single-station waveform input, which effectively improved the magnitude estimation results. Furthermore, a convolutional neural network model with multiple fully connected layers is established, and the frequency domain data of P wave at different periods of 3-9 s, the source information (epicentral distance and source depth) and site information of the corresponding earthquake event are selected as full data set, to train the proposed model and evaluate the estimation performance. Lomax (2019) used three-component acceleration time history data of a single station as input data of a first convolutional layer, normalized the acceleration peak value in the three-component band, and used the normalized value as auxiliary input parameters of the fully connected layer, then using the above two parts of data to iteratively train the model to achieve the purpose of automatic detection and rapid characterization of earthquake magnitude.
- However, at present, there is no magnitude calculation method which uses microseismic strain data monitored from offset-wells by DAS (Distributed Acoustic Sensor). After searching patent databases, there is only a magnitude calculation method using conventional geophone seismic data (generally displacement, velocity or acceleration data, not strain data) and a magnitude calculation method using conventional geophone microseismic data.
- The purpose of the present disclosure is to at least partly overcome the defects of the prior art, and provide a real-time microseismic magnitude calculation method and a corresponding device based on deep learning.
- In an embodiment, the present disclosure provides a real-time microseismic magnitude calculation method based on deep learning, in which the strain data collected by a DAS system may be inputted directly without converting the DAS strain data to conventional data (generally displacement, velocity or acceleration data), thereby improving the ability to quickly estimate the microseismic magnitude.
- In an embodiment, the present disclosure provides a real-time microseismic magnitude calculation method based on deep learning to improve the accuracy of the microseismic magnitude estimation.
- In an embodiment, the present disclosure provides a real-time microseismic magnitude calculation device based on deep learning, which is used for implementing the aforementioned real-time microseismic magnitude calculation method based on deep learning.
- In order to achieve the above-mentioned embodiments, the technical solutions of the present disclosure are as follows:
- A real-time microseismic magnitude calculation method based on deep learning, comprising:
- constructing a DAS-based horizontal well microseismic monitoring system;
- constructing a training data set, comprising DAS microseismic strain data and theoretical magnitude Gaussian distribution label;
- constructing a magnitude calculation module, wherein the magnitude calculation module comprises two input branches of frequency spectrum and time waveform, the two input branches use a 3-layer convolution structure to extract frequency characteristic and waveform characteristic of a microseismic event, and then a model fusion is performed, and then 2 fully connected layers are used, and finally a calculated magnitude is outputted; wherein an input of the magnitude calculation module is a single DAS microseismic strain data, and an output is the calculated magnitude, and an error between the calculated magnitude and a theoretical magnitude Gaussian distribution is used to update network parameters; wherein the network uses a mean square error as a loss function to calculate an error of a neural network;
- training the magnitude calculation module by using the training data set; and
- analyzing and processing field data, comprising: preprocessing the field data; and inputting the preprocessed data into the magnitude calculation module to obtain the calculated magnitude.
- According to some embodiments of the present disclosure, the training data set comprises Nk DAS microseismic strain data, each DAS microseismic strain data comprises Nt channels, each channel comprises Ns sampling points, forming a data body of Nk*Nt*Ns, where Nk, Nt and Ns are all 1000.
- According to some embodiments of the present disclosure, the theoretical magnitude Gaussian distribution label is a Gaussian distribution centered on the theoretical magnitude, the magnitude ranges from -3 to 3, the resolution is 0.1, and it is a vector with 61 values.
- According to some embodiments of the present disclosure, the constructing a training data set comprises:
- generating initial microseismic data;
- data-augmenting the initial microseismic data; and
- preprocessing all data.
- According to some embodiments of the present disclosure, the single initial microseismic data is generated by:
- setting acquisition parameters according to the real microseismic monitoring system, comprising setting source parameters which comprises source mechanism parameters, source mechanism range and source location parameters;
- using a analytical Green’s function under a uniform velocity model to obtain synthetic DAS microseismic strain data, with a size of Nt*Ns;
- collecting monitored background noise data, adding the background noise data to the above-mentioned synthetic DAS microseismic strain data, simulating a signal-to-noise ratio and characteristics of collected filed data, to obtain the synthetic microseismic data.
- According to some embodiments of the present disclosure, the magnitude distribution is a Gaussian distribution of the current microseismic theoretical magnitude, with a distribution calculation formula of
-
- where z0 is a theoretical magnitude, σ is a Gaussian distribution value; the magnitude ranges from -3 to 3, the resolution is 0.1, the number of type labels is 61, and the displacement M(z) corresponds to the current magnitude.
- According to some embodiments of the present disclosure, the data-augmenting the initial microseismic data comprises:
- magnifying or shrinking the synthetic DAS microseismic strain data according to different amplitude scales;
- constructing abnormal traces; and
- magnifying or shrinking the background noise data according to different amplitude scales.
- According to some embodiments of the present disclosure, the preprocessing all data comprises:
- performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- According to some embodiments of the present disclosure, the magnitude calculation module uses a 2D convolutional layer, the convolutional layer group comprises a convolutional layer, a max pooling layer, and a Dropout layer, the convolution operation uses a ReLu activation function to extract characteristic from the inputted data, the max pooling layer is after the convolutional layer, a rectified linear unit is used for activation in each hidden convolutional layer, and the kernel size of the convolution is set to 64×3×3; the final output layer is a vector of 1*61, corresponding to magnitude -3 to
magnitude 3 of seism magnitude, then a maximum value is taken as the current calculated magnitude. - According to some embodiments of the present disclosure, the 3-layer convolutional structure comprises a convolutional layer, a pooling layer and a fully connected layer.
- According to some embodiments of the present disclosure, the training the magnitude calculation module by using the training data set comprises:
- dividing the data set into a training set and a test set with a ratio of 8:2;
- letting a neural network model adopt a stochastic gradient descent optimization method; setting a dynamic learning rate, the initial value being set to 0.0001, reduced by half every 50 times; the batch size being set to 40; the number of iterations being 200; and
- training the neural network model on a graphics processing unit (GPU).
- According to some embodiments of the present disclosure, the DAS-based horizontal well microseismic monitoring system comprises a metal sleeve and an armored optical cable fixedly arranged outside the metal sleeve, and a single-mode optical fiber is arranged in the armored optical cable;
- wherein constructing a DAS-based horizontal well microseismic monitoring system comprises:
- inserting the metal sleeve and the armored optical cable into a drilled wellbore synchronously;
- pumping cement slurry from a bottom of a well, letting the cement slurry return to a wellhead from the bottom of the well along an annular space between an outer wall of the metal sleeve and the wellbore, allowing the metal sleeve, the armored optical cable and formation rock to be fixed together after the cement slurry is solidified; and
- providing a DAS modem, and connecting the single-mode optical fiber in the armored optical cable to a DAS signal input port of the DAS modem at the wellhead.
- According to some embodiments of the present disclosure, the preprocessing the field data comprises:
- performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- According to another aspect of the present disclosure, there is provided a real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to any one of the foregoing embodiments.
- The present disclosure proposes a real-time microseismic magnitude calculation method and a corresponding device based on deep learning, which integrates waveform and frequency spectrum characteristics with a deep model network, and trains the model through continuous DAS strain records, without the need for converting the DAS strain data to conventional seismic data (generally displacement, velocity or acceleration data), improving the ability to quickly estimate the microseismic magnitude and improving the accuracy of the microseismic magnitude estimation.
-
FIG. 1 is a flowchart of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure; -
FIG. 2 shows a software data processing flow of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure; and -
FIG. 3 is a schematic diagram of a network structure of a real-time microseismic magnitude calculation method based on deep learning according to an embodiment of the present disclosure. - In order to make the above objectives, features, and advantages of the present disclosure more understandable, the present disclosure will be described in further detail below with reference to the accompanying drawings and the specific embodiments.
- Distributed Acoustic Sensing (DAS) is an emerging data acquisition technology that has developed rapidly in recent years. DAS generally consists of a demodulator (including an optical system and a signal acquisition system) and a normal optical fiber or a special optical fiber for sensing. DAS realizes the measurement of strain signal along the axial direction of the optical fiber by detecting a phase change of backward Rayleigh scattering light generated by the laser pulse in the internal scatterer of the optical fiber.
- In DAS hydraulic fracturing monitoring, DAS optical fibers are generally arranged along the horizontal well. Compared with conventional geophone detection, DAS can obtain the measurement data of the whole well section, which greatly increases the azimuth of data collection; DAS has 1 m spatial sampling, which significantly improves the density and spatial continuity of the monitoring data.
- The deep learning method is a good tool to solve difficult problems in the field of seismology and geophysical research. The application of deep learning in the field of seismology continues to expand, for example, to seismic phase identification, first arrival picking, source location, etc.
- The present disclosure proposes a method for quickly estimating seismic magnitude based on deep learning and feature fusion, which integrates waveform and frequency spectrum characteristics with deep model network, and trains the model through continuous DAS strain records, improving the ability to quickly estimate the seismic magnitude. According to a specific embodiment of the present disclosure, a real-time microseismic magnitude calculation method based on deep learning is proposed, and the calculation method is carried out in accordance with
FIG. 1 (FIG. 1 ): - Firstly, a DAS-based horizontal well microseismic monitoring system is constructed. The DAS-based horizontal well microseismic monitoring system comprises a metal sleeve and an armored optical cable fixedly arranged outside the metal sleeve, and a single-mode optical fiber is arranged in the armored optical cable; the DAS modem is placed near a wellhead, the DAS signal port of the DAS modem is connected to the special single-mode optical fiber outside the metal sleeve.
- Specifically, the construction of a DAS-based horizontal well microseismic monitoring system includes:
- slowly and synchronously inserting the metal sleeve and the armored optical cable into a drilled wellbore;
- installing a ring-shaped metal clamp at the joint of two metal sleeves at the wellhead to fix and protect the armored optical cable from moving and/or being damaged during the process of inserting the metal sleeve;
- pumping cement slurry from the bottom of the well by a high-pressure pump truck, letting the cement slurry return to the wellhead from the bottom of the well along an annular space between an outer wall of the metal sleeve and the wellbore, allowing the metal sleeve, the armored optical cable and formation rock to be permanently fixed together after the cement slurry is solidified;
- providing the DAS modem, and connecting the single-mode optical fiber in the armored optical cable to a DAS signal input port of the DAS modem at the wellhead.
- During hydraulic fracturing operation, the DAS-based horizontal well microseismic monitoring system may use the armored optical cables permanently arranged outside the metal sleeve to perform hydraulic fracturing microseismic monitoring. The single-mode optical fibers outside the metal sleeve arranged in the well are used to collect data, the data is transmitted to the DAS modem near the wellhead for demodulation, then the microseismic events are obtained, wherein the microseismic events are generated from the stimulated fracture of the underground reservoir of the offset well or in-well when the continuously- hydraulic fracturing operation is performed. Combining with the longitudinal wave velocity field of the underground reservoir obtained in the preceding step, relevant information of the microseismic events can be obtained, and such information includes the time, three-dimensional spatial position and energy magnitude of the microseismic events.
- Furthermore, focal mechanism analysis and magnitude analysis can be performed according to the recorded longitudinal wave signal characteristics of the microseismic events, so as to obtain fracture mechanism of the most microseismic events, and the distribution characteristic and regularity of fractures after the hydraulic fracturing operation may be analyzed. Based on the above information, the effective and reliable, qualitative and quantitative evaluation of the hydraulic fracturing effect of reservoir can be carried out.
- The real-time microseismic magnitude calculation method based on deep learning further includes:
- Constructing a training data set:
- Overview of the data set: The structure of the data set is like this: (2D microseismic data, 1D theoretical magnitude Gaussian distribution label). For the single DAS microseismic strain data, a total of Nt (1000) channels are included, each channel includes Ns (1000) sampling points, and the training data set includes Nk (1000) DAS microseismic strain data. The theoretical magnitude Gaussian distribution label is a 1D Gaussian distribution centered on the theoretical magnitude. The magnitude ranges from -3 to 3 Ml, and the resolution is 0.1, and it is a vector with 61 values. Therefore, the data body is Nk*Nt*Ns and the number of the type labels is Nk*61.
- Data generation: generating initial microseismic data: the single initial microseismic data is generated by:
- setting acquisition parameters according to the real microseismic monitoring system, comprising: setting source parameters, which includes source mechanism parameters and range (only consider the case of double-couple source): dip, rake, strike, wherein these parameters cover the range of the main source mechanism; source location parameters: sx, sy, sz, wherein the location should not be located near the fracturing section, to simulate possible seismic source locations;
- using a analytical Green’s function under a uniform velocity model (Nx*Ny*Nz) to obtain synthetic DAS microseismic strain data, with a size of Nt*Ns;
- collecting monitored background noise data, adding the background noise data to the above-mentioned synthetic DAS microseismic strain data, simulating a signal-to-noise ratio and characteristics of collected field data, to obtain the initial microseismic data.
- The magnitude distribution is a Gaussian distribution of the current microseismic theoretical magnitude, with a distribution calculation formula of
-
- where z0 is a theoretical magnitude, σ is a Gaussian distribution value; the magnitude ranges from -3 to 3, the resolution is 0.1, so the network has a total of 61 labels, and the displacement M(z) corresponds to the current magnitude.
- Data augmentation: data-augmenting the initial microseismic data: magnifying or shrinking the synthetic DAS microseismic strain data according to different amplitude scales; constructing abnormal traces (bad traces, missing traces); and magnifying or shrinking the background noise data according to different amplitude scales.
- Data preprocessings: performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- Constructing a magnitude calculation module:
- 2.1. Module architecture: As shown in
FIG. 3 , the magnitude calculation module consists of two parts, including two input branches of frequency spectrum and time waveform, the two input branches use a three-layer convolution structure (convolution layer, pooling layer, activation function, fully connected layer), to extract frequency characteristic and waveform characteristic of the microseismic event, and then a model fusion is performed, and then 2 fully connected layers are used, and finally a calculated magnitude is outputted. The magnitude calculation module uses a 2D convolutional layer, the convolutional layer group includes a convolutional layer, a max pooling layer, and a Dropout layer, the convolution operation uses a ReLu activation function to extract characteristic from the inputted data, the max pooling layer is right after the convolutional layer, its main functions are dimensionality reduction, feature compression, and network complexity simplification. The Dropout layer refers to discarding neurons from the network according to a certain probability, thereby reducing the risk of over-fitting of the model. A rectified linear unit is used for activation in each hidden convolutional layer, and the kernel size of the convolution is set to (64×3×3); padding and stride are 1 and 2, respectively. The final output layer is a vector of 1*61, corresponding to magnitude -3 tomagnitude 3 of seism magnitude, then a maximum value is taken as the current calculated magnitude. - 2.2. Module IO: The input of the computational neural network model is the single DAS microseismic strain data, and the output is the calculated magnitude. An error between the calculated magnitude and theoretical magnitude Gaussian distribution is used to update network parameters. The network uses a mean squared error (MSE) as a loss function to calculate the error of the neural network.
-
FIG. 3 shows a network structure of the real-time microseismic magnitude calculation method based on deep learning according to the embodiment of the present disclosure. The network consists of two sub-networks, which are respectively inputted with the waveform and the corresponding frequency spectrum of 2D data, and the output is a 1D Gaussian function of the magnitude, and the peak value corresponds to the accurate value of the magnitude. - Training the magnitude calculation module by using the training data set:
- 3.1. Dividing the data set into training set and test set with a ratio of 8:2.
- 3.2. Letting a neural network model adopt a stochastic gradient descent optimization method; setting a dynamic learning rate, the initial value being set to 0.0001, reduced by half every 50 times; the batch size being set to 40; the number of iterations being 200.
- 3.3. Training the neural network model on a graphics processing unit (GPU).
- Analyzing and processing field data:
- 4.1. Preprocessing the field data: performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
- 4.2. Inputting the preprocessed data into the magnitude calculation module to obtain the calculated magnitude.
- According to another aspect of the present disclosure, there is provided a real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to any one of the foregoing embodiments.
- The present disclosure proposes a real-time microseismic magnitude calculation method and a corresponding device based on deep learning, which integrates waveform and frequency spectrum characteristics with deep model network, and trains the model through continuous DAS strain records, without the need for converting the DAS strain data, improving the ability to quickly estimate the microseismic magnitude and improving the accuracy of the microseismic magnitude estimation.
- While the embodiments of the present disclosure have been shown and described above, it is to be understood that the above embodiments are illustrative and are not to be construed as limiting the scope of the present disclosure. Changes, modifications, substitutions and variations may be made to the above embodiments by those skilled in the art within the scope of the present disclosure.
Claims (14)
1. A real-time microseismic magnitude calculation method based on deep learning, comprising:
constructing a DAS-based horizontal well microseismic monitoring system;
constructing a training data set, comprising DAS microseismic strain data and theoretical magnitude Gaussian distribution label;
constructing a magnitude calculation module, wherein the magnitude calculation module comprises two input branches of frequency spectrum and time waveform, the two input branches use a 3-layer convolution structure to extract frequency characteristic and waveform characteristic of a microseismic event, and then a model fusion is performed, and then 2 fully connected layers are used, and finally a calculated magnitude is outputted; wherein an input of the magnitude calculation module is the DAS microseismic strain data, and its output is the calculated magnitude, and error between the calculated magnitude and a theoretical magnitude Gaussian distribution is used to update network parameters; wherein the network uses a mean square error as a loss function to calculate an error of a neural network;
training the magnitude calculation module by using the training data set; and
analyzing and processing field data, comprising: preprocessing the field data; and inputting the preprocessed data into the magnitude calculation module to obtain the calculated magnitude,
wherein the constructing a training data set comprises: generating initial microseismic data; data-augmenting the initial microseismic data; and preprocessing all data;
wherein the magnitude calculation module uses a 2D convolutional layer, the convolutional layer group comprises a convolutional layer, a max pooling layer, and a Dropout layer, the convolution operation uses a ReLu activation function to extract characteristic from the inputted data, the max pooling layer is after the convolutional layer, a rectified linear unit is used for activation in each hidden convolutional layer, and the kernel size of the convolution is set to 64x3 x3; the final output layer is a vector of 1*61, corresponding to magnitude -3 to magnitude 3 of seism magnitude, then a maximum value is taken as the current calculated magnitude.
2. The real-time microseismic magnitude calculation method based on deep learning according to claim 1 , wherein the training data set comprises Nk DAS microseismic strain data, each DAS microseismic strain data comprises Nt channels, each channel comprises Ns sampling points, forming a data body of Nk*Nt*Ns, where Nk, Nt and Ns are all 1000.
3. The real-time microseismic magnitude calculation method based on deep learning according to claim 2 , wherein the single initial microseismic data is generated by:
setting acquisition parameters according to the real microseismic monitoring system, comprising setting source parameters which comprises source mechanism parameters, source mechanism range and source location parameters;
using a analytical Green’s function under a uniform velocity model to obtain synthetic DAS microseismic strain data, with a size of Nt*Ns;
collecting monitored field background noise data, adding the background noise data to the above-mentioned synthetic DAS microseismic strain data, simulating a signal-to-noise ratio and characteristics of collected field data, to obtain the initial microseismic data.
4. The real-time microseismic magnitude calculation method based on deep learning according to claim 3 , wherein the data-augmenting the initial microseismic data comprises:
magnifying or shrinking the synthetic DAS microseismic strain data according to different amplitude scales;
constructing abnormal traces; and
magnifying or shrinking the background noise data according to different amplitude scales.
5. The real-time microseismic magnitude calculation method based on deep learning according to claim 2 , wherein the preprocessing all data comprises:
performing interpolation and replacement to damaged channel data by adjacent channels, removing mean value, removing abnormally large value, and replacing the removed abnormally large value by interpolation.
6. The real-time microseismic magnitude calculation method based on deep learning according to claim 1 , wherein the training the magnitude calculation module by using the training data set comprises:
dividing the data set into a training set and a test set with a ratio of 8:2;
letting a neural network model adopt a stochastic gradient descent optimization method; setting a dynamic learning rate, the initial value being set to 0.0001, reduced by half every 50 times; the batch size being set to 40; the number of iterations being 200; and
training the neural network model on a graphics processing unit (GPU).
7. The real-time microseismic magnitude calculation method based on deep learning according to claim 1 , where:
the DAS-based horizontal well microseismic monitoring system comprises a metal sleeve and an armored optical cable fixedly arranged outside the metal sleeve, and a single-mode optical fiber is arranged in the armored optical cable;
wherein constructing a DAS-based horizontal well microseismic monitoring system comprises:
inserting the metal sleeve and the armored optical cable into a drilled wellbore synchronously;
pumping cement slurry from a bottom of a well, letting the cement slurry return to a wellhead from the bottom of the well along an annular space between an outer wall of the metal sleeve and the wellbore, allowing the metal sleeve, the armored optical cable and formation rock to be fixed together after the cement slurry is solidified; and
providing a DAS modem, and connecting the single-mode optical fiber in the armored optical cable to a DAS signal input port of the DAS modem at the wellhead.
8. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 1 .
9. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 2 .
10. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 3 .
11. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 4 .
12. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 5 .
13. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 6 .
14. A real-time microseismic magnitude calculation device based on deep learning, wherein, the calculation device is used for implementing the real-time microseismic magnitude calculation method based on deep learning according to claim 7 .
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211206267.7A CN115291281B (en) | 2022-09-30 | 2022-09-30 | Real-time micro-earthquake magnitude calculation method and device based on deep learning |
CN202211206267.7 | 2022-09-30 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20230324577A1 true US20230324577A1 (en) | 2023-10-12 |
US11789173B1 US11789173B1 (en) | 2023-10-17 |
Family
ID=83834949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US18/304,026 Active US11789173B1 (en) | 2022-09-30 | 2023-04-20 | Real-time microseismic magnitude calculation method and device based on deep learning |
Country Status (2)
Country | Link |
---|---|
US (1) | US11789173B1 (en) |
CN (1) | CN115291281B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20240027639A1 (en) * | 2022-07-19 | 2024-01-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Seismometer with high sensitivity, broadband and all-dip |
CN117949687A (en) * | 2024-03-27 | 2024-04-30 | 山东省科学院激光研究所 | Wind speed measurement method and system based on distributed optical fiber sensing and deep learning |
CN118033744A (en) * | 2024-04-11 | 2024-05-14 | 山东省地震局 | Method and system for calculating seismic source parameters in real time |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116340757B (en) * | 2023-04-25 | 2023-08-08 | 中国地震局地震研究所 | Characteristic self-adaptive earthquake early-warning magnitude prediction method and system thereof |
CN117313557B (en) * | 2023-11-29 | 2024-01-23 | 中国石油大学(华东) | Deep learning-based fracturing magnetic method real-time monitoring middle fracture diagnosis method |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA3120493A1 (en) * | 2018-11-29 | 2020-06-04 | Bp Exploration Operating Company Limited | Das data processing to identify fluid inflow locations and fluid type |
CN110058305A (en) * | 2019-05-24 | 2019-07-26 | 吉林大学 | A kind of DAS seismic data noise-reduction method based on convolutional neural networks |
US20210089897A1 (en) * | 2019-09-24 | 2021-03-25 | Quantico Energy Solutions Llc | High-resolution earth modeling using artificial intelligence |
CN111458748B (en) * | 2020-03-30 | 2021-03-05 | 青岛理工大学 | Performance earthquake motion risk analysis method based on three-layer data set neural network |
CN111538076B (en) * | 2020-05-13 | 2021-08-03 | 浙江大学 | Earthquake magnitude rapid estimation method based on deep learning feature fusion |
CN111983676A (en) * | 2020-07-07 | 2020-11-24 | 东华理工大学 | Earthquake monitoring method and device based on deep learning |
CN111665568A (en) * | 2020-07-09 | 2020-09-15 | 中油奥博(成都)科技有限公司 | Micro-logging device and measuring method based on distributed optical fiber acoustic wave sensing technology |
US11699099B2 (en) * | 2020-10-28 | 2023-07-11 | Quantico Energy Solutions Llc | Confidence volumes for earth modeling using machine learning |
CN112946727B (en) * | 2021-01-28 | 2021-09-14 | 成都高新减灾研究所 | Earthquake monitoring method, earthquake early warning method and system |
CN112782762A (en) * | 2021-01-29 | 2021-05-11 | 东北大学 | Earthquake magnitude determination method based on deep learning |
CN113514877B (en) * | 2021-07-07 | 2022-05-31 | 浙江大学 | Self-adaptive quick earthquake magnitude estimation method |
US20230054654A1 (en) * | 2021-08-20 | 2023-02-23 | Silixa Ltd. | Fracture detection using distributed optical fiber sensing |
CN114167487B (en) * | 2021-12-02 | 2022-09-27 | 中国地震局工程力学研究所 | Seismic magnitude estimation method and device based on characteristic waveform |
CN115032693B (en) * | 2022-07-11 | 2023-02-03 | 中国地震局地球物理研究所 | Method and device for automatically identifying pre-earthquake of strong earthquake |
CN114879252B (en) * | 2022-07-11 | 2022-09-13 | 中国科学院地质与地球物理研究所 | DAS (data acquisition system) same-well monitoring real-time microseism effective event identification method based on deep learning |
-
2022
- 2022-09-30 CN CN202211206267.7A patent/CN115291281B/en active Active
-
2023
- 2023-04-20 US US18/304,026 patent/US11789173B1/en active Active
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20240027639A1 (en) * | 2022-07-19 | 2024-01-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Seismometer with high sensitivity, broadband and all-dip |
US11953635B2 (en) * | 2022-07-19 | 2024-04-09 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Seismometer with high sensitivity, broadband and all-dip |
CN117949687A (en) * | 2024-03-27 | 2024-04-30 | 山东省科学院激光研究所 | Wind speed measurement method and system based on distributed optical fiber sensing and deep learning |
CN118033744A (en) * | 2024-04-11 | 2024-05-14 | 山东省地震局 | Method and system for calculating seismic source parameters in real time |
Also Published As
Publication number | Publication date |
---|---|
US11789173B1 (en) | 2023-10-17 |
CN115291281B (en) | 2022-12-20 |
CN115291281A (en) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11789173B1 (en) | Real-time microseismic magnitude calculation method and device based on deep learning | |
CN110109895B (en) | Surrounding rock grading combined prediction method suitable for TBM tunneling tunnel and application | |
CN106194159B (en) | Mine inclination measurement while drilling exploration system and measuring method thereof | |
CN102373923B (en) | Reservoir stratum identification method | |
US11899154B2 (en) | DAS same-well monitoring real-time microseismic effective event identification method based on deep learning | |
CN102736124B (en) | Tunnel excavation surrounding rock dynamic refined classification method based on integrated parameters | |
CN106855636A (en) | Based on the prototype geological model Seismic forward method that carbonate reservoir is appeared | |
CN116591777B (en) | Multi-field multi-source information fusion rock burst intelligent monitoring and early warning device and method | |
CN116378648B (en) | Near-bit stratum detection method and device based on while-drilling acoustic wave forward looking | |
CN105719433A (en) | In-hole seismic wave based advanced prediction method | |
CN109854303A (en) | Rock burst monitoring and early warning method for underground mining | |
CN106646661A (en) | Comprehensive survey system for hydrogeology of mineral deposit | |
US20080112263A1 (en) | System and method for determining seismic event location | |
CN113189644A (en) | Microseismic source positioning method and system | |
US20240230942A9 (en) | Method and System for Real-Time Calculating a Microseismic Focal Mechanism Based on Deep Learning | |
CN117590466A (en) | Micro fault detection method and fault detection device | |
CN116721223A (en) | Rock mass stability analysis method based on oblique photography and measurement while drilling technology | |
CN116794733A (en) | Tunnel magnetic resonance and transient electromagnetic joint inversion method | |
CN105929445B (en) | Micro rupture vector scan method | |
CN105824042A (en) | Maximum longitudinal distance design method based on optimal lighting energy | |
Wang et al. | Application of Deep Learning in Microseismic Waveform Classification: A Case Study of The Yebatan Hydropower Station Project | |
Paap et al. | A neural network approach for improved seismic event detection in the Groningen gas field, the Netherlands | |
CN109164482A (en) | Underground microseism method for solving, apparatus and system based on fibre optical sensor | |
Chakravarty et al. | Hydraulic Fracturing-Driven Infrasound Signals–a New Class of Signal for Subsurface Engineering | |
CN107797138A (en) | The improvement of forecast precision influence factor and optimization method in a kind of advance geologic prediction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |