The pressure-based acoustic approximation of the elastic wave equations in anisotropic media has advantages corresponding to computational efciency and numerical stability. However, the numerical scattering potentials from the anisotropic parameter perturbations for the pressure wavefeld are not consistent with the elastic scattering theory. In multiparameter full-waveform inversion (FWI), choosing a suitable parameterization, considering the acquisition parameters (e.g., the ofsetto-depth ratio and frequency band) and the accuracy of the anisotropy information in the background initial velocity model, is an important component to a successful anisotropic parameter estimation, because the parameterization determines the trade-of between inverted model parameters and their resolution power. However, because it is difcult to perform multiparameter FWI with various types of parameterization using the pressure-based acoustic wave equation, inaccurate scattered wavefelds cause the gradient direction to lose its unique properties with respect to each model parameter. To overcome these issues, we adopt the combination of pressure- and vector-based acoustic wave equations converted vector virtual sources, which preserves the computational efciency and the angular dependency of the partial derivative wavefelds in elastic media. With the proposed method, we generate the numerical PP scattering patterns for various parameterizations, which are consistent with the elastic scattering theory. Through the numerical tests using the synthetic anisotropic Marmousi-II models and a real ocean bottom cable dataset from the North Sea, we conduct acoustic FWI with three diferent parameterizations using the proposed method and verify that the modifed scattering patterns accurately refect the characteristics of the anisotropic parameter perturbations.
Sejsmiczność indukowana działalnością górniczą jest stale monitorowana i analizowana w celu zwiększenia bezpieczeństwa wydobycia kopalin. Jedną z procedur, chociaż niestosowanych rutynowo w polskim górnictwie węgla i miedzi, jest wyznaczanie mechanizmów ogniskowych silniejszych wstrząsów. Najprostszą i jak dotąd jedyną stosowaną metodą ich obliczania jest inwersja pierwszych wstąpień fali P w domenie czasu. Metoda ta pomimo swojej prostoty i możliwości bezpośredniej kontroli analizowanego wejścia fali sejsmicznej, jest procedurą bardzo czułą na geometrię sieci pomiarowej. Jest to szczególnie kłopotliwe w przypadku dołowych sieci kopalnianych, gdzie geometria sieci ograniczona jest istniejącą infrastrukturą podziemną. W niniejszym opracowaniu proponujemy i testujemy metodę inwersji pełnego pola falowego, przeprowadzaną w domenie częstotliwości oraz czasu, która do tej pory nie była powszechnie stosowana w polskim górnictwie. W oparciu o sejsmogramy syntetyczne obliczone dla lokalnego modelu prędkościowego oraz z wykorzystaniem geometrii rzeczywistej sieci monitoringu sejsmicznego, przeprowadziliśmy badanie przydatności inwersji pełnego pola falowego do obliczeń w warunkach górniczych. W pracy prezentujemy analizy 10 różnych modeli mechanizmów ogniskowych testowanych w różnych lokalizacjach kopalni dla zmiennej siły wstrząsu oraz z dodatkowym rzeczywistym szumem sejsmicznym. Na podstawie otrzymanych wyników możemy stwierdzić, iż głównym czynnikiem warunkującym jakość rozwiązania jest siła wstrząsu i wynikająca z niego zależność amplituda/szum. Natomiast metoda ta nie jest czuła na geometrię sieci pomiarowej.
EN
In order to improve the safety procedures of georesources exploitation, the induced seismicity is constantly monitored and examined. One of the procedures method, although even though not routinely used in Polish coal and copper mining mines, is the calculation of focal mechanisms of strong seismic events. The simplest and so far the only used method for source mechanisms estimation is the inversion of the first P wave onsets in the time domain. Despite its simplicity and the ability to direct control of analyzed seismic waves used during inversion, it is a procedure very sensitive to the geometry of the monitoring seismic network. This is particularly troublesome in the case of underground seismic monitoring system, where the geometry of the network is limited by the existing underground infrastructure. In this study, we propose and test the full waveform inversion method, performed both in the frequency and time domain. This kind of method has not been widely used in Polish miningmines. Basing on synthetic seismograms prepared for the local velocity model and using geometry of the real seismic monitoring network, we conducted a study of the suitability of the full waveform inversion for calculations in underground mining conditions. In this paper, we present analyzes of 10 different source models of focal mechanisms tested in various mine locations within mining area. The tests were performed for variable events strength and with additional real seismic noise. On the grounds of obtained results, we can conclude that the main factor determining the solution quality is the strength of mining tremor and amplitude/noise relationship resulting from it. However, this method is not sensitive to the geometry of the measurement network.
3
Dostęp do pełnego tekstu na zewnętrznej witrynie WWW
The conventional full-waveform inversion (FWI) often minimizes the objective function using some local optimization algorithms. As a result, when the initial model is not good enough, the inversion process will drop into a local minimum. The low-frequency components contained in seismic data are of vital importance for reducing the initial model dependence and mitigating the cycle-skipping phenomenon of FWI. In this research, a frequency extension method using the nth power operation is proposed, which compresses the seismic data in time domain and extends their frequency band. Based on this, we construct a new objective function using the nth power wavefeld and derive the corresponding gradient formula. The new objective function shows better property to overcome local minimum than the conventional one. When conduct inversion, we can invert from high-order to low-order successively, which is a new multiscale strategy. Since seismic data is more sensitive to source wavelet errors after high-order operation, we make the method more robust by proposing a source-independent method to mitigate the efects of source wavelet inaccuracy. After that, we extend the proposed method to encoded multisource waveform inversion. The numerical examples on the Marmousi model demonstrate that the proposed method can efectively mitigate the cycle-skipping of FWI, and it also has good anti-noise property.
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.