Automatic cross-well tomography: an application of the differential semblance optimization to two real examples

An automatic tomography algorithm, based on differential semblance optimization (DSO), has been used to invert real cross-well seismic data for the background velocity. The method relies on the first-arrival transmitted waves. Given a background velocity model, the traveltimes between the sources and the receivers are computed, then semblance panels are created by back-propagating the data traces. If the velocity model is correct all the first-arrival transmitted waves will be aligned in the semblance panels. The DSO method consists of finding the background velocity by minimizing the L ₂-norm of the difference between adjacent back-propagated traces. Thanks to the good behaviour of this DSO cost function about the solution, a local (gradient) optimization can be performed. This provides a relatively fast algorithm when ray tracing and analytic computation of the gradient are used.
Unfortunately the method fails in the presence of caustics in the data. However, this difficulty can be circumvented by applying suitable masks to the data. This approach is first applied to a synthetic example then to two real data sets: the McElroy data set recorded in West Texas and the NIMR data set recorded in Oman. The results are quite encouraging and similar to those obtained with classical tomography.     

Inversion of reflection seismograms by differential semblance analysis: algorithm structure and synthetic examples1

Seismograms predicted from acoustic or elastic earth models depend very non-linearly on the long wavelength components of velocity. This sensitive dependence demands the use of special variational principles in waveform-based inversion algorithms. The differential semblance variational principle is well-suited to velocity inversion by gradient methods, since its objective function is smooth and convex over a large range of velocity models. An extension of the adjoint state technique yields an accurate estimate of the differential semblance gradient. Non-linear conjugate gradient iteration is quite successful in locating the global differential semblance minimum, which is near the ordinary least-squares global minimum when coherent data noise is small. Several examples, based on the 2D primaries-only acoustic model, illustrate features of the method and its performance.     

A theory of atmospheric tides generated by the differential heating between land and sea

Summary The differential heating between land and sea is incorporated into the theory of atmospheric tides. This involves the representation of the land and sea distribution by a set of Fourier series.The theory postulates the existence of waves of angular speeds different from the angular velocity of the earth with some of the waves travelling from West to East instead of the usual East to West.By considering the diurnal variation of eddy diffusion of heat energy absorbed close to the surface, the semidiurnal standing waves at the poles was calculated by the application of this theory. The order of magnitude of the calculated tides agreed well with observation, though, contrary to observation, the calculated amplitude at the South pole is larger than that at the North pole.It is also shown that the polar standing oscillation is caused mainly by the land and sea distribution between 75°N and 45°S.     

Solute transport by surface runoff from low‐angle slopes: theory and application

The removal of chemicals in solution by overland flow from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface‐mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright © 2000 John Wiley & Sons, Ltd.     

约束优化的微分进化算法在波动方程反问题中的应用

本文提出了一种约束优化的微分进化算法,该算法使得种群在演化过程中能保持较好的多样性,且参数设置简单,不容易陷入局部最优,并能在较短时间内找到约束优化问题的最优解.结合有限元法对波动方程反问题进行研究,在对一维及二维波动方程反问题的数值模拟中都得到了较好的结果,表明了该算法具有较好的稳定性和较强的抗噪能力.     

Automatic,unstructured mesh optimization for simulation and assessment of tide- and surge-driven hydrodynamics in a longitudinal estuary: St. Johns River

A localized truncation error analysis with complex derivatives (LTEA+CD) is applied recursively with advanced circulation (ADCIRC) simulations of tides and storm surge for finite element mesh optimization. Mesh optimization is demonstrated with two iterations of LTEA+CD for tidal simulation in the lower 200 km of the St. Johns River, located in northeast Florida, and achieves more than an over 50% decrease in the number of mesh nodes, relating to a twofold increase in efficiency, at a zero cost to model accuracy. The recursively generated meshes using LTEA+CD lead to successive reductions in the global cumulative truncation error associated with the model mesh. Tides are simulated with root mean square error (RMSE) of 0.09–0.21 m and index of agreement (IA) values generally in the 80s and 90s percentage ranges. Tidal currents are simulated with RMSE of 0.09–0.23 m s^?1 and IA values of 97% and greater. Storm tide due to Hurricane Matthew 2016 is simulated with RMSE of 0.09–0.33 m and IA values of 75–96%. Analysis of the LTEA+CD results shows the M2 constituent to dominate the node spacing requirement in the St. Johns River, with the M4 and M6 overtides and the STEADY constituent contributing some. Friction is the predominant physical factor influencing the target element size distribution, especially along the main river stem, while frequency (inertia) and Coriolis (rotation) are supplementary contributing factors. The combination of interior- and boundary-type computational molecules, providing near-full coverage of the model domain, renders LTEA+CD an attractive mesh generation/optimization tool for complex coastal and estuarine domains. The mesh optimization procedure using LTEA+CD is automatic and extensible to other finite element-based numerical models. Discussion is provided on the scope of LTEA+CD, the starting point (mesh) of the procedure, the user-specified scaling of the LTEA+CD results, and the iteration (termination) of LTEA+CD for mesh optimization.     

Effects of petrophysical parameters on attenuation and dispersion of seismic waves in the simplified poroelastic theory: Comparison between the Biot's theory and viscosity‐extended theory

The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.     

Sonification of geophysical data through time–frequency analysis: theory and applications

In this paper, we introduce a new method of geophysical data interpretation based on simultaneous analysis of images and sounds. The final objective is to expand the interpretation workflow through multimodal (visual–audio) perception of the same information. We show how seismic data can be effectively converted into standard formats commonly used in digital music. This conversion of geophysical data into the musical domain can be done by applying appropriate time–frequency transforms. Using real data, we demonstrate that the Stockwell transform provides a very accurate and reliable conversion. Once converted into musical files, geophysical datasets can be played and interpreted by using modern computer music tools, such as sequencers. This approach is complementary and not substitutive of interpretation methods based on imaging. It can be applied not only to seismic data but also to well logs and any type of geophysical time/depth series. To show the practical implications of our integrated visual–audio method of interpretation, we discuss an application to a real seismic dataset in correspondence of an important hydrocarbon discovery.     

Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography

We report a new model of the upper mantle structure beneath Italy obtained by means of P-wave teleseismic tomography. Besides the recent and remarkable development of the Italian Seismic Network, a high model resolution has been achieved improving the inversion method upon the ACH method used in previous investigations and picking high quality arrival times with the Multi-Channel Cross-Correlation technique. The finer details of our Vp model yield new insights into the heterogeneous structure of the Adria continental lithosphere involved in the collision between the Africa and Europe plates. A wide low Vp anomaly located in the northern Adria mantle, facing the Alpine high Vp slab, supports the idea that the Adria lithosphere has been hydrated and thinned during the Alpine subduction. We argue that this mantle softening may have played a key role in favoring the subsequent delamination of the Adria lithosphere in the northern Apennines. We hypothesize that delamination of continental lithosphere previously thinned in a back-arc setting may be considered a key process to favor subduction polarity reversal and recycling of continental material into the mantle circulation. Conversely, in the central-southern Apennines, the velocity structure is consistent with the existence of a deeper oceanic slab that flattens at the base of the upper mantle, in agreement with the widely accepted geodynamic evolution of the central Mediterranean by slab retreat and back-arc spreading. The oceanic slab is discontinuously detached from the surface plate, suggesting a different structure of the Adria lithosphere, which resists subduction instead of favoring delamination.     

Subspace system identification of support‐excited structures—part I: theory and black‐box system identification

This paper reviews the theoretical principles of subspace system identification as applied to the problem of estimating black‐box state‐space models of support‐excited structures (e.g., structures exposed to earthquakes). The work distinguishes itself from past studies by providing readers with a powerful geometric interpretation of subspace operations that relates directly to theoretical structural dynamics. To validate the performance of subspace system identification, a series of experiments are conducted on a multistory steel frame structure exposed to moderate seismic ground motions; structural response data is used off‐line to estimate black‐box state‐space models. Ground motions and structural response measurements are used by the subspace system identification method to derive a complete input–output state‐space model of the steel frame system. The modal parameters of the structure are extracted from the estimated input–output state‐space model. With the use of only structural response data, output‐only state‐space models of the system are also estimated by subspace system identification. The paper concludes with a comparison study of the modal parameters extracted from the input–output and output‐only state‐space models in order to quantify the uncertainties present in modal parameters extracted from output‐only models. Copyright © 2012 John Wiley & Sons, Ltd.     

On the ecohydrology of the Yucatan Peninsula: Evapotranspiration and carbon intake dynamics across an eco‐climatic gradient

The hydrology and productivity of the ecosystems of the Yucatan Peninsula (YP) are highly constrained by two factors: (a) the lack of surface drainage networks due to the existence of a highly permeable and connected karstic aquifer roughly the size of the peninsula and (b) a climatic gradient that leads to a transition from seasonally dry deciduous and sub‐deciduous tropical forests, in the north‐western and central parts of the Peninsula, to evergreen forests, in the southern and eastern parts. As a result, surface water fluxes of the YP are restricted to evapotranspiration (ET) that are tightly coupled to ecosystems health and gross primary productivity (GPP). The magnitude and seasonal variation of these fluxes are sensitive to climatic variability and perturbations caused by extreme events such as droughts and tropical storms that are frequent in the YP. In this study, we assess the spatio‐temporal dynamics of ET and GPP above average dry and wet conditions through time series analyses of 15 years of remotely sensed data from both Moderate Resolution Imaging Spectroradiometer and Tropical Rainfall Measuring Mission satellite products. Our results show that ET and GPP follow a regional moisture and temperature gradient that highly controls the distribution of ecosystems within the peninsula. We observe that ET and GPP are in phase with the rainy season in the deciduous forests, but for the evergreen forests, only the GPP is in phase. Additionally, and with the exception of droughts on deciduous ecosystems of the northern part of the YP, the productivity of these ecosystems shows a legacy effect, responding more to a defined trajectory (wetting or drying on the previous years), rather than to punctual extreme climatic events. This has implications on the resilience of these ecosystems to natural perturbations of climate. Comparisons between deciduous and evergreen forest indicate that both types of ecosystems have different plant water use strategies in response to hydrologic variability.     

Mechanical compaction in heterogeneous clastic formations from plastic–poroelastic deformation principles: theory and applications

Mechanical compaction or loss of porosity due to increase in effective stress is a fundamental geological process that governs many of the rock elastic and transport parameters, all of great importance in exploring and developing subsurface reservoirs. The ability to model the compaction process enables us to improve our understanding of the seismic signature of the basin and better relate the geology of deposition to current porosity, velocity, pore pressure, and other mechanical parameters that depend on the state of compaction of the sediment. In this paper, a set of mathematical equations that can be used to model the plastic deformation associated with primary and secondary loading curves is presented. Compaction laws are posed in terms of natural strain increment formulation often used in plasticity theory to model large deformation. Laboratory and field estimates of constitutive plastic deformation relations for sand–shale mixtures are used in a numerical model that generates estimates of porosity under various pore pressures, shale content, and loading scenarios. These estimates can be used in a variety of settings to predict various basin and reservoir properties associated with different loading conditions and/or sedimentation processes.     

Quantitative imaging of the Permo‐Mesozoic complex and its basement by frequency domain waveform tomography of wide‐aperture seismic data from the Polish Basin

In this study we present the workflow and results of 2D frequency domain waveform tomography applied to the global‐offset seismic data acquired in central Poland along a 50‐km long profile during the GRUNDY 2003 experiment. The waveform tomography method allows full exploitation of the wide‐aperture content of these data and produces in a semi‐automatic way both the detailed P‐wave velocity model and the structural image (i.e., perturbations in respect to the starting model). Thirteen frequencies ranging from 4 to 16 Hz were inverted sequentially, gradually introducing higher wavenumbers and more details into the velocity models. Although the data were characterised by relatively large shot spacings (1.5 km), we obtained clear images both of the Mesozoic and Permian sedimentary cover. Velocity patterns indicated facies changes within the Jurassic and Zechstein strata. A high velocity layer (ca. 5500 m/s) was found near the base of Triassic (Scythian), which made the imaging of a deeper layer difficult. Nevertheless, we were able to delineate the base of the Permian (i.e., base of the Rotliegend), which was not possible to derive from conventional common‐depth‐point processing, as well as some deeper events, attributed to the Carboniferous. The sub‐Permian events formed a syn‐form which favoured our previous interpretation of a depression filled with Upper Carboniferous molasse. The validity of the waveform tomography‐derived model was confirmed by well‐log data. Forward ray‐tracing modelling and synthetic seismograms calculations provided another justification for the key structures present in the waveform tomography model.     

Concentration profiles for fine and coarse sediments suspended by waves over ripples: An analytical study with the 1-DV gradient diffusion model

Field and laboratory measurements of suspended sediments over wave ripples show, for time-averaged concentration profiles in semi-log plots, a contrast between upward convex profiles for fine sand and upward concave profiles for coarse sand. Careful examination of experimental data for coarse sand shows a near-bed upward convex profile beneath the main upward concave profile. Available models fail to predict these two profiles for coarse sediments. The 1-DV gradient diffusion model predicts the main upward concave profile for coarse sediments thanks to a suitable β(y)$\beta (y)$-function (where β$\beta$ is the inverse of the turbulent Schmidt number and y   is the distance from the bed). In order to predict the near-bed upward convex profile, an additional parameter α$\alpha$ is needed. This parameter could be related to settling velocity (α$\alpha$ equal to inverse of dimensionless settling velocity) or to convective sediment entrainment process. The profiles are interpreted by a relation between second derivative of the logarithm of concentration and derivative of the product between sediment diffusivity and α$\alpha$.     

Seismic imaging of the middle and upper crust by double-difference tomography: the Haicheng earthquake (Ms 7.3) in Liaoning Province

The Haicheng earthquake (Ms 7.3) occurred in Liaoning Province (39°N–43°N, 120°E–126°E ), China on February 4, 1975. The mortality rate was only 0.02% owing to the first timely and accurate prediction, although the area affected by the earthquake was 9200 km² and covered cities with a population density of 1000 p/km². In this study, the doubledifference (DD) tomography method was used to obtain high-resolution three-dimensional (3D) P- and S-wave velocity (Vp and Vs) structures and Vp/Vs as well as the earthquake locations. Tomography results suggest that velocity structure at shallow depth coincides well with topography and sediment thickness. The earthquake locations form a northwest-striking zone associated with the Jinzhou(JZ) Fault and a northeast-striking zone associated with the Haichenghe-Dayanghe (HD) Fault, and suggest that the JZ Fault consists of three faults and the Ms 7.3 Haicheng earthquake originated at the intersection of the JZ and the Faults. Lowvelocity zones (LVZs) with low Vp/Vs are observed at 15–20 km depth beneath the Haicheng (HC) region. We interpret the LVZs in the middle crust as regions of fluids, suggesting rock dehydration at high temperatures. The LVZs and low Vp/Vs in the upper crust are attributed to groundwater-filled cracks and pores. We believe that large crustal earthquakes in this area are caused by the combination of faulting and fluid movement in the middle crust.     

Joint up/down decomposition and reconstruction using three‐component streamers with or without ghost model: the sampling theory

This paper addresses two artefacts inherent to marine towed‐streamer surveys: 1) ghost reflections and 2) too sparse a sampling in the crossline direction. A ghost reflection is generated when an upcoming reflection bounces off the sea surface back into the sensors and can, in principle, be removed by decomposing the measured wavefield into its up‐ and downgoing constituents. This process requires a dense sampling of the wavefield in both directions along and perpendicular to the streamers. A dense sampling in the latter direction is, however, often impossible due to economical and operational constraints. Recent multi‐component streamers have been designed to record the spatial gradients on top of the pressure, which not only benefits the wavefield decomposition but also facilitates a lower‐than‐Nyquist sampling rate of the pressure. In this paper, wavefield reconstruction and deghosting are posed as a joint inverse problem. We present two approaches to establish a system matrix that embeds both a deghosting and an interpolation operator. The first approach is derived with a ghost model, whereas the second approach is derived without a ghost model. The embodiment of a ghost model leads to an even lower sampling rate but relies on a more restrictive assumption on the sea surface.