非对称声源多极子随钻声波测井实验室测量研究

现有的随钻声波测井通常采用三种测量方式:单极子、偶极子及四极子测量模式.分别采用单极子源、偶极子源和四极子源激发单极子、偶极子和四极子信号对地层速度信息进行测量.由于随钻测井环境比较复杂,钻杆在井孔中占据大部分空间,钻杆速度高于地层速度,导致接收到的信号中仪器波能量占主导,无法有效地提取地层的纵横波速度信息.本文探讨一种采用非对称偏心的点声源,一次激发,利用方位分布的四个接收器接收,经过合成后同时分别得到单极子、偶极子和四极子模式波信号的测量方法,为建立同时利用多极模式波信号测量地层纵横波速度的随钻声波测井实用技术提供基础.本文利用1 ∶12的缩比模型在实验室对该方法进行了考察,实验结果证明,可以通过非对称点声源得到有效的地层单极子、偶极子和四极子模式波信号.通过分析和对比实测数据和理论模拟的频散曲线,本文进一步提出利用三种模式波的频散曲线联合反演确定地层横波速度的方法.     

3D resistivity inversion using 2D measurements of the electric field

Field and 'noisy' synthetic measurements of electric-field components have been inverted into 3D resistivities by smoothness-constrained inversion. Values of electrical field can incorporate changes in polarity of the measured potential differences seen when 2D electrode arrays are used with heterogeneous 'geology', without utilizing negative apparent resistivities or singular geometrical factors. Using both the X - and Y -components of the electric field as measurements resulted in faster convergence of the smoothness-constrained inversion compared with using one component alone. Geological structure and resistivity were reconstructed as well as, or better than, comparable published examples based on traditional measurement types. A 2D electrode grid (20 × 10), incorporating 12 current-source electrodes, was used for both the practical and numerical experiments; this resulted in 366 measurements being made for each current-electrode configuration. Consequently, when using this array for practical field surveys, 366 measurements could be acquired simultaneously, making the upper limit on the speed of acquisition an order of magnitude faster than a comparable conventional pole–dipole survey. Other practical advantages accrue from the closely spaced potential dipoles being insensitive to common-mode noise (e.g. telluric) and only 7% of the electrodes (i.e. those used as current sources) being susceptible to recently reported electrode charge-up effects.     

Improved understanding of spatio‐temporal controls on regional scale groundwater flooding using hydrograph analysis and impulse response functions

Controls on the spatio‐temporal extent of groundwater flooding are poorly understood, despite the long duration of groundwater flood events and distinct social and economic impacts. We developed a novel approach using statistical analysis of groundwater level hydrographs and impulse response functions (IRFs) and applied it to the 2013/2014 Chalk groundwater flooding in the English Lowlands. We proposed a standardized index of groundwater flooding which we calculated for monthly groundwater levels for 26 boreholes in the Chalk. We grouped these standardized series using k‐means cluster analysis and cross‐correlated the cluster centroids with the Standardized Precipitation Index accumulated over time intervals between 1 and 60 months. This analysis reveals 2 spatially coherent groups of standardized hydrographs that responded to precipitation over different timescales. We estimated IRF models of the groundwater level response to effective precipitation for 3 boreholes in each group. The IRF models corroborate the Standardized Precipitation Index analysis showing different response functions between the groups. We applied identical effective precipitation inputs to each of the IRF models and observed differences between the hydrographs from each group. It is suggested this is due to the hydrogeological properties of the Chalk and of overlying relatively low permeability superficial deposits (recent unconsolidated sediments overlying the bedrock, such as clays and tills), which are extensive over 1 of the groups. The overarching controls on groundwater flood response are concluded to be a complex combination of antecedent conditions, rainfall, and catchment hydrogeological properties. These controls should be taken into consideration when anticipating and managing future groundwater flood events. The approach presented is generic and parsimonious and can be easily applied where sufficient groundwater level and rainfall data are available.     

3D acoustic waveform inversion in the Laplace domain using an iterative solver

In this paper we propose a 3D acoustic full waveform inversion algorithm in the Laplace domain. The partial differential equation for the 3D acoustic wave equation in the Laplace domain is reformulated as a linear system of algebraic equations using the finite element method and the resulting linear system is solved by a preconditioned conjugate gradient method. The numerical solutions obtained by our modelling algorithm are verified through a comparison with the corresponding analytical solutions and the appropriate dispersion analysis. In the Laplace‐domain waveform inversion, the logarithm of the Laplace transformed wavefields mainly contains long‐wavelength information about the underlying velocity model. As a result, the algorithm smoothes a small‐scale structure but roughly identifies large‐scale features within a certain depth determined by the range of offsets and Laplace damping constants employed. Our algorithm thus provides a useful complementary process to time‐ or frequency‐domain waveform inversion, which cannot recover a large‐scale structure when low‐frequency signals are weak or absent. The algorithm is demonstrated on a synthetic example: the SEG/EAGE 3D salt‐dome model. The numerical test is limited to a Laplace‐domain synthetic data set for the inversion. In order to verify the usefulness of the inverted velocity model, we perform the 3D reverse time migration. The migration results show that our inversion results can be used as an initial model for the subsequent high‐resolution waveform inversion. Further studies are needed to perform the inversion using time‐domain synthetic data with noise or real data, thereby investigating robustness to noise.     

一种基于阵列声波测井的胶结指数计算方法在饱和度评价中的应用

利用饱和度公式评价储层含油性时,胶结指数一般在某一层位取某一固定值.在复杂岩性储层中,储层非均质性强,孔隙结构复杂,连通性差,导致胶结指数变化较大,利用常数胶结指数计算的含油饱和度精度降低,给储层测井评价带来了极大的困难.基于弹性介质理论,利用阵列声波数据得到地层骨架的剪切模量和体积模量,并计算出骨架的纵波速度.分析发现,地层和骨架的纵波速度比与地层孔隙度和胶结指数具有一定的函数关系,通过地层和骨架的纵波速度比和孔隙度可以计算出连续的胶结指数,进而得到地层饱和度.实验结果表明,本文方法计算得到的地层饱和度与岩心分析饱和度吻合较高,方法一定程度上克服了传统的饱和度模型在复杂岩性储层中的不适用性,为复杂岩性储层含油性准确评价提供了依据.     

4D tropospheric tomography using GPS slant wet delays

Tomographic techniques are successfully applied to obtain 4D images of the tropospheric refractivity in a local dense network of global positioning system (GPS) receivers. We show here how GPS data are processed to obtain the tropospheric slant wet delays and discuss the validity of the processing. These slant wet delays are the observables in the tomographic processing. We then discuss the inverse problem in 4D tropospheric tomography making extensive use of simulations to test the system and define the resolution and the impact of noise. Finally, we use data from the Kilauea network in Hawaii for February 1, 1997, and a local 4 × 4 × 40 voxel grid on a region of 400 km² and 15 km in height to produce the corresponding 4D wet refractivity fields, which are then validated using forecast analysis from the European Center for Medium Range Weather Forecast (ECMWF). We conclude that tomographic techniques can be used to monitor the troposphere in time and space.     

A new method for quantification of liquid saturation in 2D translucent porous media systems using light transmission

Five physically based models for predicting liquid saturation from light transmission in 2D laboratory systems containing translucent porous media were developed and tested (Models A–E). The models were based upon various simplifying assumptions concerning pore geometry, wettability, and drainage. Models A–D assumed uniform-sized pores, and liquid saturation was an explicit function of light transmission. Model E considered a distribution of pore sizes whose drainage characteristics were inferred from the Laplace equation. Mass balances were calculated using data from drainage and infiltration experiments, in four textures of silica sand with water as the fluid. Model E performed the best overall, with systematic errors of less than 2.3% saturation. Model E represents a robust easily applied new method for the determination of liquid saturation by light transmission. The other four models are presented, and compared, for reasons of historical interest and to investigate the impact of the various simplifying assumptions.     

COSMIC数据验证FY-3 MWTS通道4的观测结果

利用无线电掩星探测在平流层具有较高精度和稳定性的优势,本文采用COSMIC资料验证我国风云3号（FY-3）系列卫星上微波温度探测仪（MWTS）通道4表征的平流层低层测量结果.两年的FY-3A/3B上MWTS通道4观测亮温与掩星数据模拟的亮温相比较分析表明：两颗卫星平台上MWTS通道4观测亮温都偏高,尤其在热带地区和极区夏季;FY3A/3B上MWTS通道4亮温偏差的月均值变化趋势较为一致,两个年份的变化趋势也较为一致.MWTS通道4表征的平流层低层观测亮温在不同纬度带呈现不同偏差分布特征：热带地区明显系统偏高2~4 K,中纬度地区偏高1 K;而高纬度地区则出现随季节显著变化的偏差,尤其在南极地区季节差异达到5 K,这种依赖于环境温度出现的较大亮温偏差不完全是频率偏移引起.热带地区观测与模拟亮温差异显著,并且热带地区样本对于全球总体亮温偏差影响程度达到20%,这表明COSMIC数据在热带地区的验证结果需谨慎使用.

     

Improved design of sliding mode control for civil structures with saturation problem

A systematic and improved design procedure for sliding mode control (SMC) of seismically excited civil structures with saturation problem is provided in this paper. In order to restrict the control force to a certain level, a procedure for determining the upper limits of the control forces for single or multiple control units is proposed based on the design response spectrum of external loads. Further, an efficient procedure using the LQR method for determining sliding surfaces appropriate for different controller types is provided through the parametric evaluation of the dynamic characteristics of sliding surfaces in terms of SMC controller performance. Finally, a systematic design procedure for SMC required to achieve a given performance level is provided and its effectiveness is verified by applying it to multi‐degree‐of‐freedom (MDOF) systems. Copyright © 2004 John Wiley & Sons, Ltd.     

Extreme flood‐driven fluvial bank erosion and sediment loads: direct process measurements using integrated Mobile Laser Scanning (MLS) and hydro‐acoustic techniques

This methods paper details the first attempt at monitoring bank erosion, flow and suspended sediment at a site during flooding on the Mekong River induced by the passage of tropical cyclones. We deployed integrated mobile laser scanning (MLS) and multibeam echo sounding (MBES), alongside acoustic Doppler current profiling (aDcp), to directly measure changes in river bank and bed at high (~0.05 m) spatial resolution, in conjunction with measurements of flow and suspended sediment dynamics. We outline the methodological steps used to collect and process this complex point cloud data, and detail the procedures used to process and calibrate the aDcp flow and sediment flux data. A comparison with conventional remote sensing methods of estimating bank erosion, using aerial images and Landsat imagery, reveals that traditional techniques are error prone at the high temporal resolutions required to quantify the patterns and volumes of bank erosion induced by the passage of individual flood events. Our analysis reveals the importance of cyclone‐driven flood events in causing high rates of erosion and suspended sediment transport, with a c. twofold increase in bank erosion volumes and a fourfold increase in suspended sediment volumes in the cyclone‐affected wet season. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimating uncertainty of streamflow measurements with moving-boat acoustic Doppler current profilers

A theoretical and a semi-empirical model are presented for estimating the uncertainty of streamflow measurements made with an acoustic Doppler current profiler (ADCP) mounted on a moving platform. Both models are based on the statistical analysis of ADCP ensemble discharge time series collected during a transect; therefore, they account for all random error sources associated with the measurement at a site. The theoretical model is developed based on the law of propagation of variance; it explores the theoretical relationship between the variables involved in the problem. The semi-empirical model is developed based on the theory of dimensional analysis; it explores the empirical relationship between the variables. The semi-empirical model is calibrated using 205 transect datasets and verified with an additional 382 transect datasets. It provides a useful tool for the uncertainty analysis and uncertainty-based measurement quality control of moving-boat ADCP streamflow measurements.     

Gas saturation prediction and effect of low frequencies on acoustic impedance images at Foinaven Field

Low frequencies are necessary in seismic data for proper acoustic impedance imaging and for petrophysical interpretation. Without lower frequencies, images can be distorted leading to incorrect reservoir interpretation and petrophysical predictions. As part of the Foinaven Active Reservoir Management (FARM) project, a Towed Streamer survey and an Ocean Bottom Hydrophone (OBH) survey were shot in both 1995 and 1998. The OBH surveys contain lower frequencies than the streamer surveys, providing a unique opportunity to study the effects that low frequencies have on both the acoustic impedance image along with petrophysical time‐lapse predictions. Artefacts that could easily have been interpreted as high‐resolution features in the streamer data impedance volumes can be distinguished by comparison with the impedance volumes created from the OBH surveys containing lower frequencies. In order to obtain results from the impedance volumes, impedance must be related to saturation. The mixing of exsolved gas, oil and water phases involves using the Reuss (uniform) or Voigt (patchy approximation) mixing laws. The Voigt average is easily misused by assuming that the end‐points correspond to 0% and 100% gas saturation. This implies that the patches are either 0% gas saturation or 100% gas saturation, which is never the case. Here, the distribution of gas as it comes out of solution is assumed to be uniform until the gas saturation reaches a sufficiently high value (critical gas saturation) to allow gas to flow. Therefore, at low gas saturations the distribution is uniform, but at saturations above critical, it is patchy, with patches that range from critical gas saturation to the highest gas saturation possible (1 minus residual oil and irreducible water saturation).     

Determination of porosity and saturation using seismic waveform inversion

Characterization of a reservoir model requires determination of its petrophysical parameters, such as porosity and saturation. We propose a new method to determine these parameters directly from seismic data. The method consists of the computation and inversion of seismic waveforms. A high frequency method is presented to model wave propagation through an attenuative and dispersive poroelastic medium. The high frequency approximation makes it possible to efficiently compute sensitivity functions. This enables the inversion of seismic waveforms for porosity and saturation. The waveform inversion algorithm is applied to two laboratory crosswell datasets of a water saturated sand. The starting models were obtained using travel time tomography. The first dataset is inverted for porosity. The misfit reduction for this dataset is approximately 50%. The second dataset was obtained after injection of a nonaqueous-phase liquid (NAPL), possibly with some air, which made the medium more heterogeneous. This dataset was inverted for NAPL and air saturation using the porosity model obtained from the first inversion. The misfit reduction of the second experiment was 70%. Regions of high NAPL and high air saturation were found at the same location. These areas correlate well with the position of one of the injection points as well as regions of higher NAPL concentrations found after excavation of the sand. It is therefore possible to directly invert waveforms for pore fluid saturation by taking into account the attenuation and dispersion caused by the poroelasticity.     

Controlled‐source electromagnetic monitoring of reservoir oil saturation using a novel borehole‐to‐surface configuration

To advance and optimize secondary and tertiary oil recovery techniques, it is essential to know the areal propagation and distribution of the injected fluids in the subsurface. We investigate the applicability of controlled‐source electromagnetic methods to monitor fluid movements in a German oilfield (Bockstedt, onshore Northwest Germany) as injected brines (highly saline formation water) have much lower electrical resistivity than the oil within the reservoir. The main focus of this study is on controlled‐source electromagnetic simulations to test the sensitivity of various source–receiver configurations. The background model for the simulations is based on two‐dimensional inversion of magnetotelluric data gathered across the oil field and calibrated with resistivity logs. Three‐dimensional modelling results suggest that controlled‐source electromagnetic methods are sensitive to resistivity changes at reservoir depths, but the effect is difficult to resolve with surface measurements only. Resolution increases significantly if sensors or transmitters can be placed in observation wells closer to the reservoir. In particular, observation of the vertical electric field component in shallow boreholes and/or use of source configurations consisting of combinations of vertical and horizontal dipoles are promising. Preliminary results from a borehole‐to‐surface controlled‐source electromagnetic field survey carried out in spring 2014 are in good agreement with the modelling studies.     

Virtual signals combination – phase analysis for 2D and 3D data representation (acoustic medium)

We formulate the Kirchhoff‐Helmholtz representation theory for the combination of seismic interferometry signals synthesized by cross‐correlation and by cross‐convolution in acoustic media. The approach estimates the phase of the virtual reflections from the boundary encompassing a volume of interest and subtracts these virtual reflections from the total seismic‐interferometry wavefield. The reliability of the combination result, relevant for seismic exploration, depends on the stationary‐phase and local completeness in partial coverage regions. The analysis shows the differences in the phase of the corresponding seismic interferometry (by cross‐correlation) and virtual reflector (by cross‐convolution) signals obtained by 2D and 3D formulations, with synthetic examples performed to remove water layer multiples in ocean bottom seismic (OBS) acoustic data.     

Three‐dimensional modeling of coastal boulders using multi‐view image measurements

Coarse clastic sediments (boulders) on coastlines have seen a groundswell in geomorphic research interest over recent years, associated in part with the potential of boulder evidence for interpreting characteristics of high‐energy wave processes. Yet, the fundamental property of boulder volume is normally difficult to measure accurately owing to complex clast morphology and irregular surface texture. To tackle this problem, this paper concentrates on creating precise, measurable and textured three‐dimensional (3D) models of coastal boulders without physical contact with the object, based on multi‐view image measurement techniques. This method has several advantages over traditional measurements that are inaccurate or alternative solutions using costly techniques such as terrestrial laser scanning. Our methods propose the use of low‐cost equipment (digital cameras) that can be used in various coastal environments to easily acquire numerous images of the object of interest. Initial results can be rapidly assessed in the field for immediate quality control. Resulting 3D models, built from overlapping multi‐view digital photographs, allow the reconstruction of realistic‐looking and textured boulder surfaces. A particular interest in this task is the family of algorithms known as structure from motion (SFM). The work presents analysis of SFM techniques by examining 3D models of boulders observed at a coastal field site on Lu Dao Island in south‐eastern Taiwan. Copyright © 2013 John Wiley & Sons, Ltd.     

2.5D finite-difference solution of the acoustic wave equation

The finite‐difference method applied to the full 3D wave equation is a rather time‐consuming process. However, in the 2.5D case, we can take advantage of the medium symmetry. By taking the Fourier transform with respect to the out‐of‐plane direction (the symmetry axis), the 3D problem can be reduced to a repeated 2D problem. The third dimension is taken into account by a sum over the corresponding wave‐vector component. A criterion for where to end this theoretically infinite sum derives from the stability conditions of the finite‐difference schemes employed. In this way, the computation time of the finite‐difference calculations can be considerably reduced. The quality of the modelling results obtained with this 2.5D finite‐difference scheme is comparable to that obtained using a standard 3D finite‐difference scheme.     

Monitoring CO2 saturation using time‐lapse amplitude versus offset analysis of 3D seismic data from the Ketzin CO2 storage pilot site,Germany

The injection of CO₂ at the Ketzin pilot site commenced in June 2008 and was terminated in August 2013 after 67 kT had been injected into a saline formation at a depth of 630–650 m. As part of the site monitoring program, four 3D surface seismic surveys have been acquired to date, one baseline and three repeats, of which two were conducted during the injection period, and one during the post‐injection phase. The surveys have provided the most comprehensive images of the spreading CO₂ plume within the reservoir layer. Both petrophysical experiments on core samples from the Ketzin reservoir and spectral decomposition of the 3D time‐lapse seismic data show that the reservoir pore pressure change due to CO₂ injection has a rather minor impact on the seismic amplitudes. Therefore, the observed amplitude anomaly is interpreted to be mainly due to CO₂ saturation. In this study, amplitude versus offset analysis has been applied to investigate the amplitude versus offset response from the top of the sandstone reservoir during the injection and post‐injection phases, and utilize it to obtain a more quantitative assessment of the CO₂ gaseous saturation changes. Based on the amplitude versus offset modelling, a prominent decrease in the intercept values imaged at the top of the reservoir around the injection well is indeed associated solely with the CO₂ saturation increase. Any change in the gradient values, which would, in case it was positive, be the only signature induced by the reservoir pressure variations, has not been observed. The amplitude versus offset intercept change is, therefore, entirely ascribed to CO₂ saturation and used for its quantitative assessment. The estimated CO₂ saturation values around the injection area in the range of 40%–60% are similar to those obtained earlier from pulsed neutron‐gamma logging. The highest values of 80% are found in the second seismic repeat in close vicinity to the injection and observation wells.     

Study on the simulation of acoustic logging measurements in horizontal and deviated wells

The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.