叠前逆时深度偏移中的激发时间成像条件

与其他偏移方法相比,逆时偏移基于精确的波动方程而不是对其近似,用时间外推来代替深度外推.因此,它具有良好的精度,不受地下构造倾角和介质横向速度变化的限制.激发时间成像条件的求取是叠前逆时偏移的难点之一,本文采用求解程函方程的方法得到地下各点的初至波走时,以此作为叠前逆时偏移的成像条件.基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法精度较高并适用于强纵横向变速的复杂介质.试算结果表明,在复杂介质模型中利用叠前逆时深度偏移收到了很好的成像效果.     

两种类型共炮点域相关型叠前深度偏移方法分析

研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。     

盐丘模型 推覆体模型 波动方程 叠前深度偏移

三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节.由于其技术实现不仅涉及波场延拓理论的创新，而且需要大规模计算，因而研究难度较大. 本文以实验效果的取得为目的，完整地实现了SEG/EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算. 文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力. 具体包括：1）对于盐丘模型，文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较；2）对推覆体模型，文中进行了脉冲响应测试；3）由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力.     

用GPU提速地震资料单程波有限差分叠前深度偏移(英文)

复杂介质情况下,地震波延多路径传播,此时基于波动方程延拓的深度成像方法,相对于Kirchhoff方法能够获得更为精确的成像效果,但是,该深度偏移方法由于高昂的计算消耗阻碍了它在生产中的应用。譬如,叠前深度偏移计算需要大规模的计算机集群,占地面积和电能消耗大。本文介绍了应用一种新的GPU计算架构来辅助CPU进行偏移计算。基于新架构的波动方程深度偏移提高了计算效率,而且机器占地面积和电能消耗也大幅度减少。本文以有限差分波动方程深度偏移为例,介绍了其编程模型和程序优化环节,提高了深度偏移计算效率。2D和3D测试表明,与相同单位个AMD2.5GHz CPU计算相比,该架构下的有限差分波动方程叠前深度偏移计算效率提高至少35倍。     

任意复杂介质中主能量法地震波走时计算

积分法叠前深度偏移及层析成像的核心是复杂介质情况下的地震波走时计算. 复杂构造的高精度地震成像需要有稳健的走时计算方法。本文把 Nichols提出的用地震波主能量计算走时的方法由二维推广到三维,并推导出三维波动方程Helmholtz形式在球坐标系下用因式分解法求解的差分表达式.三维SEG/EAGE盐丘模型的理论走时计算和积分法叠前深度偏移的实践都验证了本文方法的正确性.     

2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

面炮成像、控制照明与AVA道集

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点.面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题.本文采用面炮成像技术进行叠前深度偏移,通过对面炮震源下行波场的质量控制和优选射线个数和范围,以达到最佳的成像效果.采用控制照明技术,较大地提高了目标地层的成像精度.与此同时,得到振幅随入射角变化(AVA)道集,有利于叠前振幅解释和储层岩性预测.数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算.     

任意广角波动方程叠前逆时深度偏移

从任意广角声波方程出发,在深度与时间方向采用二阶差分代替微分,在水平方向上采用高阶差分代替微分,得到任意广角声波方程逆时延拓的有限差分格式. 通过差分求解程函方程得到叠前逆时深度偏移的成像条件,在此基础上实现了二维各向同性介质中任意广角方程的叠前逆时深度偏移. 对理论模型的偏移处理表明,该方法能够有效地消除波场下行传播过程中的层间反射,减小偏移噪声. 同时,该方法还具有保真性高、 适用于大倾角地层的偏移等优点.      

深水崎岖海底地震数据成像方法与应用

本文针对深水崎岖海底地震数据深部成像困难的问题,讨论了波动方程叠前深度域保幅偏移的基本理论,利用数值模型验证了波动方程叠前深度域保幅偏移算法的保持振幅特性,给出了该方法与非保幅偏移方法对于复杂构造成像精度的对比,证明了保幅偏移方法可以提高复杂构造的成像精度.本文还将波动方程叠前深度域保幅偏移算法应用到实际资料的处理中,处理结果表明该方法能够有效消除崎岖海底对深部地层的影响.     

双复杂介质条件下频率空间域有限差分法保幅偏移

油气勘探的重点正转向复杂地表条件和复杂地质条件的区域.双复杂条件下的叠前深度偏移是解决复杂地表条件和复杂地质构造成像的有效手段.基于“逐步累加”的“直接下延”法是解决复杂地表成像的有效手段,能够较好地消除地形起伏的影响.波动方程频率空间域有限差分(xwfd)叠前深度偏移对介质速度横向变化有较强的适应性,适宜于复杂构造的偏移成像,同其他常规波动方程深度偏移一样,常规的xwfd偏移方法,主要也是针对相位进行波场延拓,没有对振幅做任何处理.我们基于保幅单程波方程,推导出了基于xwfd的保幅波场延拓算子,针对xwfd求解时引入误差的影响,我们在xwfd保幅波场延拓过程中加入了误差补偿,实现了带误差补偿的xwfd保幅偏移.基于带误差补偿的xwfd保幅算子,应用适合起伏地表的直接下延法,对双复杂介质模型和实际资料进行了试算,改善了双复杂介质的成像效果.其中,误差补偿可以在若干个外推步长上进行,所以相对于保幅傅里叶有限差分(ffd)法偏移来说,该方法在改善成像质量的同时,也具有较高的运算效率.     

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

正~~     

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.     

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.     

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.