轴对称各向异性介质波动方程深度偏移的对称非平稳相移方法

由所建立的三维qP波相速度表示式出发，导出并解析求解各向异性介质中的频散方程，得到三维各向异性介质中的相移算子，进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质，发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致，将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化；文中给出的算例虽然是针对二维VTI介质的，但所提出的算法同样适用于三维TI介质.     

大步长波场深度延拓的理论

波场延拓是地震偏移成像的基础. 快速进行目标区波场延拓对石油勘探中急需发展的深部地震勘探和无组合海量地震数据的成像有重要意义. 在目标区成像中，目前已有的波场延拓方法，包括基于走时计算的Dix方法和射线追踪方法，以及基于小步长波场递推的方法，在适应复杂介质、计算精度和计算效率的某一方面还不能完全满足实际需要. 本文提出一种基于“算子相位”李代数积分的快速计算延拓算子的方法，称为大步长波场延拓方法. 在该方法中，指向目标区的波场延拓算子象征的复相位被表示成波数的线性组合. 线性组合的系数是层速度函数及其导数的深度积分，计算和存储较为方便. 波场延拓算子通过相移算子加校正的方法，利用快速Fourier变换在空间域和波数域予以实现. 利用动力学等价关系导出了便于计算的表达式. 本文比较了算子主象征函数用一步法展开和用两步法展开的精度，从而说明大步长方法的精度要高于递推方法. 在横向和纵向线性变化介质中，将大步长方法的脉冲响应与递推法做了比较，说明大步长延拓算子的走时精度主要取决于相移因子中的横向变速校正项；且在各种近似下，大步长算子发生的频散都非常小.     

时-空局域化地震波传播方法:Dreamlet叠前深度偏移

提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.     

波场延拓短算子构造方法

在频率-空间域显式叠前深度偏移中，波场深度延拓是通过显 式差分短算子与波场的空间褶积完成的. 基于对显式差分短算子的设计方法的研究，提出了 一种基于相移算子约束的离散光滑插值的构造一维显式短算子的方法. 通过离散光滑插值法 ，在频率-波数域中，以传播区内的相移算子为约束，在传播区外的算子两端处以零点为约 束，进行离散光滑插值，使得所得算子具有二阶光滑可导性，则其对应的频率-空间域中的 算子就可以取得很短. 该方法设计简单，精度高，能够满足波场深度延拓的需要.     

基于多参量的模拟退火全局优化傅里叶有限差分算子

傅里叶有限差分法(FFD)能够处理复杂地质构造中的波传播问题,但对陡倾角成像仍有明显的误差.优化参数的方法能够在保持计算效率的前提下进一步提高陡倾角的成像精度.本文在有理近似的基础上,将FFD算子展开式中的常系数由两个拓展为四个,然后采用模拟退火算法对这四个参数进行全局优化.本方法除了考虑速度对比度以外,还考虑了频率和延拓步长等参量的影响.理论误差分析和脉冲响应测试均表明该方法能极大地提高FFD算子的精确传播角度.二维SEG/EAGE盐丘模型实验表明本文方法对陡倾角以及盐下构造的成像精度明显高于未优化的FFD法.将本文的方法与交替方向加插值的方法结合应用于三维脉冲响应测试更进一步证实了本文方法的有效性.     

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀...     

基于局部余弦基小波束的观测系统沉降法叠前深度偏移

将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点.     

波动方程深度偏移的频率相关变步长延拓方法

发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据（Marmousi模型）和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法.     

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

裂步法最小二乘偏移

最小二乘偏移的核心思想是用正向传播算子的逆算子代替其共轭转置算子进行偏移.它克服了常规偏移方法用共轭转置算子偏移的缺点.本文详细地描述了裂步法最小二乘偏移的原理,实现算法,并通过两个模型展示了裂步法最小二乘偏移的可行性和有效性.并且通过研究发现最小二乘较容易被扩展到基于其它波场延拓算子的最小二乘偏移.     

SCIENCE CHINA Earth Sciences SUBJECT INDEX TO VOLUME 57(2014)

正~~     

Distribution of throughfall and stemflow in multi‐strata agroforestry,perennial monoculture,fallow and primary forest in central Amazonia,Brazil

The partitioning of rain water into throughfall, stemflow and interception loss when passing through plant canopies depends on properties of the respective plant species, such as leaf area and branch angles. In heterogeneous vegetation, such as tropical forest or polycultural systems, the presence of different plant species may consequently result in a mosaic of situations with respect to quantity and quality of water inputs into the soil. As these processes influence not only the water availability for the plants, but also water infiltration and nutrient leaching, the understanding of plant effects on the repartitioning of rain water may help in the optimization of land use systems and management practices. We measured throughfall and stemflow in a perennial polyculture (multi‐strata agroforestry), monocultures of peach palm (Bactris gasipaes) for fruit and for palmito, a monoculture of cupuaçu (Theobroma grandiflorum), spontaneous fallow and primary forest during one year in central Amazonia, Brazil. The effect on rain water partitioning was measured separately for four useful tree species in the polyculture and for two tree species in the primary forest. Throughfall at two stem distances, and stemflow, differed significantly between tree species, resulting in pronounced spatial patterns of water input into the soil in the polyculture system. For two tree species, peach palm for fruit (Bactris gasipaes) and Brazil nut trees (Bertholletia excelsa), the water input into the soil near the stem was significantly higher than the open‐area rainfall. This could lead to increased nutrient leaching when fertilizer is applied close to the stem of these trees. In the primary forest, such spatial patterns could also be detected, with significantly higher water input near a palm (Oenocarpus bacaba) than near a dicotyledonous tree species (Eschweilera sp.). Interception losses were 6·4% in the polyculture, 13·9 and 12·3% in the peach palm monocultures for fruit and for palmito, respectively, 0·5% in the cupuaçu monoculture and 3·1% in the fallow. With more than 20% of the open‐area rainfall, the highest stemflow contributions to the water input into the soil were measured in the palm monocultures and in the fallow. Copyright © 1999 John Wiley & Sons, Ltd.     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Anthropogenic Eutrophication and Red Tide Outbreak in Lacustrine Systems of the Kashmir Himalaya

Red tide, a recurrent phenomenon has become conspicuous in several Kashmir lake ecosystems since 1991. The responsible organism (Euglena pedunculata), a rare flagellate rediscovered in the Kashmir Himalaya (Khan 1993) caused first and unprecedented red tide outbreak, constituting a maximum of 96% of resident numerical phytoplankton density in Dal Lake. At present, conflicting hypotheses exist on the generation of causal assemblage(s) imparting redness to waters: Jeeji Bai (1991) linked its origin to acid precipitation – a fallout of burning oil‐fields during the Gulf War – whilst Khan (1993) holds local factor(s) responsible. Field/experimental studies support the latter contention that the influx of untreated sewage, in unison with warm temperatures, high levels of PhAR, iron and interruption to hydrological flow‐pattern together with absence/or reduction in grazing activity created conducive environmental milieu for red tide outbreak. Dal Lake “red tide” drifted the bloom‐inoculum to other waters, including Lake Wular, where additional ecological niches were carved out, threatening the aesthetic value and biological diversity of Kashmir lakes. Ecological monitoring indicates frequent seasonal red tide occurrence in Dal Lake (including summer‐autumn event of 1998) which testifies its unabated eutrophication status. Further studies are needed on ecological adaptability and biogeographic distribution of this rare and unique red tide‐causing flagellate.