An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb...An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.展开更多
Seismic illumination plays an important role in subsurface imaging. A better image can be expected either through optimizing acquisition geometry or introducing more advanced seismic mi- gration and/or tomographic inv...Seismic illumination plays an important role in subsurface imaging. A better image can be expected either through optimizing acquisition geometry or introducing more advanced seismic mi- gration and/or tomographic inversion methods involving illumination compensation. Vertical cable survey is a potential replacement of traditional marine seismic survey for its flexibility and data quality. Conventional vertical cable data processing requires separation of primaries and multiples before migration. We proposed to use multi-scale full waveform inversion （FWI） to improve illumination coverage of vertical cable survey. A deep water velocity model is built to test the capability of multi-scale FWI in detecting low velocity anomalies below seabed. Synthetic results show that multi-scale FWI is an effective model building tool in deep-water exploration. Geometry optimization through target ori- ented illumination analysis and multi-scale FWI may help to mitigate the risks of vertical cable survey. The combination of multi-scale FWI, low-frequency data and multi-vertical-cable acquisition system may provide both high resolution and high fidelity subsurface models.展开更多
An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 199...An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.展开更多
Integrated data and energy transfer(IDET)is capable of simultaneously delivering on-demand data and energy to low-power Internet of Everything(Io E)devices.We propose a multi-carrier IDET transceiver relying on superp...Integrated data and energy transfer(IDET)is capable of simultaneously delivering on-demand data and energy to low-power Internet of Everything(Io E)devices.We propose a multi-carrier IDET transceiver relying on superposition waveforms consisting of multi-sinusoidal signals for wireless energy transfer(WET)and orthogonal-frequency-divisionmultiplexing(OFDM)signals for wireless data transfer(WDT).The outdated channel state information(CSI)in aging channels is employed by the transmitter to shape IDET waveforms.With the constraints of transmission power and WDT requirement,the amplitudes and phases of the IDET waveform at the transmitter and the power splitter at the receiver are jointly optimised for maximising the average directcurrent(DC)among a limited number of transmission frames with the existence of carrier-frequencyoffset(CFO).For the amplitude optimisation,the original non-convex problem can be transformed into a reversed geometric programming problem,then it can be effectively solved with existing tools.As for the phase optimisation,the artificial bee colony(ABC)algorithm is invoked in order to deal with the nonconvexity.Iteration between the amplitude optimisation and phase optimisation yields our joint design.Numerical results demonstrate the advantage of our joint design for the IDET waveform shaping with the existence of the CFO and the outdated CSI.展开更多
Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse ...Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.展开更多
The absence of low-frequency information in seismic data is one of the most difficult problems in elastic full waveform inversion. Without low-frequency data, it is difficult to recover the long-wavelength components ...The absence of low-frequency information in seismic data is one of the most difficult problems in elastic full waveform inversion. Without low-frequency data, it is difficult to recover the long-wavelength components of subsurface models and the inversion converges to local minima. To solve this problem, the elastic envelope inversion method is introduced. Based on the elastic envelope operator that is capable of retrieving low- frequency signals hidden in multicomponent data, the proposed method uses the envelope of multicomponent seismic signals to construct a misfit function and then recover the long- wavelength components of the subsurface model. Numerical tests verify that the elastic envelope method reduces the inversion nonlinearity and provides better starting models for the subsequent conventional elastic full waveform inversion and elastic depth migration, even when low frequencies are missing in multicomponent data and the starting model is far from the true model. Numerical tests also suggest that the proposed method is more effective in reconstructing the long-wavelength components of the S-wave velocity model. The inversion of synthetic data based on the Marmousi-2 model shows that the resolution of conventional elastic full waveform inversion improves after using the starting model obtained using the elastic envelope method. Finally, the limitations of the elastic envelope inversion method are discussed.展开更多
In most effective bits evaluation of waveform recorders, the prerequisite is that there is no signal source distortion, or the distortion can be neglected. But when the distortion can be neglected or how it affects th...In most effective bits evaluation of waveform recorders, the prerequisite is that there is no signal source distortion, or the distortion can be neglected. But when the distortion can be neglected or how it affects the evaluation when it can't be neglected it is not determined yet. In this paper, the influence of signal source distortion to the evaluation of the effective bits of waveform recorders is discussed, then, the correction method of the effective bits error caused by the distortion influence is given. Finally , the error limit of the effective bits is given and how to selecte the calibrator is introduced. In the end , some simulation results of the new method in test are described.展开更多
Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to a...Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to address this issue and present a fast automatic seismic waveform classification method to pick densely-sampled FBs directly from common-shot gathers with sparsely distributed traces.Through feeding a large number of representative shot gathers with missing traces and the corresponding binary labels segmented by manually interpreted fully-sampled FBs,we can obtain a welltrained Seg Net model.When any unseen gather including the one with irregular trace spacing is inputted,the Seg Net can output the probability distribution of different categories for waveform classification.Then FBs can be picked by locating the boundaries between one class on post-FBs data and the other on pre-FBs background.Two land datasets with each over 2000 shots are adopted to illustrate that one well-trained 25-layer Seg Net can favorably classify waveform and further pick fully-sampled FBs verified by the manually-derived ones,even when the proportion of randomly missing traces reaches50%,21 traces are missing consecutively,or traces are missing regularly.展开更多
Seismic waveform clustering is a useful technique for lithologic identification and reservoir characterization.The current seismic waveform clustering algorithms are predominantly based on a fixed time window,which is...Seismic waveform clustering is a useful technique for lithologic identification and reservoir characterization.The current seismic waveform clustering algorithms are predominantly based on a fixed time window,which is applicable for layers of stable thickness.When a layer exhibits variable thickness in the seismic response,a fixed time window cannot provide comprehensive geologic information for the target interval.Therefore,we propose a novel approach for a waveform clustering workfl ow based on a variable time window to enable broader applications.The dynamic time warping(DTW)distance is fi rst introduced to effectively measure the similarities between seismic waveforms with various lengths.We develop a DTW distance-based clustering algorithm to extract centroids,and we then determine the class of all seismic traces according to the DTW distances from centroids.To greatly reduce the computational complexity in seismic data application,we propose a superpixel-based seismic data thinning approach.We further propose an integrated workfl ow that can be applied to practical seismic data by incorporating the DTW distance-based clustering and seismic data thinning algorithms.We evaluated the performance by applying the proposed workfl ow to synthetic seismograms and seismic survey data.Compared with the the traditional waveform clustering method,the synthetic seismogram results demonstrate the enhanced capability of the proposed workfl ow to detect boundaries of diff erent lithologies or lithologic associations with variable thickness.Results from a practical application show that the planar map of seismic waveform clustering obtained by the proposed workfl ow correlates well with the geological characteristics of wells in terms of reservoir thickness.展开更多
文摘An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.
基金the financial support by the National Natural Science Foundation of China (Nos.41304109 and 41230318)the Fundamental Research Funds for the Central Universities,China University of Geosciences (Wuhan) (Nos.CUG130103 and CUG110803)
文摘Seismic illumination plays an important role in subsurface imaging. A better image can be expected either through optimizing acquisition geometry or introducing more advanced seismic mi- gration and/or tomographic inversion methods involving illumination compensation. Vertical cable survey is a potential replacement of traditional marine seismic survey for its flexibility and data quality. Conventional vertical cable data processing requires separation of primaries and multiples before migration. We proposed to use multi-scale full waveform inversion （FWI） to improve illumination coverage of vertical cable survey. A deep water velocity model is built to test the capability of multi-scale FWI in detecting low velocity anomalies below seabed. Synthetic results show that multi-scale FWI is an effective model building tool in deep-water exploration. Geometry optimization through target ori- ented illumination analysis and multi-scale FWI may help to mitigate the risks of vertical cable survey. The combination of multi-scale FWI, low-frequency data and multi-vertical-cable acquisition system may provide both high resolution and high fidelity subsurface models.
文摘An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.
基金financial support of Natural Science Foundation of China(No.61971102,62132004)MOST Major Research and Development Project(No.2021YFB2900204)+1 种基金Sichuan Science and Technology Program(No.2022YFH0022)Key Research and Development Program of Zhejiang Province(No.2022C01093)。
文摘Integrated data and energy transfer(IDET)is capable of simultaneously delivering on-demand data and energy to low-power Internet of Everything(Io E)devices.We propose a multi-carrier IDET transceiver relying on superposition waveforms consisting of multi-sinusoidal signals for wireless energy transfer(WET)and orthogonal-frequency-divisionmultiplexing(OFDM)signals for wireless data transfer(WDT).The outdated channel state information(CSI)in aging channels is employed by the transmitter to shape IDET waveforms.With the constraints of transmission power and WDT requirement,the amplitudes and phases of the IDET waveform at the transmitter and the power splitter at the receiver are jointly optimised for maximising the average directcurrent(DC)among a limited number of transmission frames with the existence of carrier-frequencyoffset(CFO).For the amplitude optimisation,the original non-convex problem can be transformed into a reversed geometric programming problem,then it can be effectively solved with existing tools.As for the phase optimisation,the artificial bee colony(ABC)algorithm is invoked in order to deal with the nonconvexity.Iteration between the amplitude optimisation and phase optimisation yields our joint design.Numerical results demonstrate the advantage of our joint design for the IDET waveform shaping with the existence of the CFO and the outdated CSI.
基金Supported by the National Natural Science Foundation of China (69983005)
文摘Frequency-domain waveform seismic tomography includes modeling of wave propagation and full waveform inversion of correcting the initial velocity model. In the forward modeling, we use direct solution based on sparse matrix factorization, combined with nine-point finite-difference for the linear system of equations. In the waveform inversion, we use preconditioned gradient method where the preconditioner is provided by the diagonal of the approximate Hessian matrix. We successfully applied waveform inversion method from low to high frequency in two sets of Marmousi data. One is the data set generated by frequencydomain finite-difference modeling, and the other is the original Marmousi shots data set. The former result is very close to the true velocity model. In the original shots data set inversion, we replace the prior source with estimated source; the result is also acceptable, and consistent with the true model.
文摘The absence of low-frequency information in seismic data is one of the most difficult problems in elastic full waveform inversion. Without low-frequency data, it is difficult to recover the long-wavelength components of subsurface models and the inversion converges to local minima. To solve this problem, the elastic envelope inversion method is introduced. Based on the elastic envelope operator that is capable of retrieving low- frequency signals hidden in multicomponent data, the proposed method uses the envelope of multicomponent seismic signals to construct a misfit function and then recover the long- wavelength components of the subsurface model. Numerical tests verify that the elastic envelope method reduces the inversion nonlinearity and provides better starting models for the subsequent conventional elastic full waveform inversion and elastic depth migration, even when low frequencies are missing in multicomponent data and the starting model is far from the true model. Numerical tests also suggest that the proposed method is more effective in reconstructing the long-wavelength components of the S-wave velocity model. The inversion of synthetic data based on the Marmousi-2 model shows that the resolution of conventional elastic full waveform inversion improves after using the starting model obtained using the elastic envelope method. Finally, the limitations of the elastic envelope inversion method are discussed.
文摘In most effective bits evaluation of waveform recorders, the prerequisite is that there is no signal source distortion, or the distortion can be neglected. But when the distortion can be neglected or how it affects the evaluation when it can't be neglected it is not determined yet. In this paper, the influence of signal source distortion to the evaluation of the effective bits of waveform recorders is discussed, then, the correction method of the effective bits error caused by the distortion influence is given. Finally , the error limit of the effective bits is given and how to selecte the calibrator is introduced. In the end , some simulation results of the new method in test are described.
基金financially supported by the National Key R&D Program of China(2018YFA0702504)the Fundamental Research Funds for the Central Universities(2462019QNXZ03)+1 种基金the National Natural Science Foundation of China(42174152 and 41974140)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX 2020-03)。
文摘Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to address this issue and present a fast automatic seismic waveform classification method to pick densely-sampled FBs directly from common-shot gathers with sparsely distributed traces.Through feeding a large number of representative shot gathers with missing traces and the corresponding binary labels segmented by manually interpreted fully-sampled FBs,we can obtain a welltrained Seg Net model.When any unseen gather including the one with irregular trace spacing is inputted,the Seg Net can output the probability distribution of different categories for waveform classification.Then FBs can be picked by locating the boundaries between one class on post-FBs data and the other on pre-FBs background.Two land datasets with each over 2000 shots are adopted to illustrate that one well-trained 25-layer Seg Net can favorably classify waveform and further pick fully-sampled FBs verified by the manually-derived ones,even when the proportion of randomly missing traces reaches50%,21 traces are missing consecutively,or traces are missing regularly.
基金supported by the National Science and Technology Major Project (No. 2017ZX05001-003)。
文摘Seismic waveform clustering is a useful technique for lithologic identification and reservoir characterization.The current seismic waveform clustering algorithms are predominantly based on a fixed time window,which is applicable for layers of stable thickness.When a layer exhibits variable thickness in the seismic response,a fixed time window cannot provide comprehensive geologic information for the target interval.Therefore,we propose a novel approach for a waveform clustering workfl ow based on a variable time window to enable broader applications.The dynamic time warping(DTW)distance is fi rst introduced to effectively measure the similarities between seismic waveforms with various lengths.We develop a DTW distance-based clustering algorithm to extract centroids,and we then determine the class of all seismic traces according to the DTW distances from centroids.To greatly reduce the computational complexity in seismic data application,we propose a superpixel-based seismic data thinning approach.We further propose an integrated workfl ow that can be applied to practical seismic data by incorporating the DTW distance-based clustering and seismic data thinning algorithms.We evaluated the performance by applying the proposed workfl ow to synthetic seismograms and seismic survey data.Compared with the the traditional waveform clustering method,the synthetic seismogram results demonstrate the enhanced capability of the proposed workfl ow to detect boundaries of diff erent lithologies or lithologic associations with variable thickness.Results from a practical application show that the planar map of seismic waveform clustering obtained by the proposed workfl ow correlates well with the geological characteristics of wells in terms of reservoir thickness.
