Browsing by Author "Tao, Kai"

Now showing 1 - 5 of 5
  • Thumbnail Image
    Item
    Estimating sedimentary and crustal structure using wavefield continuation: theory, techniques and applications
    (Oxford University Press, 2014) Tao, Kai; Liu, Tianze; Ning, Jieyuan; Niu, Fenglin
    Receiver function techniques are widely used in imaging crustal and mantle structure beneath a seismic station. The weak P-to-S conversions at deep seismic structures are usually masked by strong shallow reverberations when unconsolidated sediments are present below the station, making it nearly impossible to utilize receiver function techniques. We develop a method to estimate sediment and crustal structures beneath a seismic station based on wavefield downward continuation and decomposition method. The method parametrizes velocity structure beneath the station with a stack of constant velocity layers overlying a homogeneous half-space, and approximates the teleseismic P wave and its coda by the structural response to an incoming plane P wave. Our method is based on the principle that the upgoing S wavefield is absent in the half-space, and searches for the optimum velocity and thickness of the layers that give the minimum S-wave energy flux from the half-space to the layers. An iterative grid-search algorithm from the top to the bottom layers is employed to implement the search. In this study, we only use models comprising either only one crust layer or two layers (sediment+crust) with a half-space mantle, although models with more layers are also implementable. The method is especially useful in resolving seismic structure beneath a station sitting on unconsolidated sediments. It not only can be used to determine the sediment thickness and velocity structure, but also provides an effective way to generate subsurface receiver functions, which are formed by deconvolving the upgoing P wavefield from the upgoing S waves at the top of hardrock crust, and thus are free from shallow reverberations. The technique is applied to various synthetic data generated with different types of velocity model and noise levels, and appears to have good capability in recovering the input models. We further applied this method to teleseismic data recorded at a station inside the Songliao Basin in northeast China. The estimated sediment thickness and velocity agrees well with the results of previous activesource studies. The subsurface receiver functions also show a superior power in exposing the Moho Ps conversions, resulting in a well-defined peak in the H-? domain, which are absent in the regular receiver function data.
  • Thumbnail Image
    Item
    Full-waveform inversion of triplicated data using a normalized-correlation-coefficient-based misfit function
    (Oxford University Press, 2017) Tao, Kai; Grand, Stephen P.; Niu, Fenglin
    In seismic full-waveform inversion (FWI), the choice of misfit function determines what information in data is used and ultimately affects the resolution of the inverted images of the Earth's structure. Misfit functions based on traveltime have been successfully applied in global and regional tomographic studies. However, wave propagation through the upper mantle results in multiple phases arriving at a given receiver in a narrow time interval resulting in complicated waveforms that evolve with distance. To extract waveform information as well as traveltime, we use a misfit function based on the normalized correlation coefficient (CC). This misfit function is able to capture the waveform complexities in both phase and relative amplitude within the measurement window. It is also insensitive to absolute amplitude differences between modeled and recorded data, which avoids problems due to uncertainties in source magnitude, radiation pattern, receiver site effects or even miscalibrated instruments. These features make the misfit function based on normalized CC a good candidate to achieve high-resolution images of complex geological structures when interfering phases coexist in the measurement window, such as triplication waveforms. From synthetic tests, we show the advantages of this misfit function over the cross-correlation traveltime misfit function. Preliminary inversion of data from an earthquake in Northeast China images a sharper and stronger amplitude slab stagnant in the middle of the transition zone than FWI of cross-correlation traveltime.
  • Thumbnail Image
    Item
    Hydraulic Injection‐Induced Velocity Changes Revealed by Surface Wave Coda and Polarization Data at a Shale Play Site in Southwest China
    (Wiley, 2020) Zhang, Yan; Niu, Fenglin; Tao, Kai; Ning, Jieyuan; Chen, Haichao; Tang, Youcai
    We investigated temporal variations of seismic wave velocity associated with hydraulic fracturing using Green's functions computed from ambient noise data. In October and November of 2014, we set up a broadband array at a shale play site inside the Sichuan basin where a pilot horizontal drilling and hydraulic injections were conducted. We first computed cross‐correlation functions using continuous data recorded by 21 three‐component broadband sensors deployed around the treatment well. We then employed a running window correlation‐based coda wave interferometry technique to measure apparent velocity changes from the daily Green's functions of all the station pairs in the frequency range of 1 to 3 Hz. We found significant velocity changes right after the hydraulic fracturing, which exhibited a clear direction‐dependent pattern. S wave velocity along raypaths parallel to the well trajectory showed a clear increase while those perpendicular exhibited a small decrease. The anisotropic changes in seismic velocity observed here were also confirmed from surface wave horizontal particle motion data. By comparing our observations with normal stress changes calculated with a half‐space elastic model, we speculate that stress changes induced by the hydraulic fracturing were likely to be responsible for the observed anisotropic changes in seismic velocity. Our results suggest that time‐lapse seismic imaging with ambient noise data provides a promising probe for monitoring geomechanical changes related to exploitation of unconventional oil and gas resources.
  • Thumbnail Image
    Item
    Seismic Structure of the Upper Mantle Beneath Eastern Asia From Full Waveform Seismic Tomography
    (Wiley, 2018) Tao, Kai; Grand, Stephen P.; Niu, Fenglin
    To better understand the subsurface behavior of subducting slabs and their relation to the tectonic evolution of the overriding plate, we conduct a full waveform inversion on a large data set to determine a high‐resolution seismic model, FWEA18 (Full Waveform inversion of East Asia in 2018), of the upper mantle beneath eastern Asia. FWEA18 reveals sharper, more intense high‐velocity slabs in the upper mantle under the southern Kuril, Japan, and Ryukyu arcs, than previous studies have found. The subducting Pacific plate is imaged as a roughly 100 km thick high‐velocity slab to near 550 km depth indicating relatively little deformation. Stagnation near 600 km depth is observed over horizontal distances of 600 km or less. The Pacific plate we image accounts for roughly 25 Myr of subduction with older slab likely located in the lower mantle. The Philippine plate, subducting beneath the Ryukyu Islands, has a clear termination at about 450 km depth. This may indicate a tearing event in the past or that less Philippine Sea plate has subducted than previously thought. We found a double‐layer high‐velocity anomaly above and below 660 km under the Yellow Sea and eastern coast of North China. This may correspond to parts of the Philippine Sea plate that detached in the past and Pacific plate that have intersected at depth or a complicated behavior of the Pacific plate at that depth. Slow cylindrical anomalies cross the entire upper mantle are imaged beneath major Holocene volcanoes, which are likely upwellings associated with the edges of deep slabs that are entering the lower mantle.
  • Thumbnail Image
    Item
    Temporal Variations of Near‐Surface Anisotropy Induced by Hydraulic Fracturing at a Shale Play Site in Southwest China
    (Wiley, 2018) Zuo, Qiankun; Tang, Youcai; Niu, Fenglin; Li, Guoliang; Chen, Haichao; Tao, Kai; Chen, Benchi
    Knowledge of the geometric properties of fractures and cracks in a petroleum reservoir is important to reservoir exploitation. When aligned and partially connected, fractures and cracks can act as conduits for fluid flow and thus can significantly increase the permeability of the reservoir. The aligned fractures and cracks, on the other hand, are an effective means to generate seismic anisotropy. In this study, we utilize the seismic data recorded by a vertical array installed in a shallow borehole at a shale play site in southwest China. By applying seismic interferometry to the ambient noise data recorded by 12 three‐component geophones, we extract P and S waves propagating vertically along the borehole. The S waves show up to 20% velocity variations with respect to their polarization directions. Such large S wave anisotropy can be explained by the horizontal transverse isotropic model and is likely caused by natural fractures that are widely present in the area and align approximately in the NE‐SW direction. During the 13‐day period of hydraulic fracking treatment, we also observe large and systematic temporal variations in S wave velocity, degree of S wave polarization anisotropy, and fast polarization direction. By comparing our observations with normal strain changes calculated with a half‐space elastic model, we speculate that strain changes induced by hydraulic injection and fracturing are likely to be responsible for the observed temporal variations in seismic anisotropy. As such, seismic interferometry with shallow borehole acquisition might provide an alternative means to monitor hydraulic fracturing and wastewater injection in the future.