随钻声波测井固井质量评价理论与数值模拟研究

大斜度井、水平井的固井质量评价需求,推动了基于随钻测井仪器的固井质量评价技术的发展.本文基于弹性波传播理论,建立柱状分层介质的套管井随钻测井模型,模拟了不同水泥胶结情况的随钻声场响应,特别是其波形特征.结果表明,套管井中随钻与电缆波形特征最大的区别在于随钻波列中存在隔声衰减后的钻铤波.特别是当水泥胶结较好时,钻铤波与套管波幅度大小相近,甚至钻铤波幅度更大,而时间域上两者波形重叠,该特点使得常规电缆声波（声幅/变密度）固井质量评价方法不适用于随钻测井环境.本文考察了不同水泥胶结程度下随钻波列中套管波幅度和衰减与胶结指数的关系,形成一种同时利用幅度和衰减的随钻声波胶结指数计算方法,为随钻声波固井质量评价提供理论指导.

     

Shaking-table tests and numerical simulations on a subway structure in soft soil

Shaking table tests were performed to investigate the damage mechanisms of a subway structure in soft soil while experiencing strong ground motions. The seismic responses of the structure and soil were found to be more sensitive to input motions with richer low-frequency components. The excess pore pressure ratio of soil increased slightly, and the distribution of the excess pore pressure surrounding the structure showed clear spatial effects. The frequency spectrum characteristics of input ground motions clearly influenced the lateral displacement of the structure. In addition, the structure was most severely damaged at the top or the bottom of the interior columns. Finite element analyses were conducted by using the modified Martin–Seed–Davidenkov viscoelastic and the rate-independent plastic-damage constitutive models for soil and concrete, respectively. Satisfactory agreement was observed between the simulation and test results, the difference between these results was discussed in detail. The results provide insight into how the characteristics of strong ground motion might influence and present a simplified analysis method to quantitatively evaluate the damage of subway structures in soft soil.     

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

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

Stability of finite difference numerical simulations of acoustic logging-while-drilling with different perfectly matched layer schemes

In acoustic logging-while-drilling （ALWD） finite difference in time domain （FDTD） simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns （radius -27 mm）. Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole （the difference is 〉30 times）, the stability and efficiency of the perfectly matched layer （PML） scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain （FDTD） in the ALWD simulation, including field-splitting PML （SPML）, multiaxial PML（M- PML）, non-splitting PML （NPML）, and complex frequency-shifted PML （CFS-PML）. The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor （a0） and the scaling factor （β0） depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.     

Measuring water velocity in highly turbulent flows: field tests of an electromagnetic current meter (ECM) and an acoustic Doppler velocimeter (ADV)

Two field tests were completed to compare the performance of an electromagnetic current meter (ECM) with that of an acoustic Doppler velocimeter (ADV) in gravel‐bed rivers. Research was particularly motivated by the need to measure flow properties in highly energetic turbulent flows. Measurements were made at two field sites, one at moderate velocities (up to 70 cm/s) and with moderate turbulence intensities (10–20% of mean flow), and the other in an area of non‐uniform flow that included locations with fast mean velocities (up to 1.75 m/s) and high turbulent intensities (up to 50% of mean flow). Comparison of means, standard deviations, turbulent kinetic energy and Reynolds shear stress confirm the general agreement between the ECMs and ADVs. The general agreement is subject to limitations associated with the sample volume and frequency response of the instruments, and only applies within restricted velocity (up to ≈1.25 m/s) and turbulence intensity ranges (up to ≈0·125 m/s). At higher turbulence intensities, spectral analysis showed anomalous behavior of the ADV signal, especially in the vertical velocity component. Quadrant analysis of the Reynolds stress suggests that these problems occur predominantly in quadrants 1 and 3. Errors in ADV measurements were estimated using four different methods: one that utilized the characteristic noise floor in spectral plots, one based on internal ADV measurements of signal correlation and two techniques that aggregate errors related to various sub‐factors. Estimates were divergent at high flows. Techniques that rely on sub‐factors appeared to underestimate the impact of high turbulence on signal quality. The key conclusion for future field applications is that the older ECM technology provides the more reliable estimates of flow parameters in high turbulence. Copyright © 2007 John Wiley & Sons, Ltd.     

利用阵列声波进行固井质量定量评价的方法及其应用

声波法固井质量测井技术一般采用短源距测量方式测量套管波的幅度、衰减或套后声阻抗来进行固井质量评价,然而在油气田开发生产期间,用户往往提出地层声波信息和固井质量同时评价的要求,此时短源距测量方式已不能满足用户需求.在固井质量好的情况下,阵列声波测井能够得到地层声波信息,但对固井质量评价涉及较少.为满足客户需求,充分利用阵列波形中的信息,降低测井成本,本文基于阵列声波测井资料,开展了固井质量定量评价方法研究,首先数值模拟研究了不同固井质量情况下的阵列声波声场响应特征,然后利用阵列仪器多接收器的特点,引入了一种定量计算套管波衰减的方法,并结合常规固井质量CBL/VDL测井资料,建立了新的固井质量评价标准,实现了阵列声波固井质量的定量评价.通过处理现场测井资料,验证了该方法的正确性和有效性.     

Shaking table tests and numerical analyses of an RC coupled wall structure with replaceable coupling beams

The replaceable coupling beam (RCB) is an innovative structural component developed to increase the seismic resilience of reinforced concrete (RC) shear wall structures. In this study, two 1/5‐scale 5‐story 3‐dimensional RC shear wall structures—one with conventional RC coupling beams and the other with RCBs—were designed, constructed, and tested on a shaking table. The failure pattern, dynamic properties, and structural responses, including the acceleration, displacement, story force, and strain responses, of the 2 structures are compared under earthquake excitations. The test results demonstrate that the seismic performance of the structure with RCBs was improved when RCBs were working compared with the structure with conventional RC coupling beams. In addition, the replaceable devices suffering the severe damage during an earthquake can be conveniently replaced after the earthquake. However, after the sudden failure of RCBs during the severe earthquakes, the inter‐story drift and floor acceleration of the structure with RCBs became larger. The design and manufacture quality of RCBs should be improved to avoid the sudden failure. Then, numerical models for the test structures were established using the commercial software PERFORM‐3D. Numerical simulations of the tests were conducted. The simulation results correspond well with the experimental results, thus verifying the accuracy of the numerical models. The RC shear wall structure installed with RCBs can be applied as a new type of earthquake‐resilient structure in engineering practice.     

Seasonal and decadal subsurface thaw dynamics of an Aufeis feature investigated through numerical simulations

Aufeis (also known as icings) are large sheet-like masses of layered ice that form in river channels in arctic environments in the winter as groundwater discharges to the land surface and subsequently freezes. Aufeis are important sources of water for Arctic river ecosystems, bolstering late summer river discharge and providing habitat for caribou escaping insect harassment. The aim of this research is to use numerical simulations to evaluate a conceptual model of subsurface hydrogeothermal conditions that can lead to the formation of aufeis. We used a conceptual model based on geophysical data from the Kuparuk aufeis field on the North Slope of Alaska to develop a two-dimensional heterogeneous vertical profile model of groundwater flow, heat transport, and freeze/thaw dynamics. Modelling results showed that groundwater can flow to the land surface through subvertical high permeability pathways during winter months when the lower permeability soils near the land surface are frozen. The groundwater discharge can freeze on the surface, contributing to aufeis formation throughout the winter. We performed sensitivity analyses on subsurface properties and surface temperature and found that aufeis formation is most sensitive to the volume of unfrozen water available in the subsurface and the rate at which the subsurface water travels to the land surface. Although a trend of warming air temperatures will lead to a greater volume of unfrozen subsurface water, the aufeis volume can be reduced under warming conditions if the period of time for which air temperatures are below freezing is reduced.     

Shake-table tests and numerical simulation of an innovative isolation system for highway bridges

Damage investigation of small to medium-span highway bridges in Wenchuan earthquake revealed that typical damage of these bridges included: sliding between laminated-rubber bearings and bridge girders, concrete shear keys failure, excessive girder displacements and even span collapse. However, the bearing sliding could actually act as a seismic isolation for piers, and hence, damage to piers for these bridges was minor during the earthquake. Based on this concept, an innovative solation system for highway bridges with laminated-rubber bearings is developed. The system is comprised of typical laminated-rubber bearings and steel dampers. Bearing sliding is allowed during an earthquake to limit the seismic forces transmitting to piers, and steel dampers are applied to restrict the bearing displacements through hysteretic energy dissipation. As a major part of this research, a quarter-scale, two-span bridge model was constructed and tested on the shake tables to evaluate the performance of this isolation system. The bridge model was subjected to a Northridge and an artificial ground motion in transverse direction. Moreover, numerical analyses were conducted to investigate the seismic performance of the bridge model. Besides the test bridge model, a benchmark model with the superstructure fixed to the substructure in transverse direction was also included in the numerical analyses. Both the experimental and the numerical results showed high effectiveness of this proposed isolation system in the bridge model. The system was found to effectively control the pier-girder relative displacements, and simultaneously, protect the piers from severe damage. Numerical analyses also validated that the existing finite element methods are adequate to estimate the seismic response of bridges with this isolation system.     

Identification and mapping of riverbed sediment facies in the Columbia River through integration of field observations and numerical simulations

In the Hanford Reach of the Columbia River, a thin layer of recent alluvium overlies the sedimentary formations that comprise the unconfined groundwater aquifer. Experimental and modelling studies have demonstrated that this alluvial layer exerts significant control on the exchange of groundwater and surface water (hydrologic exchange flux), and is associated with elevated levels of biogeochemical activity. This layer is also observed to be strongly heterogeneous, and quantifying the spatial distribution of properties over the range of scales of interest is challenging. Facies are elements of a sediment classification scheme that groups complex geologic materials into a set of discrete classes according to distinguishing features. Facies classifications have been used as a framework for assigning heterogeneous material properties to grid cells of numerical models of flow and reactive transport in subsurface media. The usefulness of such an approach hinges on being able to relate facies to quantitative properties needed for flow and reactive transport modelling, and on being able to map facies over the domain of interest using readily available information. Although aquifer facies have been used in various modelling contexts, application of this concept to riverbed sediments is relatively new. Here, we describe an approach for categorizing and mapping recent alluvial (riverbed) sediments based on the integration of diverse observations with numerical simulations of river hydrodynamics. The facies have distinct distributions of sediment texture that correspond to variations in hydraulic properties, and therefore provide a useful framework for assigning heterogeneous properties in numerical simulations of hydrologic exchange flows and biogeochemical processes.     

Multifluid flow in bedded porous media: laboratory experiments and numerical simulations

Understanding light nonaqueous-phase liquid (LNAPL) movement in heterogeneous vadose environments is important for effective remediation design. We investigated LNAPL movement near a sloping fine- over coarse-grained textural interface, forming a capillary barrier. LNAPL flow experiments were performed in a glass chamber (50 cm×60 cm×1.0 cm) using two silica sands (12/20 and 30/40 sieve sizes). Variable water saturations near the textural interface were generated by applying water uniformly to the sand surface at various flow rates. A model LNAPL (Soltrol® 220) was subsequently released at two locations at the sand surface. Visible light transmission was used to quantitatively determine water saturations prior to LNAPL release and to observe LNAPL flow paths. Numerical simulations were performed using the Subsurface Transport Over Multiple Phases (STOMP) simulator, employing two nonhysteretic relative permeability–saturation–pressure (k–S–P) models. LNAPL movement strongly depended on the water saturation in the fine-grained sand layer above the textural interface. In general, reasonable agreement was found between observed and predicted water saturations near the textural interface and LNAPL flow paths. Discrepancies between predictions based on the van Genuchten/Mualem (VGM) and Brooks–Corey/Burdine (BCB) k–S–P models existed in the migration speed of the simulated LNAPL plume and the LNAPL flow patterns at high water saturation above the textural interface. In both instances, predictions based on the BCB model agreed better with experimental observations than predictions based on the VGM model. The results confirm the critical role water saturation plays in determining LNAPL movement in heterogeneous vadose zone environments and that accurate prediction of LNAPL flow paths depends on the careful selection of an appropriate k–S–P model.     

Remediation of DNAPL source zones with granular iron: laboratory and field tests

Degradation of dissolved chlorinated solvents using granular iron is an established in situ technology. This paper reports on investigations into mixing iron and bentonite with contaminated soil for in situ containment and degradation of dense nonaqueous phase liquid source zones. In the laboratory, hypovials containing soil, water, bentonite, iron, and free-phase trichloroethene (TCE) were assembled. Periodic measurement of TCE, chloride, and degradation products showed progressive degradation of TCE to nondetectable levels. Subsequently, a demonstration was conducted at Canadian Forces Base Borden near Alliston, Ontario, Canada, where, in 1991, a portion of the surficial aquifer was isolated and free-phase tetrachloroethene (PCE) was introduced. Using a drill rig equipped with large-diameter mixing blades, three mixed zones were prepared containing 0%, 5%, and 10% granular iron by volume. The bentonite was added to serve as a lubricant to facilitate injection of the iron and to isolate the contaminated zone. Analysis of core samples showed reasonably uniform distributions of iron through the mixed zones. Monitoring over a 13-month period following installation showed, relative to the control, a decline in PCE concentrations to virtually nondetectable values. Reaction rates in the laboratory tests were similar to those reported in the literature, while the rate in the field test was substantially lower. The lower rate may be a consequence of mass transfer limitations under the static conditions of the field test. Results indicate that mixing iron and bentonite into source zones may be an effective means of source-zone remediation, with the particular advantage of being relatively immune to effects of geologic heterogeneity.     

Solar flare model: Comparison of the results of numerical simulations and observations

The electrodynamic flare model is based on numerical 3D simulations with the real magnetic field of an active region. An energy of ∼10³² erg necessary for a solar flare is shown to accumulate in the magnetic field of a coronal current sheet. The thermal X-ray source in the corona results from plasma heating in the current sheet upon reconnection. The hard X-ray sources are located on the solar surface at the loop foot-points. They are produced by the precipitation of electron beams accelerated in field-aligned currents. Solar cosmic rays appear upon acceleration in the electric field along a singular magnetic X-type line. The generation mechanism of the delayed cosmic-ray component is also discussed.     

HDR devices for the seismic protection of frame structures: Experimental results and numerical simulations

As part of a national research programme an experimental campaign was carried out on a real scale mock‐up consisting of a steel–concrete composite frame equipped with dissipative bracings, based on high damping rubber (HDR) devices. Free vibration tests, followed by force‐controlled and displacement‐controlled cyclic tests were performed. The experimental tests were aimed at studying the dynamic response of the coupled system in order to demonstrate the effectiveness of HDR devices in increasing the stiffness and dissipation capacity of the frame and investigating the ability of the constitutive HDR model proposed by the authors to predict the dynamic response of the coupled system. Copyright © 2009 John Wiley & Sons, Ltd.     

The emplacement of pahoehoe toes: field observations and comparison to laboratory simulations

We observed active pahoehoe lobes erupted on Kilauea during May-June 1996, and found a range of emplacement styles associated with variations in local effusion rate, flow velocity, and strain rate. These emplacement styles were documented and quantified for comparison with earlier laboratory experiments.At the lowest effusion rates, velocities, and strain rates, smooth-surfaced lobes were emplaced via swelling, where new crust formed along an incandescent lip at the front of the lobe and the rest of the lobe was covered with a dark crust. At higher effusion rates, strain rates and velocities, lobes were emplaced through tearing or cracking. Tearing was characterized by ripping of the ductile crust near the initial breakout point, and most of the lobe surface was incandescent during its emplacement. This mechanism was observed to generate both smooth-surfaced lobes, and, when the lava encountered an obstacle, folded lobes. Cracking lobes were similar to those emplaced via tearing, but involved breaking of a thicker, brittle crust at the initial breakout of the lobe and therefore required somewhat higher flow rates than did tearing. Cracking lobes typically formed ropy folds in the center of the lobe, and smooth margins. At the highest effusion rates, strain rates, and flow velocities, the lava formed open channels with distinct levees.The final lobe morphologies were compared to results from laboratory simulations, which were designed to infer effusion rate from final flow morphology, to quantitatively test the laboratory results on the scale of individual natural pahoehoe lobes. There is general agreement between results from laboratory simulations and natural lavas on the scale of individual pahoehoe lobes, but there are disparities between laboratory flows and lava flows on the scale of an entire pahoehoe lava flow field.Editorial responsibility: A. Woods     

Pulled array continuous electrical sounding with an additional inductive source: an experimental design study

Electrical and electromagnetic profiling methods are used extensively in environmental geophysical investigations for many different purposes. The pulled array continuous electrical sounding (PACES) method, where a tail of electrodes is towed behind a small vehicle while continuously and simultaneously measuring several electrode configurations, has been used extensively for mapping the vulnerability of aquifers in Denmark. Measurements are taken every 1 m, and 10–15 km of profile can be achieved in one day.
This paper presents a theoretical study of the resolution capabilities of PACES measurements as they are now performed, and an experimental design study for including an inductive source in the measuring equipment. The joint interpretation of the galvanic data set of ordinary PACES measurements with inductive data from a horizontal magnetic dipole source will enhance the resolution capabilities of the data set significantly. The study is carried out as an analysis of the uncertainty of the model parameters of one-dimensional three-layer models using the estimation error variances of the inversion problem.
The results indicate that the addition of only two frequency data from a magnetic dipole source will substantially improve the resolution of the subsurface resistivity structure. The improvement is model dependent, but reduction in the relative error of model parameters by an order of magnitude is observed.     

Island-enhanced cooling mechanism in typhoon events revealed by field observations and numerical simulations for a coral reef area,Sekisei Lagoon,Japan

While widely known for their destructive power, typhoon events can also bring benefit to coral reef ecosystems through typhoon-induced cooling which can mitigate against thermally stressful conditions causing coral bleaching. Sensor deployments in Sekisei Lagoon, Japan’s largest coral reef area, during the summer months of 2013, 2014, and 2015 were able to capture local hydrodynamic features of numerous typhoon passages. In particular, typhoons 2015-13 and 2015-15 featured steep drops in near-bottom temperature of 5 °C or more in the north and south sides of Sekisei Lagoon, respectively, indicating local cooling patterns which appeared to depend on the track and intensity of the passing typhoon. This was further investigated using Regional Ocean Modeling System (ROMS) numerical simulations conducted for the summer of 2015. The modeling results showed a cooling trend to the north of the Yaeyama Islands during the passage of typhoon 2015-13, and a cooling trend that moved clockwise from north to south of the islands during the passage of typhoon 2015-15. These local cooling events may have been initiated by the Yaeyama Islands acting as an obstacle to a strong typhoon-generated flow which was modulated and led to prominent cooling of waters on the leeward sides. These lower temperature waters from offshore may then be transported to the shallower inner parts of the lagoon area, which may partly be due to density-driven currents generated by the offshore-inner area temperature difference.     

Magnetohydrodynamics (MHD) numerical simulations on the interaction of the solar wind with the magnetosphere: A review

The magnetosphere is the outermost layer of the geospace, and the interaction of the solar wind with the magnetosphere is the key element of the space weather cause-and-effect chain process from the Sun to Earth, which is one of the most challenging scientific problems in the geospace weather study. The nonlinearity, multiple component, and time-dependent nature of the geospace make it very difficult to describe the physical process in geospace using traditional analytic analysis approach. Numerical simulations, a new research tool developed in recent decades, have a deep impact on the theory and application of the geospace. MHD simulations started at the end of the 1970s, and the initial study was limited to two-dimensional (2D) cases. Due to the intrinsic three-dimensional (3D) characteristics of the geospace, 3D MHD simulations emerged in the 1980s, in an attempt to model the large-scale structures and fundamental physical processes in the magnetosphere. They started to combine with the space exploration missions in the 1990s and make comparisons with observations. Physics-based space weather forecast models started to be developed in the 21st century. Currently only a few space-power countries such as USA and Japan have developed 3D magnetospheric MHD models. With the rapid advance of space science in China, we have developed a new global MHD model, namely PPMLR-MHD, which has high order spatial accuracy and low numerical dissipation. In this review, we will briefly introduce the global 3D MHD modeling, especially the PPMLR-MHD code, and summarize our recent work based on the PPMLR-MHD model, with an emphasis on the interaction of interplanetary shocks with the magnetosphere, large-scale current systems, reconnection voltage and transpolar potential drop, and Kelvin-Helmholtz (K-H) instability at the magnetopause.