Inversion of normal moveout for monoclinic media1

Multiple vertical fracture sets, possibly combined with horizontal fine layering, produce an equivalent medium of monoclinic symmetry with a horizontal symmetry plane. Although monoclinic models may be rather common for fractured formations, they have hardly been used in seismic methods of fracture detection due to the large number of independent elements in the stiffness tensor. Here, we show that multicomponent wide-azimuth reflection data (combined with known vertical velocity or reflector depth) or multi-azimuth walkaway VSP surveys provide enough information to invert for all but one anisotropic parameters of monoclinic media. In order to facilitate the inversion procedure, we introduce a Thomsen-style parametrization for monoclinic media that includes the vertical velocities of the P-wave and one of the split S-waves and a set of dimensionless anisotropic coefficients. Our notation, defined for the coordinate frame associated with the polarization directions of the vertically propagating shear waves, captures the combinations of the stiffnesses responsible for the normal-moveout (NMO) ellipses of all three pure modes. The first group of the anisotropic parameters contains seven coefficients (ε^(1,2), δ^(1,2,3) and γ^(1,2)) analogous to those defined by Tsvankin for the higher-symmetry orthorhombic model. The parameters ε^(1,2), δ^(1,2) and γ^(1,2) are primarily responsible for the pure-mode NMO velocities along the coordinate axes x₁ and x₂ (i.e. in the shear-wave polarization directions). The remaining coefficient δ⁽³⁾ is not constrained by conventional-spread reflection traveltimes in a horizontal monoclinic layer. The second parameter group consists of the newly introduced coefficients ζ^(1,2,3) which control the rotation of the P-, S₁- and S₂-wave NMO ellipses with respect to the horizontal coordinate axes. Misalignment of the P-wave NMO ellipse and shear-wave polarization directions was recently observed on field data by Pérez et al. Our parameter-estimation algorithm, based on NMO equations valid for any strength of the anisotropy, is designed to obtain anisotropic parameters of monoclinic media by inverting the vertical velocities and NMO ellipses of the P-, S₁- and S₂-waves. A Dix-type representation of the NMO velocity of mode-converted waves makes it possible to replace the pure shear modes in reflection surveys with the PS₁- and PS₂-waves. Numerical tests show that our method yields stable estimates of all relevant parameters for both a single layer and a horizontally stratified monoclinic medium.     

GRAVITY GRADIENT TENSORS DUE TO A POLYHEDRON WITH POLYGONAL FACETS1

Using the conjugate complex variables formulation, closed-form formulae for the gravity gradient tensors of the gravitational potential due to a homogeneous polyhedral body composed of polygonal facets are derived. The treatise considers the cases of the observation point being inside the polyhedron, on the surface of a facet, or outside the polyhedron.     

A generalized Biot–Gassmann model for the acoustic properties of shaley sandstones1

We obtain the wave velocities of clay-bearing sandstones as a function of clay content, porosity and frequency. Unlike previous theories, based simply on slowness and/or moduli averaging or two-phase models, we use a Biot-type three-phase theory that considers the existence of two solids (sand grains and clay particles) and a fluid. The theory, which is consistent with the critical porosity concept, uses three free parameters that determine the dependence of the dry-rock moduli of the sand and clay matrices as a function of porosity and clay content.
Testing of the model with laboratory data shows good agreement between predictions and measurements. In addition to a rock physics model that can be useful for petrophysical interpretation of wave velocities obtained from well logs and surface seismic data, the model provides the differential equation for computing synthetic seismograms in inhomogeneous media, from the seismic to the ultrasonic frequency bands.     

Non-equilibrium compaction and abnormal pore-fluid pressures: effects on rock properties1

Knowledge of pore pressure using seismic data will help in planning the drilling process to control potentially dangerous abnormal pressures. Various physical processes cause anomalous pressures on an underground fluid. Non-equilibrium compaction is a significant process of overpressure generation. This occurs when the sedimentation rate is so rapid that the pore fluids do not have a chance to 'escape' from the pore space.
The model assumes a closed system and that the pore space is filled with water and hydrocarbon in a liquid state. Balancing mass and volume fractions yields the fluid pressure versus time of deposition and depth of burial. Thermal effects are taken into account. The pore pressure, together with the confining pressure, determines the effective pressure which, in turn, determines the bulk moduli of the rock matrix.
We assume a sandstone saturated with hydrocarbons and water, for which calibration of the model with experimental data is possible. The seismic velocities and attenuation factors are computed by using Biot's theory of dynamic poroelasticity and the generalized linear solid. The example shows that the formation can be overpressured or underpressured depending on the properties of the saturating fluid. Wave velocities and quality factors decrease with decreasing differential pressure. The effect is important below approximately 20 MPa. The model is in good agreement with experimental data for Berea sandstone and provides a tool for predicting pore pressure from seismic attributes.     

Petrogenetic grid for ultrahigh-pressure metamorphism in the model system CaO-MgO-SiO2-CO2-H2O

Abstract Petrogenetic grids for ultrahigh-pressure (UHP) metamorphism were calculated at different X_co2 conditions in the model system CaO-MgO-SiO₂-CO₂-H₂O involving coesite (Co), diopside (Di), dolomite (Do), enstatite (En), forsterite (Fo), magnesite (Ms), quartz (Qz), talc (Tc), tremolite (Tr) using a published internally consistent thermodynamic data set. Two P-T grids at X_co2= 0.01 and 0.5 are described. In the calculated P-T grid at X_co2= 0.01, four out of 10 stable invariant points, Co-En-Ms-Tc, Co-Di-En-Tc-Tr, Co-Di-Ms-Tc-Tr and Di-En-Ms-Tc-Tr lie within the stability field of coesite. If the fluid phase has X_co2= 0.5, no invariant point is stable under UHP conditions. Some magnesite-bearing assemblages are stabilized by the following three reactions: Di + Ms = Do + Fo + CO₂, Ms + Tr = Do + Fo + CO₂+ H₂O and Ms + Tc = Fo+ CO₂+ H₂O at X_co2= 0.01 and by reaction Ms + Tc = Fo + CO₂+ H₂O together with these three at X_co2= 0.5. Ten possible UHP assemblages for mafic and ultramafic compositions at very low X_co2 conditions include the following: Co-Do-Ms, Co-Di-Ms, Co-Di-Tc, Di-Ms-Tc, Di-En-Tc-, Di-En-Ms, Co-Di-En, Di-En-Fo, Di-Fo-Ms, Di-Do-Fo. Among them, talc-bearing assemblages are restricted to X_co2 < 0.02 and their high-P limit is 31.7 kb (749°C) at X_co2= 0.01. Dolomite-magnesite-silica assemblages have large P-T stability fields even if X_co2 is as low as 0.1, and could occur in cold subduction zones with very low geothermal gradients. Reported UHP coesite-dolomite assemblage is restricted only to a calc-silicate rock interlayered with marble where X_co2 is relatively higher; no such assemblage appears for mafic and ultramafic rocks with low X_co2 evidenced by the occurrence of diopside (or omphacite) at the expense of dolomite + coesite. The effect of X_co2 on the stability of coesite-dolomite-magnesite, diopside-enstatite-magnesite, diopside-talc assemblages is examined and the occurrence of coesite-dolomite, magnesite-bearing and talc-bearing assemblages in the Dabie UHP rocks are interpreted by employing the calculated P-T grids.     

40Ar–39Ar analysis of phlogopite in the Horoman Peridotite Complex, Hokkaido, Japan and implications for its origin

Abstract ⁴⁰Ar–³⁹Ar analysis of phlogopite separated from a plagioclase lherzolite of the Horoman Peridotite Complex, Hokkaido, Japan, has yielded a plateau age of 20.6 ± 0.5 Ma in an environment where the metamorphic fluid was characterized by an almost atmospheric Ar isotopic ratio. The age spectrum is slightly saddle-shaped, implying some incorporation of excess ⁴⁰Ar during the formation of the phlogopite at a depth. As the phlogopite has been inferred to have formed in veins and/or interstitials during exhumation of the peridotite body, metasomatic fluids, to which ground- and sea water might have contributed, were probably involved in the formation of phlogopite in the crustal environment. A total ⁴⁰Ar–³⁹Ar age of 129 Ma of a whole rock sample of the plagioclase lherzolite, from which the phlogopite was separated and is representative of the main lithology of the Horoman Peridotite Complex, indicates the occurrence of excess ⁴⁰Ar. Hence, the age has no geological meaning.     

230Th/234U and 14C dating of a lowstand coral reef beneath the insular shelf off Irabu Island, Ryukyus, southwestern Japan

Abstract   High-resolution seismic reflection profiles delineated the distribution of mound-shaped reflections, which were interpreted as reefs, beneath the insular shelf western off Irabu Island, Ryukyus, southwestern Japan. A sediment core through one of the mounded structures was recovered from the sea floor at a depth of −118.2 m by offshore drilling and was dated by radiometric methods. The lithology and coral fauna of the core indicate that the mounded structure was composed of coral–algal boundstone suggesting a small-scaled coral reef. High-precision α-spectrometric ²³⁰Th/²³⁴U dating coupled with calibrated accelerator mass spectrometric ¹⁴C ages of corals obtained reliable ages of this reef ranging from 22.18 ± 0.63 to 30.47 ± 0.98 ka. This proves that such a submerged reef was formed during the lowstand stage of marine oxygen isotope stages 3–2. The existence of low-Mg calcite in the aragonitic coral skeleton of 22.18 ± 0.63 ka provides evidence that the reef had once been exposed by lowering of the relative sealevel to at least −126 m during the last glacial maximum in the study area. There is no room for doubt that a coral reef grew during the last glacial period on the shelf off Irabu Island of Ryukyus in the subtropical region of western Pacific.     

K–Ar and 40Ar/39Ar phengite ages of Sanbagawa schist clasts from the Kuma Group, central Shikoku, southwest Japan

Conglomerates of the Kuma Group, central Shikoku, southwest Japan contain Sanbagawa schist clasts with a variety of metamorphic grades and lithologies. K–Ar and ⁴⁰Ar/³⁹Ar dating of phengite show all the pelitic schist clasts from low- to high-grade zones have similar phengite ages (82–84 Ma) that are significantly older than those from the in situ Sanbagawa sequence of central Shikoku. This is because the Kuma–Sanbagawa sequence was exhumed earlier than the in situ Asemi sequence with an exhumation process intermediate between those for the Kanto Mountains and the in situ Asemi sequences. ⁴⁰A/³⁹Ar plateau ages (103 and 117 Ma) of phengite in amphibolites indicate the timing of the early stage of the exhumation of the metamorphic pile, probably close to the peak metamorphic age.     

Characterization of dipping fractures in a transverselyisotropic background

Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities.     

40Ar/ 39Ar and K–Ar geochronological age constraints for the inception and early evolution of the Izu–Bonin – Mariana arc system

K–Ar and ⁴⁰Ar/³⁹Ar dates are presented for locations in the Izu–Bonin – Mariana (IBM) forearc (Ocean Drilling Program (ODP) sites 786 & 782, Chichijima, Deep Sea Drilling Program (DSDP) sites 458 & 459, Saipan), and Palau on the remnant arc of the Kyushu–Palau Ridge. For a number of these locations, the ⁴⁰Ar/³⁹Ar plateau and ³⁶Ar/⁴⁰Ar versus ³⁹Ar/⁴⁰Ar isochrons give older ages than the K–Ar results. The most important results are: (i) at site 786, initial construction of the proto-IBM (now forearc) basement occurred at least by ca 47–45 Ma, consistent with the age of the immediately overlying sediments (middle Eocene nannofossil Zone CP13c); the younger pulse of construction dated at ca 35 Ma by K–Ar could not be confirmed by ⁴⁰Ar/³⁹Ar analysis; (ii) ⁴⁰Ar/³⁹Ar ages for the initial construction of the Mariana portion of the IBM system are as old as those of the Izu–Bonin portion, for example at site 458, initial construction commenced at least by ca 49 Ma and at ca 47 Ma at Saipan (Sankakayuma Formation); and (iii) a combination of K–Ar and ⁴⁰Ar/³⁹Ar ages indicate continued boninite magmatism in the Izu–Bonin forearc (and remnant arc at Palau) until ca 35 Ma. Subduction inception including boninite series rocks along most of the exposed length of the IBM system, clearly preceded by some 5 million years the Middle Eocene (ca 43.5 Ma) change in Pacific plate motion. Boninitic series magmatism persisted at locations now exposed in the forearc for ～ 15 million years after arc inception concurrently with low-K tholeiitic series eruptions from a subaerial arc system, established at ≥ 40 Ma, on the Kyushu–Palau Ridge. For the Mariana portion of the IBM system, reconstruction of the proto-arc places this activity adjacent to the concurrent but orthogonally spreading Central Basin Ridge of the West Philippine Basin. It is possible that a combination of subduction of a young North New Guinea Plate beneath newly created back-arc basin crust may account for some of the features of the Mariana system. It is clear, however, that the understanding of the processes of subduction initiation and early IBM arc development is incomplete.     

Stable Isotopes (δ18O, δ2H) of Pore Waters in Clay-Rich Aquitards: A Comparison and Evaluation of Measurement Techniques

Pore water collected from piezometers installed in a thick clay-rich till were used to compare and evaluate four techniques for obtaining δD and δ¹⁸O values in these media. The techniques included mechanical squeezing, centrifugation, azeotropic distillation, and a direct soil-water equilibration technique. Direct CO₂-core equilibration yielded sufficiently accurate and reproducible δ¹⁸O results of pore water in clay-rich tills. In addition, this method eliminated the need for labor-intensive complete extraction of water from the geologic media. Mechanical squeezing and centrifugation produced results similar to direct equilibration. However, both of these methods exhibited a greater degree of variability and were laborious and more time consuming. Small differences in δ¹⁸O values between piezometer water and equilibrated, squeezed, and centhfuged samples suggested that each method collected different fractions of the clay-water reservoir. Although these subtle differences were not conclusive, they did suggest the presence of weakly bound water and highlighted the difference between these three techniques for determining the stable isotopic composition of pore water in clay-rich aquitards. Azeotropic distillation produced a high level of discrepancy in δD andδ¹⁸O results compared to the other methods. Incomplete extraction was considered the most probable cause of this error. The results of this study suggested that direct equilibration is the best method for determining detailed δD and δ¹⁸O values of pore water in clay-rich aquitards.     

微囊藻水华的资源化利用:吸附重金属离子Cu2+、Cd2+和Ni2+的实验研究

为将捞取的微囊藻水华资源化,用作重金属生物吸附剂,本文研究了在不同金属离子浓度和不同pH条件下微囊藻水华对重金属Cu²⁺、Cd²⁺和Ni²⁺的吸附效果.研究结果显示,重金属离子浓度分别为20.00、60.00、100.00和140.00μg/ml时,150 mg干重微囊藻水华对Cu²⁺的去除率分别为67.79%、37.47%、35.93%和34.70%,对Cd²⁺的去除率分别为73.31%、65.87%、60.09%和42.37%,对Ni²⁺的去除率分别为47.89%、31.87%、28.46%和21.61%.在较低的重金属浓度(20.00μg/ml)下,微囊藻对3种重金属的去除率最高.但从单位微囊藻生物量吸附重金属量来看,金属离子初始浓度越高,吸附重金属的总量越高;在相同金属浓度下,微囊藻水华对3种金属的吸附效率:Cd²⁺ > Cu²⁺ > Ni²⁺.微囊藻吸附Cu²⁺的最适pH为5.0,Cd²⁺和Ni²⁺均为6.0,20.00μg/ml浓度时去除率分别达到了71.90%、85.67%和55.43%.微囊藻裂解释放的可溶性物质对重金属吸附影响不大.研究结果显示微囊藻水华可作为重金属吸附剂有效地运用于重金属污水处理.     

40Ar/39Ar ages of blueschist facies pelitic schists from Qingshuigou in the Northern Qilian Mountains, western China

Abstract   Pelitic schists from Qingshuigou in the Northern Qilian Mountains of China contain mainly glaucophane, garnet, white mica, clinozoisite, chlorite and piemontite. Isotopic age dating of these schists provides new constraints on the formation of the high-grade blueschists at Qingshuigou. White mica ⁴⁰Ar/³⁹Ar ages range from 442.1 to 447.5 Ma (total fusion age of single grain) and from 445.7 to 453.9 Ma (integrated age of white mica concentrates). These ages (442.1–453.9 Ma) represent the peak metamorphic ages or cooling ages of the blueschists during exhumation shortly after peak metamorphism. The ⁴⁰Ar/³⁹Ar dates in the present study are similar to ages previously reported for eclogites and blueschists in the area; this suggests that both the eclogites and pelitic sediments underwent high-grade metamorphism during the same subduction event. From this chronological evidence and the presence of well-developed Silurian remnant-sea flysch and Devonian molasse, it is concluded that the Northern Qilian Ocean had closed by the end of the Ordovician, and rapid orogenic uplift followed in the Devonian.