密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

海洋声层析的基本原理和应用

从射线声学和简正波声学的角度，概述了海洋声层析的基本理论，包括射线走时反演、简正波走时反演、简正波相位反演和简正波水平折射层析。海洋声层析以反演海水温度和流速为基础。还总结了声层析在海洋学研究中的应用。     

Monitoring the Kuroshio Extension with Dynamically Constrained Synthesis of the Acoustic Tomography, Satellite Altimeter and in situ Data

A finite-difference quasigeostrophic (QG) model of an open ocean region has been employed to produce a dynamically constrained synthesis of acoustic tomography and satellite altimetry data with in situ observations. The assimilation algorithm is based upon the 4D variational data interpolation scheme controlled by the model's initial and boundary conditions. The data sets analyzed include direct and differential travel times measured at the array of five acoustic transceivers deployed by JAMSTEC in the region of the Kuroshio Extension in 1997, Topex/Poseidon altimetry, CTD soundings, and ADCP velocity profiles. The region monitored is located within the area 27.5°–36.5°N, 143°–155°. The results of assimilation show that mesoscale variability can be effectively reconstructed by five transceivers measuring direct and reciprocal travel times supported by relatively sparse in situ measurements. The misfits between model and data lie within the observational error bars for all the data types used in assimilation. We have compared the results of assimilation with the statistical inversion of travel time data and analyzed energy balances of the optimized model solution. Energy exchange between the depth-averaged and shear components of the observed currents reveals a weak decay of the barotropic mode at the rate of 0.2 ± 0.7⋅10⁻⁵ cm²/s³ due to topographic interaction. Mean currents in the region are unstable with an estimate of the available potential energy flux from the mean current to the eddies of 4.7 ± 2.3⋅10⁻⁵ cm²/s³. Kinetic energy transition has the same sign and is estimated as 2.8 ± 2.5⋅10⁻⁵ cm²/s³. Potential enstrophy is transferred to the mesoscale at a rate of 5.5 ± 2.7⋅10⁻¹⁸ s⁻³. These figures provide observational evidence of the properties of free geostrophic turbulence which were predicted by theory and observed in numerical experiments. This revised version was published online in July 2006 with corrections to the Cover Date.     

垦71区块井间地震资料处理研究与应用

由于避开了地表低降速带和风化层的影响,井间地震能够采集到高频地震信号,经过资料处理可以得到较高分辨率的速度和上行反射波成像剖面。2004～2005年,胜利油田在垦71 区块实施了8对井间地震资料采集。针对原始资料特点,研究了以层析成像、管波衰减、VSP-CDP技术和有限角度叠加成像等关键处理技术,较好地完成了该区块的井间地震资料处理。初步分析解释表明,井间地震可以识别低序级断层、小幅度构造和薄储层。     

New basis for the constitutive modelling of aggregated soils

Natural and compacted soils are usually characterized by aggregation of particles. The mechanical behaviour of these materials depends on soil structure. The oedometric compression tests performed on aggregated samples presented here showed that these materials exhibit a yield limit depending not only on stress history and stress state but also on soil structure. Evidence is provided using the neutron tomography technique. These results revealed that soil structure modification occurs together with plastic deformations. The experimental results are used to propose a new state parameter to quantify the soil structure. Based on pore-scale experimental observations, an evolution law for this parameter is proposed as a function of associated plastic strains. Considering both soil fabric and inter-particle bonding effects, a new yield limit depending on stress state, stress history and soil structure is introduced for the aggregated soils. Accordingly, a new constitutive framework consistent with strain hardening plasticity is proposed to consider soil structure effects in the modelling of aggregated soils.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

Mantle plumes from top to bottom

Hotspots include midplate features like Hawaii and on-axis features like Iceland. Mantle plumes are a well-posed hypothesis for their formation. Starting plume heads provide an explanation of brief episodes of flood basalts, mafic intrusions, and radial dike swarms. Yet the essence of the hypothesis hides deep in the mantle. Tests independent of surface geology and geochemistry to date have been at best tantalizing. It is productive to bare the current ignorance, rather than to dump the plume hypothesis. One finds potentially fruitful lines of inquiry using simple dynamics and observations. Ancient lithospheric xenoliths may reveal heating by plumes and subsequent thermal equilibration in the past. The effect at the base of the chemical layer is modest 50-100 K for transient heating by plume heads. Thinning of nonbuoyant platform lithosphere is readily observed but not directly attributable to plumes. The plume history in Antarctica is ill constrained because of poor geological exposure. This locality provides a worst case on what is known about surface evidence of hotspots. Direct detection of plume tail conduits in the mid-mantle is now at the edge of seismic resolution. Seismology does not provide adequate resolution of the deep mantle. We do not know the extent of a chemically dense dregs layer or whether superplume regions are cooler or hotter than an adiabat in equilibrium with the asthenosphere. Overall, mid-mantle seismology is most likely to give definitive results as plume conduits are the guts of the dynamic hypothesis. Finding them would bring unresolved deep and shallow processes into place.