Data-driven retrieval of primary plane-wave responses

Seismic images provided by reverse time migration can be contaminated by artefacts associated with the migration of multiples. Multiples can corrupt seismic images, producing both false positives, that is by focusing energy at unphysical interfaces, and false negatives, that is by destructively interfering with primaries. Multiple prediction/primary synthesis methods are usually designed to operate on point source gathers and can therefore be computationally demanding when large problems are considered. A computationally attractive scheme that operates on plane-wave datasets is derived by adapting a data-driven point source gathers method, based on convolutions and cross-correlations of the reflection response with itself, to include plane-wave concepts. As a result, the presented algorithm allows fully data-driven synthesis of primary reflections associated with plane-wave source responses. Once primary plane-wave responses are estimated, they are used for multiple-free imaging via plane-wave reverse time migration. Numerical tests of increasing complexity demonstrate the potential of the proposed algorithm to produce multiple-free images from only a small number of plane-wave datasets.     

Comparing detection and location performance of perpendicular and parallel broadside GPR antenna orientations

GPR (Ground Penetrating Radar) results are shown for perpendicular broadside and parallel broadside antenna orientations. Performance in detection and localization of concrete tubes and steel tanks is compared as a function of acquisition configuration. The comparison is done using 100 MHz and 200 MHz center frequency antennas. All tubes and tanks are buried at the geophysical test site of IAG/USP in São Paulo city, Brazil. The results show that the long steel pipe with a 38-mm diameter was well detected with the perpendicular broadside configuration. The concrete tubes were better detected with the parallel broadside configuration, clearly showing hyperbolic diffraction events from all targets up to 2-m depth. Steel tanks were detected with the two configurations. However, the parallel broadside configuration was generated to a much lesser extent an apparent hyperbolic reflection corresponding to constructive interference of diffraction hyperbolas of adjacent targets placed at the same depth. Vertical concrete tubes and steel tanks were better contained with parallel broadside antennas, where the apexes of the diffraction hyperbolas better corresponded to the horizontal location of the buried target disposition. The two configurations provide details about buried targets emphasizing how GPR multi-component configurations have the potential to improve the subsurface image quality as well as to discriminate different buried targets. It is judged that they hold some applicability in geotechnical and geoscientific studies.     

利用阵列感应测井进行储层渗透率评价

钻井过程中储层受到泥浆侵入影响的程度与储层岩性有着密切关系,其中储层渗透率对侵入深度有着较大影响,因此若可以获知泥浆侵入深度,则有望对储层渗透率进行评估.本文首先建立含泥饼增长的泥浆侵入数值模型,然后建立阵列感应测井数值模型,两者的联合正演模拟显示泥浆侵入对地层的影响可以反映在阵列感应测井响应上,利用阻尼最小二乘法对阵列感应测井响应进行反演可以得到侵入深度.对侵入深度和储层渗透率的关系进行分析发现:在渗透率为1～100mD(1mD=0.987×10~(-3)μm～2)数量级的储层中,渗透率的变化可以在侵入深度上得到反映.以储层和井数据进行二维数值模拟发现:利用阵列感应测井响应反演出来的侵入深度曲线反映了渗透率在地层上的变化趋势,采用解释图版的方法可以对储层各层段的渗透率进行粗略估算.     

Unified multi‐depth‐level field decomposition

Wavefield decomposition forms an important ingredient of various geophysical methods. An example of wavefield decomposition is the decomposition into upgoing and downgoing wavefields and simultaneous decomposition into different wave/field types. The multi‐component field decomposition scheme makes use of the recordings of different field quantities (such as particle velocity and pressure). In practice, different recordings can be obscured by different sensor characteristics, requiring calibration with an unknown calibration factor. Not all field quantities required for multi‐component field decomposition might be available, or they can suffer from different noise levels. The multi‐depth‐level decomposition approach makes use of field quantities recorded at multiple depth levels, e.g., two horizontal boreholes closely separated from each other, a combination of a single receiver array combined with free‐surface boundary conditions, or acquisition geometries with a high‐density of vertical boreholes. We theoretically describe the multi‐depth‐level decomposition approach in a unified form, showing that it can be applied to different kinds of fields in dissipative, inhomogeneous, anisotropic media, e.g., acoustic, electromagnetic, elastodynamic, poroelastic, and seismoelectric fields. We express the one‐way fields at one depth level in terms of the observed fields at multiple depth levels, using extrapolation operators that are dependent on the medium parameters between the two depth levels. Lateral invariance at the depth level of decomposition allows us to carry out the multi‐depth‐level decomposition in the horizontal wavenumber–frequency domain. We illustrate the multi‐depth‐level decomposition scheme using two synthetic elastodynamic examples. The first example uses particle velocity recordings at two depth levels, whereas the second example combines recordings at one depth level with the Dirichlet free‐surface boundary condition of zero traction. Comparison with multi‐component decomposed fields shows a perfect match in both amplitude and phase for both cases. The multi‐depth‐level decomposition scheme is fully customizable to the desired acquisition geometry. The decomposition problem is in principle an inverse problem. Notches may occur at certain frequencies, causing the multi‐depth‐level composition matrix to become uninvertible, requiring additional notch filters. We can add multi‐depth‐level free‐surface boundary conditions as extra equations to the multi‐component composition matrix, thereby overdetermining this inverse problem. The combined multi‐component–multi‐depth‐level decomposition on a land data set clearly shows improvements in the decomposition results, compared with the performance of the multi‐component decomposition scheme.     

Developments in triaxial testing technique

Summary A simple triaxial cell is described together with details of the technique employed for triaxial testing. Triaxial strength results for 254 specimens of eight rock types are tabulated.

---

Zusammenfassung Ein einfaches Triaxial-Prüfgerät wird, zusammen mit Einzelheiten des experimentellen Verfahrens, beschrieben. Ergebnisse dreiachsiger Festigkeitsprüfung sind für 254 Proben aus acht verschiedenen Gesteinsarten in Tabellen wiedergegeben.

---

Résumé Un appareil triaxial simple est décrit avec des détails sur la technique expérimentale. Des résultats de résistance triaxiale pour 254 échantillons pris parmi huit types de roche sont donnés.

---

---

With 2 Figures     

Effects of humic substances on Fe(II) sorption onto aluminum oxide and clay

We studied the effects of humic substances (HS) on the sorption of Fe(II) onto Al-oxide and clay sorbents at pH 7.5 with a combination of batch kinetic experiments and synchrotron Fe K-edge EXAFS analyses. Fe(II) sorption was monitored over the course of 4 months in anoxic clay and Al-oxide suspensions amended with variable HS types (humic acid, HA; or fulvic acid, FA) and levels (0, 1, and 4 wt%), and with differing Fe(II) and HS addition sequences (co-sorption and pre-coated experiments, where Fe(II) sorbate was added alongside and after HS addition, respectively). In the Al-oxide suspensions, the presence of HS slowed down the kinetics of Fe(II) sorption, but had limited, if any, effect on the equilibrium aqueous Fe(II) concentrations. EXAFS analyses revealed precipitation of Fe(II)–Al(III)-layered double hydroxide (LDH) phases as the main mode of Fe(II) sorption in both the HA-containing and HA-free systems. These results demonstrate that HS slow down Fe(II) precipitation in the Al-oxide suspensions, but do not affect the composition or stability of the secondary Fe(II)–Al(III)-LDH phases formed. Interference of HS with the precipitation of Fe(II)–Al(III)-LDH was attributed to the formation organo-Al complexes HS limiting the availability of Al for incorporation into secondary layered Fe(II)-hydroxides. In the clay systems, the presence of HA caused a change in the main Fe(II) sorption product from Fe(II)–Al(III)-LDH to a Fe(II)-phyllosilicate containing little structural Al. This was attributed to complexation of Al by HA, in combination with the presence of dissolved Si in the clay suspension enabling phyllosilicate precipitation. The change in Fe(II) precipitation mechanism did not affect the rate of Fe(II) sorption at the lower HA level, suggesting that the inhibition of Fe(II)–Al(III)-LDH formation in this system was countered by enhanced Fe(II)-phyllosilicate precipitation. Reduced rates of Fe(II) sorption at the higher HA level were attributed to surface masking or poisoning by HA of secondary Fe(II) mineral growth at or near the clay surface. Our results suggest that HS play an important role in controlling the kinetics and products of Fe(II) precipitation in reducing soils, with effects modulated by soil mineralogy, HS content, and HS properties. Further work is needed to assess the importance of layered Fe(II) hydroxides in natural reducing environments.     

A Comparison of Strategies for Seismic Interferometry

The extraction of the response from field fluctuations excited by random sources has received considerable attention in a variety of different fields. We present three methods for the extraction of the systems response that are based on cross-correlation, deconvolution, and the solution of an integral equation, respectively. For systems that are invariant for time-reversal the correlation method requires random sources on a bounding surface only, but when time-reversal invariance is broken, for example by attenuation, a volume distribution of sources is needed. For this reason the correlation method is not useful for diffusive or strongly attenuating systems. We provide examples of the three methods and compare their merits and drawbacks. We show that the extracted field may satisfy different boundary conditions than does the physical field. This can be used, for example, to suppress surface-related multiples in exploration seismology, to study the coupling of buildings to the subsurface, and to remove the airwave in controlled source electromagnetics (CSEM).