Moment tensor and stress inversion solutions of acoustic emissions during compression and tensile fracturing in crystalline rocks

We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.

Two-dimensional regularized inversion of AMT data based on rotation invariant of Central impedance tensor

Considering the uncertainty of the electrical axis for two-dimensional audo-magnetotelluric(AMT) data processing, an AMT inversion method with the Central impedance tensor was presented. First, we present a calculation expression of the Central impedance tensor in AMT, which can be considered as the arithmetic mean of TE-polarization mode and TM-polarization mode in the twodimensional geo-electrical model. Second, a least-squares iterative inversion algorithm is established, based on a smoothnessconstrained model, and an improved L-curve method is adopted to determine the best regularization parameters. We then test the above inversion method with synthetic data and field data. The test results show that this two-dimensional AMT inversion scheme for the responses of Central impedance is effective and can reconstruct reasonable two-dimensional subsurface resistivity structures. We conclude that the Central impedance tensor is a useful tool for two-dimensional inversion of AMT data.

Application of tensor CSAMT with high-power orthogonal signal sources in Jiama porphyry copper deposit,South Tibet

The Jiama porphyry copper deposit in Tibet is one of the proven supergiant copper deposits in the Qinghai-Tibet Plateau at present,with the reserves of geological resources equivalent to nearly 20×10^(6) t.However,it features wavy and steep terrain,leading to extremely difficult field operation and heavy interference.This study attempts to determine the effects of the tensor controlled-source audiomagnetotellurics(CSAMT)with high-power orthogonal signal sources(also referred to as the high-power tensor CSAMT)when it is applied to the deep geophysical exploration in plateaus with complex terrain and mining areas with strong interference.The test results show that the high current provided by the highpower tensor CSAMT not only greatly improved the signal-to-noise ratio but also guaranteed that effective signals were received in the case of a long transmitter-receiver distance.Meanwhile,the tensor data better described the anisotropy of deep geologic bodies.In addition,the tests also show that when the transmitting current reaches 60 A,it is still guaranteed that strong enough signals can be received in the case of the transmitter-receiver distance of about 25 km,sounding curves show no near field effect,and effective exploration depth can reach 3 km.The 2D inversion results are roughly consistent with drilling results,indicating that the high-power tensor CSAMT can be used to achieve nearly actual characteristics of underground electrical structures.Therefore,this method has great potential for application in deep geophysical exploration in plateaus and mining areas with complex terrain and strong interference,respectively.This study not only serves as important guidance on the prospecting in the Qinghai-Tibet Plateau but also can be used as positive references for deep mineral exploration in other areas.

Moment tensor inversion for the focal mechanism of the Dongfang (Hainan) earthquake swarm

Moment tensor inversion for the focal mechanism of the 12 earthquakes of the Dongfang (Hainan) earthquakeswarm occurred from June to August 1992 with near-source broadband data recorded by a temporal small-aperturenetwork consisting of DCS-302 digital three-component accelerographs. The results inverted indicate that thepredominant components of sources of all these 12 earthqualles were shear dislocations. The principal pressureaxis and the principal tension axis are in NW-SE direction and in NE-SW direction, respectively, and their dips arealmost horizontal. It could infer that these earthquakes occurred within the same ambient stress field.

Moment tensor inversion for focal mechanism of the Beibuwan earthquakes

Two earthquakes of Ms=6.0 and Ms=6. 1 consecutively occurred on December 31, 1994 and January 10, 1995 in Beibuwan region, China. By using the generalized reflection-transmission coefficient matrix and the discrete slowness integration method in the calculation of Green’s functions, we obtained the focal mechanisms of these earthquakes using long-period waveforms of regional body waves recorded by the China Digital Seismograph Network (CDSN) by means of moment tensor inversion method in frequency domain. The results inverted indicate that the focal mechanisms of these two earthquakes were similar to each other. Their principal compressional stresses are in NW-SE direction and principal tensional stresses are in NE-SW direction. It turns out that the occurrence of the two earthquakes was controlled by the same tectonic environment related to the collision of the Philippine Plate and the Eurasian Plates. On the other hand, the results imply that the stress field in the seismogenic region has a significant change after the Ms=6.0 earthquake. It may be proposed that the occurrence of the Ms=6. 1 earthquake could be related to the stress field adjustment caused by the Ms=6.0 earthquake.

Moment tensor inversion of the November 6, 1988 MS=7.6, Lancang-Gengma, China,earthquake using long-period body-waves data

Moment tensor inversion was carried out to myert the source mechanism and source time function of the Ms=7.6November 6. 1988, Lancang-Gengma. Yunnan Province, Chin4 earthquake. Waveforms of long-period bodywaves recorded by China Digital Seismograph Network (CDSN) were used in the inversion. The inverted resultshows one nodal plane of right-lateral strike-slip faulting and another of left-lateral strike-slip faulting and a simplesource time function of a duration of about 15 s and scalar seismic moment of 6.4x 102oN-N-m From the geologicaldata and tectonic settings and also from field observations and epicentral distribution of aftershocks, the nodalplane striking in the azimuth of 313° is preferred as the fault plane. The pressure axis lies almost horizontally innorth-south direction.

Applications of seismic moment tensor inversion in fast response to earthquakes

Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere sent as "Bulletins of Source Mechanism Parameters of Earthquakes" to the Seismic Regime Guards’ Office,China Seismological Bureau, and the relevant provincial seismological bureaus. These bulletins have played rolein the fast response to large earthquakes.

Frequency Time Analysis (FTAN) and Moment Tensor Inversion Solutions from Short Period Surface Waves in Cameroon (Central Africa)

Short period surface waves generated by a local earthquake recorded by broadband seismometers at distances of about 186 to 778 km from the earthquake’s epicenter located in Cameroon (Central Africa) were processed for group velocity maps and dispersion waveforms using the frequency time analysis (FTAN) method. The resulting group velocity fundamental modes of the extracted Rayleigh and Love waves were used for a joint amplitude spectral and P polarity inversion using moment tensor inversion. The corresponding group velocity dispersion curves, the residual as a function of depth, the amplitude spectra and the moment tensor solutions of the regions from the epicenter to the different stations up to a depth of about 10 km were obtained.

Moment tensor inversion of small to moderate earthquakes in the Capital Region in 2004

The moment tensor solutions of 51 small to moderate earthquakes occurred in the Capital Region in the year of 2004 are obtained by inverting the broadband waveform data. Accordingly, other source parameters, such as scalar seismic moments, moment magnitudes, double-couple (DC) components and compensated-linear-vector-dipole (CLVD) components, are determined as well as fault parameters and stress-axis parameters. The inverted results are evaluated by groups of numerical tests.

Joint inversion of gravity and multiple components of tensor gravity data

Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.

The effect of focal depth error on moment tensor inversion

ＴｈｅｅｆｅｃｔｏｆｆｏｃａｌｄｅｐｔｈｅｒｏｒｏｎｍｏｍｅｎｔｔｅｎｓｏｒｉｎｖｅｒｓｉｏｎＬＩＳＨＥＮＧＸＵ（许力生）ＹＵＮＴＡＩＣＨＥＮ（陈运泰）ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，ＳｔａｔｅＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕ...

Identification and Correction for MT Static Shift Using TEM Inversion Technique

The inversion of TEM data, using the observed magnetic fields instead of that of apparent resistivities data in this paper, avoids the errors caused by the definition of the apparent resistivity. The inversed results by fitting the magnetic fields of the transmitter source’s image with the observed magnetic fields are relatively less affected by the conductivity inhomogeneity. The MT apparent curve is calculated on the basis of the conductivity model constructed from the TEM inversion results. This curve is used as a reference curve for the correction of MT static shift, which makes the correction more reliable. Meanwhile, the domain transformation is also achieved from time to frequency between the two kinds of electromagnetic data. Therefore, the correction of the MT static shift is actualized using TEM inversion method. The corresponding application research shows that this method is very effective for the identification and correction of the MT static shift.

The MT inversion for conductivity anisotropy and EDA precursor，stress field and deformationbandintheEarthsdeepcrust

ＴｈｅＭＴｉｎｖｅｒｓｉｏｎｆｏｒｃｏｎｄｕｃｔｉｖｉｔｙａｎｉｓｏｔｒｏｐｙａｎｄＥＤＡｐｒｅｃｕｒｓｏｒ，ｓｔｒｅｓｆｉｅｌｄａｎｄｄｅｆｏｒ┐ｍａｔｉｏｎｂａｎｄｉｎｔｈｅＥａｒｔｈｓｄｅｅｐｃｒｕｓｔＣＨＡＮＧ－ＹＯＵＬＩＮ（林长佑），Ｃ...

Full moment tensor inversion constrained by doublecouple focal mechanism for induced seismicity

In this study,we propose a new method to determine full moment tensor solution for induced seismicity.This method generalizes the full waveform matching algorithm we have developed to determine the double-couple(DC)focal mechanism based on the neighbourhood algorithm.One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching.Through synthetic tests using a real passive seismic survey geometry,the results show the new constraint can help better recover the DC components of inverted moment tensors.We further investigate how errors in the velocity model and source location affect the moment tensor solution.The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components.We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests.While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints,both solutions are able to fit the observed waveforms at similar levels.The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.

Strong Earthquake Sequences in Greece during 2008-2014: Moment Tensor Inversions and Fault Plane Discrimination

As is well known, Greece has a significant number of earthquakes each year. Ιn recent years, several earthquakes have occurred in Greece. For this scope, a methodology was used to determine the source parameters. This methodology is based on minimizing the difference between the observed and the synthetic waveforms, using the method Source Parameters Calculation—SPCa [1]. The source parameters, using the proposed methodology, are calculated by comparing observed seismograms and synthetic by inverting data. The synthetics are calculated using the reflectivity method (Kennett, 1983) as implemented by Randall et al. (1994) for a given earth structure. This study includes inversion results for the strongest events that occurred in Greece from 2008 to 2014. For the same events calculated the main fault plane, using the method of Hypocenter Centroid-plot (HC-plot) [2] [3]. This methodology is a simple geometrical method based on the combination between the hypocentral position and the two possible fault planes.

Generalized Inversions of Hadamard and Tensor Products for Matrices

We shall give natural generalized solutions of Hadamard and tensor products equations for matrices by the concept of the Tikhonov regularization combined with the theory of reproducing kernels.

Location and moment tensor inversion of small earthquakes using 3D Green＇s functions in models with rugged topography： application to the Longmenshan fault zone

With dense seismic arrays and advanced imaging methods, regional three-dimensional(3D) Earth models have become more accurate. It is now increasingly feasible and advantageous to use a 3D Earth model to better locate earthquakes and invert their source mechanisms by fitting synthetics to observed waveforms. In this study, we develop an approach to determine both the earthquake location and source mechanism from waveform information. The observed waveforms are filtered in different frequency bands and separated into windows for the individual phases. Instead of picking the arrival times, the traveltime differences are measured by cross-correlation between synthetic waveforms based on the 3D Earth model and observed waveforms. The earthquake location is determined by minimizing the cross-correlation traveltime differences. We then fix the horizontal location of the earthquake and perform a grid search in depth to determine the source mechanism at each point by fitting the synthetic and observed waveforms. This new method is verified by a synthetic test with noise added to the synthetic waveforms and a realistic station distribution. We apply this method to a series of M_W3.4–5.6 earthquakes in the Longmenshan fault(LMSF) zone, a region with rugged topography between the eastern margin of the Tibetan plateau and the western part of the Sichuan basin. The results show that our solutions result in improved waveform fits compared to the source parameters from the catalogs we used and the location can be better constrained than the amplitude-only approach. Furthermore, the source solutions with realistic topography provide a better fit to the observed waveforms than those without the topography, indicating the need to take the topography into account in regions with rugged topography.

Moment tensor inversion of focal mechanism for the aftershock sequence of 1982 Lulong M_S =6.1 earthquake

Using the moment tensor inversion method, we calculate the focal mechanisms of the aftershock sequence of the MS=6.1 Lulong earthquake occurred on October 19, 1982 in Hebei Province. We found that the pressure axis in Lulong basin is nearly in the east-west direction with an azimuth of N74°E. However, in the north of the basin the stress axis changes to N43°E; and in some places near the center of the basin it changes to the northwest that is almost perpendicular to the P axis obtained by us from those events around the basin. This feature illuminates that in Lulong earthquake sequence, the stress direction is different in different parts of crustal structure, which shows that the tectonic movement in Lulong region is complex. This is because that Lulong is located in the eastern part of Chinese mainland and is subject to the compression of Japanese Sea Basin driven by the Pacific Plate. On the other hand, nipped by the Yanshan and North China blocks, Lulong is obviously restricted by the block boundaries.

Moment tensor inversion and source rupture process of the September 27, 2003 M_S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

We conducted moment tensor inversion and studied source rupture process for MS=7.9 earthquake occurred in the border area of China, Russia and Mongolia on September 27 2003, by using digital teleseismic P-wave seismo-grams recorded by long-period seismograph stations of the global seismic network. Considering the aftershock distribution and the tectonic settings around the epicentral area, we propose that the MS=7.9 earthquake occurred on a fault plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of MS=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M0=0.97×1020 N?m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the MS=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of after-shocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.

Research on infrared dim and small target detection algorithm based on low-rank tensor recovery

In order to rapidly and accurately detect infrared small and dim targets in the infrared image of complex scene collected by virtual prototyping of space-based downward-looking multiband detection,an improved detection algorithm of infrared small and dim target is proposed in this paper.Firstly,the original infrared images are changed into a new infrared patch tensor mode through data reconstruction.Then,the infrared small and dim target detection problems are converted to low-rank tensor recovery problems based on tensor nuclear norm in accordance with patch tensor characteristics,and inverse variance weighted entropy is defined for self-adaptive adjustment of sparseness.Finally,the low-rank tensor recovery problem with noise is solved by alternating the direction method to obtain the sparse target image,and the final small target is worked out by a simple partitioning algorithm.The test results in various spacebased downward-looking complex scenes show that such method can restrain complex background well by virtue of rapid arithmetic speed with high detection probability and low false alarm rate.It is a kind of infrared small and dim target detection method with good performance.