叠前地震反演技术的进展及其在岩性油气藏勘探中的应用

与叠后地震反演相比,叠前地震反演可得到更丰富的储层信息,可提高储层的描述精度.目前,叠前地震反演主要包括弹性阻抗反演、叠前P波阻抗和S波阻抗联合反演、叠前地震波形反演.文中概述了叠前地震反演各项技术的主要进展及其典型应用实例,提出了叠前地震反演技术在岩性油气藏勘探开发中的应用策略.叠前地震反演技术在未来的油气勘探开发中将有广阔的应用前景.     

长江口南汇湿地植被的光谱吸收特征研究

植被的光谱吸收特征与植被生长状况密切相关,同时也受土壤湿度等外界因素影响。本文以长江口南汇湿地典型植被芦苇、互花米草和海三棱藨草为研究对象,在计算其光谱特征参数的基础上,分析了这3种植被的光谱吸收特征形变PSDI值。最后结合实测的土壤湿度,分析了3种植被的PSDI值和土壤湿度的相关性。研究结果发现:互花米草的PSDI值最大,海三棱藨草的PSDI值最小,芦苇介于两者之间,表明互花米草受外界的环境影响较大;互花米草的PSDI值与土壤湿度的相关性最大,芦苇的PSDI值与土壤湿度的相关性最小,海三棱藨草介于两者之间,表明互花米草适应湿地环境能力较强。     

基于遥感光谱反射率反演黄河口海域Ⅱ类水体悬浮泥沙浓度模型

依据黄河口海域含沙水体野外实测光谱反射率为参考实验数据,结合水体表面光谱观测数据建立了经验回归模型,从而定量估测黄河口海域水体表层悬浮泥沙质量浓度,最终得出黄河口海域不同年份水体的悬浮泥沙质量浓度图,这对分析河口海岸带冲淤变化过程具有重要的意义。     

用H/V谱比法计算云南区域数字地震台站的场地响应

用H/V谱比法,计算给出了云南区域数字地震台网23个子台S波随频率变化的台站场地响应特征。结果表明:在1～10Hz频段内,23个子台S波的场地响应相对较平坦,在1.41～2.91之间变化,平均2.08;在大于10Hz的高频段,部分台站的场地响应有较明显的放大。     

The control of maturity and kerogen type on quantitative analytical pyrolysis data

Quantitative pyrolysis-gas chromatography has been performed on 96 kerogen samples isolated from 17 wells on the Norwegian Continental shelf. Petrographic and bulk geochemical measurements were also performed on the samples, and a combined data set of 117 variables for each sample was analysed using principal components analysis (PCA). This approach provides an objective and reproducible means of kerogen characterisation, which can be easily automated. In addition to objective kerogen characterisation and facile visualisation of facies and maturity related chemical trends, the method has the potential to allow objective prediction of key geochemical parameters such as maturity level from pyrogram data.     

Attenuation of compressional waves in peridotite measured as a function of temperature at 200 MPa

A technique has been developed to determine attenuation in rocks at high temperature using a gas-media, high-pressure apparatus. A pulse transmission technique and a spectral ratio method are used to study compressional seismic properties of rocks. Seismic waves are transmitted to and from the sample through buffer rods of mullite. The effect of seismic wave reflections within the sample assembly are cancelled out by taking ratios of the spectra measured at different temperatures. In order to obtain good signal-to-noise ratio for resolving the attenuation at high pressure and temperature, special care is taken in the sample assembly and the ultrasonic coupling between the sample, buffer rods and transducers. A very tight connection of the sample-buffer rod-transducer is essential for obtaining high frequency signals (>300 kHz) at high temperature. A small mass is attached to each outside end of the transducer to drive low frequency signals (<250 kHz) into the sample. Before attenuation measurements, the sample and the buffer rods are tightly compacted in a platinum tube at high pressure and room temperature to ensure pressure seal of the sample assembly. The frequency range of measurement covers 50 to 450 kHz for the sample. Attenuation is very small in the buffer rod compared to the sample for the entire temperature range of the study. Because of the small attenuation, a wide frequency band of 50 kHz to 3.2 MHz can be covered for investigating the attenuation in the buffer rod. The technique has been used to measure attenuation at high confining pressure, and temperatures including sub- and hyper-solidus of upper mantle rocks. Therefore, effects of partial melting on attenuation can be studied.The method is applied to the attenuation measurement in a peridotite as a function of temperature to 1225°C at 200 MPa confining pressure. At high temperature, signal amplitude decays more rapidly at high frequency than at low frequency, from which attenuation (andQ) can be determined using a spectral ratio method. No frequency dependence ofQ is resolved for both the sample and the buffer rod over the entire temperature and frequency ranges of the measurement. The results show thatQ decreases rapidly with increasing temperature even in the temperature range below the solidus of peridotites. Such temperature sensitivity ofQ is probably more useful to probe thermal structure in the upper mantle than that of conductivity at temperatures below the solidus. The results in this study are compared with available seismic velocity, electrical conductivity and solidus data for peridotites, suggesting that there is no discontinuous change in both mechanical and electrical properties of peridotites at the solidus temperature. Even at hypersolidus temperatures, it appears that velocity drops and conductivity increases continuously (not abruptly) with increasing melt fraction. This implies that mechanical and electrical properties of the upper mantle will gradually change at the boundary where the geotherm crosses the solidus.     

Global spectra of energy and enstrophy and their fluxes during July 1979

Transient and stationary spectra of kinetic energy (KE), available potential energy (APE) and enstrophy (EN), and their spectral fluxes as a function of the two-dimensional wavenumbern were computed for July 1979. Triangular truncation at zonal wavenumber 42 was used for computation. The slopes of various spectra in the wavenumber range 14≤n≤25 were obtained by fitting a straight line in log-log scale by the least square method. The transientKE, APE andEN spectra in the lower (upper) troposphere had slopes −2·21 (−2·30), −2·65 (−2·64) and −0·36 (−0·46), respectively. The effect of stationary and divergent motion on the slope values was investigated. The possible correlation between the slope and percentage of transient component in the combined energy and enstrophy was examined to identify the transient motion of the atmosphere with the two-dimensional homogeneous isotropic turbulence. The vertically averaged slope of kinetic energy and enstrophy in the lower (upper) troposphere was close to the value at 700 (200) hPa level. The spectral fluxes of kinetic energy and enstrophy in the wavenumber range 14≤n≤25 satisfied, to a very rough approximation, the criteria of inertial subrange. The stationary fluxes were small. The estimated stationary-transient component of flux was larger, comparable and less than the corresponding transient flux of APE, KE and EN. Representative levels for computation of energy and enstrophy spectra and their fluxes in the lower and upper troposphere were identified.     

The Fourier Integral Method: An efficient spectral method for simulation of random fields

The Fourier Integral Method (FIM) of spectral simulation, adapted to generate realizations of a random function in one, two, or three dimensions, is shown to be an efficient technique of non-conditional geostatistical simulation. The main contribution is the use of the fast Fourier transform for both numerical calculus of the density spectral function and as generator of random finite multidimensional sequences with imposed covariance. Results obtained with the FIM are compared with those obtained by other classic methods: Shinozuka and Jan Method in 1D and Turning Bands Method in 2D and 3D, the points for and against different methodologies are discussed. Moreover, with the FIM the simulation of nested structures, one of which can be a nugget effect and the simulation of both zonal and geometric anisotropy is straightforward. All steps taken to implement the FIM methodology are discussed.     

用时间序列模拟和预报天山乌鲁木齐河源1号冰川年径流总量

本文把1958—1984年天山乌鲁木齐河源1号冰川的融水径流总量的时间序列X(t),分解为趋势项L(t)、周期项P(t)、平稳项S(t)、随机项ε(t),使该时间序列表示为X(t)=L(t)+P(t)+S(t)+ε(t)。采用非线性回归提取L(t),用谱分析和Fourier级数提取P(t),余差用自回归方程建模,用上述项的叠加作出预报,按相关指数公式计算R=0.90,效果良好。     

A univariate model for long-term streamflow forecasting

This paper, the first in a series of two, employs the principle of maximum entropy (POME) via maximum entropy spectral analysis (MESA) to develop a univariate model for long-term streamflow forecasting. Three cases of streamflow forecasting are investigated: forward forecasting, backward forecasting (or reconstruction) and intermittent forecasting (or filling in missing records). Application of the model is discussed in the second paper.