Amplitude Preserving Kirchhoff Migration: A Traveltime Based Strategy

Amplitude versus offset information is a key feature to seismic reservoir characterization. Therefore amplitude preserving migration was developed to obtain this information from seismic reflection data. For complex 3-D media, however, this process is computationally expensive. In this paper we present an efficient traveltime based strategy for amplitude preserving migration of the Kirchhoff type. Its foundations are the generation of traveltime tables using a wavefront-oriented ray-tracing technique, and a generalized moveout relation for 3-D heterogeneous media. All required quantities for the amplitude preserving migration are computed from coarsely gridded traveltime tables. The migration includes the interpolation from the coarsely gridded input traveltimes onto the fine migration grid, the computation of amplitude preserving weight functions, and, optionally, the evaluation of an optimized migration aperture. Since ray tracing is employed for the traveltime computation the input velocity model needs to be smooth, i.e. velocity variations of spatial dimensions below the wavelength of the considered reflection signals are removed. Numerical examples on simple generic models validate the technique and an application to the Marmousi model demonstrates its potential to complex media. The major advantage of the traveltime based strategy consists of its computational efficiency by maintaining sufficient accuracy. Considerable savings in storage space (10⁵ and more for 3-D data with respect to no interpolation at all) can be achieved. The computational time for the stack can be substantially reduced (up to 90% in 3-D) with the optimized migration aperture since only those traces are stacked which really contribute to the image point under consideration.     

Multiwave Amplitude Compensation and Its Sensitivity Analysis to AVO Inversion in Viscoelastic Media

We derive formulae of correction for multi-wave geometric spreading and absorption in lay-ered viscoelastic media,this provides the theoretical foundation for true amplitude compensation of field data and for our sensitivity analysis.The imaging matrix at a plane reflector between viscoelastic media can be determined in the frequency domain using linearized reflection coefficients through Born approxi-mation.We quantiqatively analyze the sensitivity by studying eigenvalues and eigenvectors of the imaging matrix.The results show that two linear combinations of petrophysical parameters can be determined from the multi-wave AVO inversion in the case of amplitude compensation.Multi-wave AVO contains the information of attenuation in the media.However,the sensitivity of multi-wave AVO inversion to at-tenuation is small.     

THE FORMATION MECHANISM OF CONCENTRIC DOUBLE EYE WALL TYPHOON—PART Ⅰ: DYNAMICAL ANALYSIS*

The strong vortex will mutually adjust the thermodynamic field and dynamic field to a state of gradient balance whilst forced by an external cold source,namely,the gradient adjustment process,in which a linearized two-layer model is dealt with in this paper.The analyses show that on account of the heterogeneous radial distribution of the cold source,the adjustment of the thermodynamic and dynamic fields results in the two-peak tangential wind feature,which is analogous to the character of concentric double eyewalls in the strong typhoon.Consequently,the gradient adjustment may be one possible mechanism for the formation of a concentric double-eye typhoon.     

Detection and Evaluation of Gas Hydrates in the Eastern Nankai Trough by Geochemical and Geophysical Methods

Abstract: Interstitial waters extracted from the sediment cores from the exploration wells, “BH‐1” and “MITI Nankai Trough”, drilled ～60 km off Omaezaki Peninsula in the eastern Nankai Trough, were analyzed for the chloride and sulfate concentrations to examine the depth profiles and occurrence of subsurface gas hydrates. Cored intervals from the seafloor to 310 mbsf were divided into Unit 1 (～70 mbsf, predominated by mud), Unit 2 (70–150 mbsf, mud with thin ash beds), Unit 3 (150–250+ mbsf, mud with thin ash and sand), and Unit 4 (275–310 mbsf, predominated by mud). The baseline level for Cl “concentrations was 540 mM, whereas low chloride anomalies (103 to 223 mM) were identified at around 207 mbsf (zone A), 234–240 mbsf (zone B), and 258–265 mbsf (zone C) in Unit 3. Gas hydrate saturation (Sh %) of sediment pores was calculated to be 60 % (zone A) to 80 % (zones B and C) in sands whereas only a few percent in clay and silt. The total amount of gas hydrates in hydrate‐bearing sands was estimated to be 8 to 10 m³ of solid gas hydrate per m², or 1.48 km³ CH₄ per 1 km². High saturation zones (A, B and C) were consistent with anomaly zones recognized in sonic and resistivity logs. 2D and high‐resolution seismic studies revealed two BSRs in the study area. Strong BSRs (BSR‐1) at ～263 mbsf were correlated to the boundary between gas hydrate‐bearing sands (zone C) and the shallower low velocity zone, while the lower BSRs (BSR‐2) at～289 mbsf corresponded to the top of the deeper low velocity zone of the sonic log. Tectonic uplift of the study area is thought to have caused the upward migration of BGHS. That is, BSR‐1 corresponds to the new BGHS and BSR‐2 to the old BGHS. Relic gas hydrates and free gas may survive in the interval between BSR‐1 and BSR‐2, and below BSR‐2, respectively. Direct measurements of the formation temperature for the top 170 m interval yield a geothermal gradient of ～4.3d?C/ 100 m. Extrapolation of this gradient down to the base of gas hydrate stability yields a theoretical BGHS at～230 mbsf, surprisingly ～35 m shallower than the base of gas hydrate‐bearing sands (zone C) and BSR‐1. As with the double BSRs, another tectonic uplift may explain the BGHS at unreasonably shallow depths. Alternatively, linear extrapolation of the geothermal gradient down to the hydrate‐bearing zones may not be appropriate if the gradient changes below the depths that were measured. Recognition of double BSRs (263 and 289 mbsf) and probable new BGHS (～230 mbsf) in the exploration wells implies that the BGHS has gradually migrated upward. Tectonically induced processes are thought to have enhanced dense and massive accumulation of gas hydrate deposits through effective methane recycling and condensation. To test the hypothetical models for the accumulation of gas hydrates in Nankai accretionary prism, we strongly propose to measure the equilibrium temperatures for the entire depth range down to the free gas zone below predicted BGHS and to reconstruct the water depths and uplift history of hydrate‐bearing area.     

Mechanical stratigraphy as a factor controlling the development of a sandbox transfer zone: a three-dimensional analysis

In thrust belts, fold–fault terminations are common features of the structural architecture and can pose complicated problems to unravel, in particular when two or more terminations are in close proximity. Such terminations usually reflect pre-existing attributes. Amongst the many factors, lateral variations in the mechanical stratigraphy can control along-strike geometry and kinematics of fault-related folds.A displacement transfer zone was produced in a compressional sandbox model by means of two adjacent, mechanically different stratigraphic domains. The experiment allowed two discrete chains to develop in the different domains, so that a complex structural setting occurred in the connecting area. Periclinal folds, oblique thrust fronts and oblique ramps developed in the resulting transfer zone. The interaction between periclines in the transfer zone produced lateral culminations in the folded structures. The analysis of displacement across the structural domains revealed that a significant loss of slip along the faults occurred in the relay zone. In this area, imbricate faulting was partially replaced by layer-parallel shortening. A linear relationship appears to exist between the bed length of the thrust sheet and the related fault slip.     

坡度阀值与坡面泥石流--以重庆市北碚区为例

通过对重庆市北碚区的 2 1个坡面泥石流的实地调查表明 ,坡度对坡面泥石流发生具有重要控制性作用。在对坡面泥石流流域大量原始坡度进行统计分析后发现 ,各坡面泥石流的坡度分布具有明显的正态分布特征 ,且平均坡度分布具有显著的规律。根据这一分布规律 ,将本区坡面泥石流发生的坡度阀值确定为三类 :第一类的坡面泥石流植被覆盖率高 ,远离居民区 ,人类影响活动微弱 ,坡度阀值为 32 4° ;第二类坡面泥石流无一例外均位于林地结合部 ,森林覆盖率很低 ,人类活动影响十分剧烈 ,坡度阀值为 2 7 3° ;第三类的两条坡面泥石流相邻 ,均位于观音峡峡谷地带 ,它们的发生完全受地形条件的控制 ,坡度阀值为 36 8°。本文从定量的角度阐述了坡度阀值与坡面泥石流发生的关系 ,为生态脆弱区的划分、工程建设、防灾减灾和政府决策提供服务     

牛顿内摩擦定律和达西定律的改进

根据物体运动的摩擦力定律和物体变形破坏的库仑定律对牛顿内摩擦定律作了改进 ,进而改进了达西定律。用改进的达西定律与实验资料进行拟合 ,阐述了达西定律与实际偏离的原因。     

New heat flow data from the immediate vicinity of the Kola super-deep borehole: Vertical variation in heat flow confirmed and attributed to advection

We present new heat flow values and other geothermal data in the upper crystalline crust in the immediate vicinity of the 12.4-km deep Kola super-deep borehole, NW Russia. Our results show a systematic vertical increase in geothermal gradient and heat flow density as deep as we could measure (1.6 km). Our results confirm earlier results on vertical heat flow trends of in the uppermost part of the Kola super-deep hole, and imply that the thermal regime is not in steady-state conductive conditions. In an area of 3-km × 5-km measurements were performed in 1–2-km deep boreholes surrounding the Kola super-deep hole and on core samples from these holes. Temperature logs are available from 36 holes. Core data exists from 23 boreholes with a total length of 11.5 km at a vertical resolution of 10 m. We carried out a very detailed study on thermal conductivity with regard to anisotropy, inhomogeneity and temperature dependence. Tensor components of thermal conductivity were determined on 1375 core samples from 21 boreholes in 3400 measurements. Additionally, we measured specific heat capacity, heat generation rate, density, porosity, and permeability on selected subsets of core samples. Heat flow from 19 boreholes varies between 31 and 45 mW m⁻² with an average value of 38 mW m⁻². In most boreholes the vertical heat flow profiles show a considerable variation with depth. This is consistent with observations in the upper part of the Kola super-deep borehole. We conclude that this variation is not caused by technical operations but reflects a natural process. It is considered to be due to a combination of advective, structural and paleoclimatic effects. Preliminary 3-D numerical modeling of heat and flow in the study area provides an indication of relative contributions of each of these factors: advective heat transfer turns out to have a major influence on the vertical variation of heat flow, although transient changes in surface temperature may also cause a significant variation. Heterogeneity of the rocks in the study area is less important.     

沾化凹陷Ro分布规律及影响有机质成熟度的因素

结合实际资料，综合确定R0=0．4％为该区生烃门限(含低熟油)，相应的深度门限平均为2250m左右，温度门限平均在90℃左右。造成R。值相对较低的相关因素如下：有机质成分、温度、压力、埋深等。通过比较，沾化凹陷中渤南洼陷和孤北洼陷生烃能力强，而富林洼陷生烃能力相对较弱。     

Effect of sill spacing and sediment size grading on scouring at grade‐control structures

Sequences of arti?cial steps are sometimes used to reproduce the natural step–pool morphology of high‐gradient streams. The depth, length and shape of the scour holes in gravel‐bed rivers can be predicted reasonably using recently developed formulae. However, the properties of the scour holes can sometimes be affected by the distance between structures. This effect is called ‘geometrical interference’ and leads to a reduction of the scour hole compared to its potential size. Geometrical interference may occur in sequences of arti?cial steps in high‐gradient torrents, where structures are sometimes built at distances of a few tens of metres apart, but may also apply to natural step–pool systems. In this paper, a series of tests have been conducted to determine the effect of bed sill spacing and sediment grading on the potential erosion by jets forming over the sills. A new formula is derived, applicable to high‐gradient streams (slope > 0·04), which can be applied to the special case of scour holes developed by interfering sills. Sediment size gradation, not accounted for in previous formulae, is found to have a signi?cant effect on the scour dimensions and is included in the new predictive formula. Copyright © 2004 John Wiley & Sons, Ltd.