Repeatability analysis of land time‐lapse seismic data: CO2CRC Otway pilot project case study

Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO₂ sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO₂ sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO₂ sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO₂ sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.     

Frequency domain modal analysis of earthquake input energy to highly damped passive control structures

A new complex modal analysis‐based method is developed in the frequency domain for efficient computation of the earthquake input energy to a highly damped linear elastic passive control structure. The input energy to the structure during an earthquake is an important measure of seismic demand. Because of generality and applicability to non‐linear structures, the earthquake input energy has usually been computed in the time domain. It is shown here that the formulation of the earthquake input energy in the frequency domain is essential for deriving a bound on the earthquake input energy for a class of ground motions and for understanding the robustness of passively controlled structures to disturbances with various frequency contents. From the viewpoint of computational efficiency, a modal analysis‐based method is developed. The importance of overdamped modes in the energy computation of specific non‐proportionally damped models is demonstrated by comparing the energy transfer functions and the displacement transfer functions. Through numerical examinations for four recorded ground motions, it is shown that the modal analysis‐based method in the frequency domain is very efficient in the computation of the earthquake input energy. Copyright © 2004 John Wiley & Sons, Ltd.     

A model‐based data‐driven dictionary learning for seismic data representation

Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real world, seismic events usually have curvature or amplitude variability, which means that their Fourier transforms are no longer strictly linear but rather occupy conic regions of the Fourier domain that are narrow at low frequencies but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised “signal cones”. In this work, we consider a space–time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting, and translation of the kernel for cone‐limited (C‐limited) functions (functions whose Fourier transform lives within a cone) or C‐Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency), which constitutes our dictionary. We find a discrete number of scaling, slanting, and translation parameters from a continuum by optimally matching the data. This is a non‐linear optimisation problem, which we address by a fixed‐point method that utilises a variable projection method with ?₁ constraints on the linear parameters and bound constraints on the non‐linear parameters. We observe that slow decay and oscillatory behaviour of the kernel for C‐limited functions constitute bottlenecks for the optimisation problem, which we partially overcome by the C‐Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity‐promoting curvelet decomposition, matching pursuit Fourier interpolation, and sparsity‐promoting plane‐wave decomposition methods.     

Growth rate is temperature‐dependent,but the proportion of resource allocation between the growth and gonads is not,in early adult stage of the lancelet branchiostoma japonicum

The pattern of resource allocation into multiple life‐history traits can vary with the environment; for example, high temperatures generally facilitate growth. Changes in investment in growth due to temperature may alter the allocation to other traits that are in trade‐off with that. Here, we investigated the trade‐off between growth and reproduction and whether the proportion of allocation to these varies with increasing seawater temperature in the lancelet Branchiostoma japonicum. We reared the lancelets of early adult stage under three different temperature treatments (control, + 1°C and + 2°C treatments) for 427 days in the laboratory. Although the total body length at initiation did not differ between treatments, the length was influenced by the interaction between elapsed days and temperature treatments, suggesting that growth rate was temperature‐dependent; the lancelets under higher temperature conditions showed a more rapid growth rate than did those under the control conditions. The frequency of individuals with gonads did not differ between treatments throughout their reproductive season. The monthly growth increments were also not significantly affected by the interaction between the presence/absence of gonads and temperature treatments, although the presence of gonads caused the growth increments to decrease. In addition, the effect of temperature treatments was significant. These results suggest that changes in temperature did not alter the pattern of allocation between growth and reproduction that were in a trade‐off, whereas high temperature would solely result in an increase in resource intake, at least in the early adult stage of this species.     

Variation of mineral composition,fabric and oxygen fugacity from massive to foliated eclogites during exhumation of subducted ocean crust in the North Qilian suture zone,NW China

Eclogites from the North Qilian suture zone are high‐pressure low‐temperature metamorphic rocks of ocean crust protolith, and occur in both massive and foliated varieties as individual blocks of tens to hundreds of metres in size. The massive type is weakly deformed and shows granoblastic texture characterized by a coarse‐grained peak mineral assemblage of Grt¹ + Omp¹ + Ph + Rt ± Lws (or retrograde Cz). In contrast, the foliated type is strongly deformed and shows a fine‐grained retrograde mineral assemblage of Grt² + Omp² + Cz + Gln + Ph. Both total FeO and aegirine contents in omphacite, as well as X_Fe[=Fe³⁺/(Fe³⁺ + Al^VI)] in clinozoisite/epidote, increase significantly from massive to foliated eclogites. Lattice preferred orientation (LPO) of omphacite, determined by electron back‐scatter diffraction analysis, is characterized by weak and strong SL‐type fabrics for massive and foliated eclogites, respectively. Clinozoisite/epidote also developed SL‐type fabric, but different from the LPOs of omphacite in <010> and <001> axes, owing to their opposite crystallographic long and short axis definitions. The transition of deformation mechanism from dislocation creep to diffusive mass transfer (DMT) creep in omphacite and the concomitant retrograde metamorphism both are efficiently facilitated when the original coarse‐grained Omp¹ + Grt¹ + Lws assemblage is dynamically recrystallized and retrogressed into the fine‐grained Fe³⁺‐rich assemblage of Omp² + Grt² + Cz + Gln. The DMT process with concomitant anisotropic growth assisted by fluids is considered to be an important deformation mechanism for most minerals in the foliated eclogite. P–T estimates yielded 2.3–2.6 GPa and 485?510 °C for the massive eclogite and 1.8–2.2 GPa and 450?480 °C for the foliated eclogite. The significant increase in total Fe and Fe³⁺ contents in omphacite and clinozoisite/epidote from massive to foliated eclogite suggests changes in mineral compositions accompanied by an increase in oxygen fugacity during ductile deformation associated with exhumation. The LPO transition of omphacite, clinozoisite and rutile from weak SL‐type in massive eclogites to strong SL‐type in foliated eclogites is interpreted to represent the increment of shear strain during exhumation along the ‘subduction channel’.     

The Characteristics and Causes of Land Subsidence in Tanggu Based on the GPS Survey System and Numerical Simulation

Land subsidence is a severe hazard threatening Tanggu, a flat lowland area, and evidences of land subsidence can be seen throughout the city. A new reasonable GPS network was set up in this area from 2008 to 2010. The monitoring data show that land subsidence was serious and two main subsidence cones were obviously formed in the region. One emerged at Hujiayuan, with the maximum rate reaching 60 mm/a, and the influence region enlarged prominently from 2005 to 2010. The other one occurred at Kaifaqu, which became obvious only after 2005, and it showed a decreasing tendency with time. To analyze the causes of ground settlement, a correlation between groundwater withdrawal and land subsidence was firstly made. The results confirmed that over-exploitation of groundwater was the major cause for the severe settlement in Hujiayuan. Meanwhile, the subsidence of Kaifaqu was also related to groundwater withdrawal before 2005. However, the relationship became unconspicuous after 2005. To find the cause of this abnormity, a three-dimensional finite element numerical model, coupled with groundwater flow and subsidence, was built. The simulation results indicate that the subsidence induced by high-rise buildings is serious, but the affected range is limited and it also shows a decreasing trend with time, corresponding to the subsidence characteristics at Kaifaqu. Therefore, more attention should be paid to this hazard induced by engineering construction besides groundwater withdrawal, as more high-rise buildings are under construction in Tanggu.     

瞬态瑞雷波法在工程检测中的应用研究

本文简要介绍了瞬态瑞雷波法的工作原理,野外观测及数据采集系统,并总结出了资料处理解释流程图,同时从工程实践出发,探讨了该项新技术应用于高速公路建设中的优点,存在问题及未来发展趋势。