Imaging quasi-vertical geological faults with earthquake data

P. Zheglova, J. R. Mclaughlin, S. W. Roecker, J. R. Yoon, Paul Renzi

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

We present a method for imaging quasi-vertically dipping faults with surface records of reflected P waves from small earthquakes. Faults are boundaries between geological structures, such as tectonic plates, and are located in earthquake active regions such as Parkfield, California. The high degree of seismic activity enables the use of multiple seismic recordings in our fault identification algorithm. Major challenges occur because of the quasi-vertical orientation of the fault and the fact that the wave reflected by the fault and recorded by the surface receivers is not well modelled by the direct arrival of the propagating wave generated by the earthquake source. Our method uses the 2-D acoustic wave equation as the model for P-wave propagation. We assume that an approximate wave speed map on the reflection side of the fault is available and the source locations are known, for example, from traveltime tomography. We also assume that the source time function is known. The new features of our method arise because earthquake sources are located very close to the fault. This has two implications: (1) the direct arrival and the reflected wave arrive almost simultaneously, so that it is impossible to separate them on a seismogram using standard techniques, and (2) most of the reflections occur above the critical angle which introduces a distortion in the reflected wave. To overcome these difficulties we use a modelled incident wave to (1) remove the direct arrival from the data, and (2) remove the post-critical distortion from the reflected wave. We justify the distortion removal using the leading-order term of an asymptotic expansion, and an optimization procedure. To complete our algorithm we utilize some features of reverse time migration: (1) the use of full acoustic wave equation for modelling and backpropagation, and (2) zero-lag correlation of the backpropagated time reversed reflected and incident fields. We present numerical examples of fault reconstructions with synthetic data.

Original languageEnglish
Pages (from-to)1584-1596
Number of pages13
JournalGeophysical Journal International
Volume189
Issue number3
DOIs
Publication statusPublished - 1 Jun 2012

Fingerprint

reflected waves
Earthquakes
earthquakes
arrivals
Imaging techniques
earthquake
P waves
wave equations
vertical orientation
plates (tectonics)
time functions
acoustics
seismograms
acoustic wave
wave equation
Wave equations
dipping
reaction time
P-wave
wave propagation

Keywords

  • Body waves
  • Wave propagation
  • Wave scattering and diffraction

ASJC Scopus subject areas

  • Geophysics
  • Geochemistry and Petrology

Cite this

Imaging quasi-vertical geological faults with earthquake data. / Zheglova, P.; Mclaughlin, J. R.; Roecker, S. W.; Yoon, J. R.; Renzi, Paul.

In: Geophysical Journal International, Vol. 189, No. 3, 01.06.2012, p. 1584-1596.

Research output: Contribution to journalArticle

Zheglova, P. ; Mclaughlin, J. R. ; Roecker, S. W. ; Yoon, J. R. ; Renzi, Paul. / Imaging quasi-vertical geological faults with earthquake data. In: Geophysical Journal International. 2012 ; Vol. 189, No. 3. pp. 1584-1596.
@article{04c79f5f5bd44c3c8203a28415b59c15,
title = "Imaging quasi-vertical geological faults with earthquake data",
abstract = "We present a method for imaging quasi-vertically dipping faults with surface records of reflected P waves from small earthquakes. Faults are boundaries between geological structures, such as tectonic plates, and are located in earthquake active regions such as Parkfield, California. The high degree of seismic activity enables the use of multiple seismic recordings in our fault identification algorithm. Major challenges occur because of the quasi-vertical orientation of the fault and the fact that the wave reflected by the fault and recorded by the surface receivers is not well modelled by the direct arrival of the propagating wave generated by the earthquake source. Our method uses the 2-D acoustic wave equation as the model for P-wave propagation. We assume that an approximate wave speed map on the reflection side of the fault is available and the source locations are known, for example, from traveltime tomography. We also assume that the source time function is known. The new features of our method arise because earthquake sources are located very close to the fault. This has two implications: (1) the direct arrival and the reflected wave arrive almost simultaneously, so that it is impossible to separate them on a seismogram using standard techniques, and (2) most of the reflections occur above the critical angle which introduces a distortion in the reflected wave. To overcome these difficulties we use a modelled incident wave to (1) remove the direct arrival from the data, and (2) remove the post-critical distortion from the reflected wave. We justify the distortion removal using the leading-order term of an asymptotic expansion, and an optimization procedure. To complete our algorithm we utilize some features of reverse time migration: (1) the use of full acoustic wave equation for modelling and backpropagation, and (2) zero-lag correlation of the backpropagated time reversed reflected and incident fields. We present numerical examples of fault reconstructions with synthetic data.",
keywords = "Body waves, Wave propagation, Wave scattering and diffraction",
author = "P. Zheglova and Mclaughlin, {J. R.} and Roecker, {S. W.} and Yoon, {J. R.} and Paul Renzi",
year = "2012",
month = "6",
day = "1",
doi = "10.1111/j.1365-246X.2012.05424.x",
language = "English",
volume = "189",
pages = "1584--1596",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "3",

}

TY - JOUR

T1 - Imaging quasi-vertical geological faults with earthquake data

AU - Zheglova, P.

AU - Mclaughlin, J. R.

AU - Roecker, S. W.

AU - Yoon, J. R.

AU - Renzi, Paul

PY - 2012/6/1

Y1 - 2012/6/1

N2 - We present a method for imaging quasi-vertically dipping faults with surface records of reflected P waves from small earthquakes. Faults are boundaries between geological structures, such as tectonic plates, and are located in earthquake active regions such as Parkfield, California. The high degree of seismic activity enables the use of multiple seismic recordings in our fault identification algorithm. Major challenges occur because of the quasi-vertical orientation of the fault and the fact that the wave reflected by the fault and recorded by the surface receivers is not well modelled by the direct arrival of the propagating wave generated by the earthquake source. Our method uses the 2-D acoustic wave equation as the model for P-wave propagation. We assume that an approximate wave speed map on the reflection side of the fault is available and the source locations are known, for example, from traveltime tomography. We also assume that the source time function is known. The new features of our method arise because earthquake sources are located very close to the fault. This has two implications: (1) the direct arrival and the reflected wave arrive almost simultaneously, so that it is impossible to separate them on a seismogram using standard techniques, and (2) most of the reflections occur above the critical angle which introduces a distortion in the reflected wave. To overcome these difficulties we use a modelled incident wave to (1) remove the direct arrival from the data, and (2) remove the post-critical distortion from the reflected wave. We justify the distortion removal using the leading-order term of an asymptotic expansion, and an optimization procedure. To complete our algorithm we utilize some features of reverse time migration: (1) the use of full acoustic wave equation for modelling and backpropagation, and (2) zero-lag correlation of the backpropagated time reversed reflected and incident fields. We present numerical examples of fault reconstructions with synthetic data.

AB - We present a method for imaging quasi-vertically dipping faults with surface records of reflected P waves from small earthquakes. Faults are boundaries between geological structures, such as tectonic plates, and are located in earthquake active regions such as Parkfield, California. The high degree of seismic activity enables the use of multiple seismic recordings in our fault identification algorithm. Major challenges occur because of the quasi-vertical orientation of the fault and the fact that the wave reflected by the fault and recorded by the surface receivers is not well modelled by the direct arrival of the propagating wave generated by the earthquake source. Our method uses the 2-D acoustic wave equation as the model for P-wave propagation. We assume that an approximate wave speed map on the reflection side of the fault is available and the source locations are known, for example, from traveltime tomography. We also assume that the source time function is known. The new features of our method arise because earthquake sources are located very close to the fault. This has two implications: (1) the direct arrival and the reflected wave arrive almost simultaneously, so that it is impossible to separate them on a seismogram using standard techniques, and (2) most of the reflections occur above the critical angle which introduces a distortion in the reflected wave. To overcome these difficulties we use a modelled incident wave to (1) remove the direct arrival from the data, and (2) remove the post-critical distortion from the reflected wave. We justify the distortion removal using the leading-order term of an asymptotic expansion, and an optimization procedure. To complete our algorithm we utilize some features of reverse time migration: (1) the use of full acoustic wave equation for modelling and backpropagation, and (2) zero-lag correlation of the backpropagated time reversed reflected and incident fields. We present numerical examples of fault reconstructions with synthetic data.

KW - Body waves

KW - Wave propagation

KW - Wave scattering and diffraction

UR - http://www.scopus.com/inward/record.url?scp=84860837514&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84860837514&partnerID=8YFLogxK

U2 - 10.1111/j.1365-246X.2012.05424.x

DO - 10.1111/j.1365-246X.2012.05424.x

M3 - Article

AN - SCOPUS:84860837514

VL - 189

SP - 1584

EP - 1596

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 3

ER -