高分辨率非线性三维整体反演研究(英文)

高分辨率非线性三维整体反演方法是基于非线性理论,在层位控制下,将工区多井(或全部井)的测井数据与井旁地震道数据输入具有多输入多输出的网络,同时进行整体训练,可获得整个工区的自适应权函数,并建立综合非线性映射关系,并根据储层在纵横方向上的地质变化特征更新这种非线性映射关系,这样,就能对反演过程及其反演结果起到约束和控制的作用,从而获得稳定且分辨率高的地震反演剖面(速度反演剖面/波阻抗反演剖面/密度反演剖面),实现整体反演,该方法通过模型试算和实际资料处理,获得较好的地质效果,证明该方法精度高、实用性强,可用于储层的定量分析。     

INVERSION OF AN UNMIGRATED STACKED SECTION TO DETERMINE AN INTERVAL VELOCITY MODEL*

A seismic inversion procedure is developed that inverts data available from an unmigrated stacked section to produce an interval velocity model. It attempts to overcome some of the limitations of existing methods by using a generalized linear inversion technique. The inversion process incorporates several features: (i) Lateral interval velocity variations are permitted, (ii) A fast accurate forward model was developed, (iii) Input data is weighted according to the accuracy with which it has been acquired. The procedure is applied to seismic data from the Gippsland Basin, an area offshore South-East Australia.     

Look-ahead vertical seismic profiling inversion approach for density and compressional wave velocity in Bayesian framework

To reduce drilling uncertainties, zero-offset vertical seismic profiles can be inverted to quantify acoustic properties ahead of the bit. In this work, we propose an approach to invert vertical seismic profile corridor stacks in Bayesian framework for look-ahead prediction. The implemented approach helps to successfully predict density and compressional wave velocity using prior knowledge from drilled interval. Hence, this information can be used to monitor reservoir depth as well as quantifying high-pressure zones, which enables taking the correct decision during drilling. The inversion algorithm uses Gauss–Newton as an optimization tool, which requires the calculation of the sensitivity matrix of trace samples with respect to model parameters. Gauss–Newton has quadratic rate of convergence, which can speed up the inversion process. Moreover, geo-statistical analysis has been used to efficiently utilize prior information supplied to the inversion process. The algorithm has been tested on synthetic and field cases. For the field case, a zero-offset vertical seismic profile data taken from an offshore well were used as input to the inversion algorithm. Well logs acquired after drilling the prediction section was used to validate the inversion results. The results from the synthetic case applications were encouraging to accurately predict compressional wave velocity and density from just a constant prior model. The field case application shows the strength of our proposed approach in inverting vertical seismic profile data to obtain density and compressional wave velocity ahead of a bit with reasonable accuracy. Unlike the commonly used vertical seismic profile inversion approach for acoustic impedance using simple error to represent the prior covariance matrix, this work shows the importance of inverting for both density and compressional wave velocity using geo-statistical knowledge of density and compressional wave velocity from the drilled section to quantify the prior covariance matrix required during Bayesian inversion.     

Shear wave velocity modelling in crustal rock for seismic hazard analysis

P-wave velocity data along with the thickness of sedimentary and crystalline layers within bedrock were collected from all global regions and presented in the Global Crustal Model CRUST2.0, published in 2001. This well-organised database provides invaluable potential contributions towards future seismic hazard modelling, particularly for stable continental regions (SCRs), where there is a scarcity of representative strong motion records for conventional modelling purposes. The P-wave velocity information presented in CRUST2.0 has been converted herein to S-wave velocity information. The latter is especially important for purposes of seismic hazard modelling. The value of the CRUST2.0 model has therefore been greatly enhanced by the important findings presented and further developed in this paper. By making the best use of available information on crustal conditions, the amplification behaviour of seismic waves affecting a region, an area or a site for any given earthquake scenario may be predicted. The developed methodology, which is intended for worldwide applications, has been illustrated by case studies in which model S-wave velocity profiles were developed for different geological regions within North America. The model profiles were found to be in excellent agreement with field measurements reported for each respective region.     

Effects of shield brine on the safe disposal of waste in deep geologic environments

The salinity of groundwater increases with depth in the Canadian Shield (up to 1.3 kg/L of density). The existence of brine can be critically important for the safe geologic disposal of radioactive wastes, as dense brine can significantly retard the upward migration of radionuclides released from repositories. Static and flushing conditions of the deep brine are analyzed using a U-tube analogy model. Velocity reduction due to the presence of dense brine is derived under flushing conditions. A set of illustrative numerical simulations in a two-dimensional cross section is presented to demonstrate that dense brine can significantly influence regional groundwater flow patterns in a shield environment. It is implied from the results that (1) the existence of Shield brine can be an indicator of a hydrogeologically stable environment, (2) activities near ground surface may not perturb the stable groundwater environment in the deep brine region, and thus, (3) the deep brine region can be considered as a candidate geologic site for the safe disposal of waste. In addition to brine, other issues associated with long-term waste disposal, such as geological, glacial and seismic events, may need to be considered for the safe storage of spent nuclear fuel in a shield environment.     

井间超声波层析成像及其在岩溶勘察中的应用

井间超声波CT技术采用网格法进行射线追踪，应用逐次线性化迭代与带阻尼最小二乘QR分解相结合的方法进行图象重建. 该技术应用于工程地质勘察领域，能够达到采用相同反演方法的井间地震波CT技术所无法达到的勘察精度要求．工程实践证明，对于井间距超过30 m的情况，超声波CT能分辨大于60 cm的地质异常体. 此方法在工程地质勘察领域具有广阔的应用前景．      

建立和发展岩性地震勘探法

岩性地震勘探法是相对于构造地震勘探法而言的。过去的地震勘探法只是研究地下构造,七十年代以来,逐步发展了亮点法和虚测井等研究地层岩性的方法,包括直接找油气。这些都还只是间接对比岩性的方法,从地震剖面直接提取各种弹性波参数的直接对比方法,还有待于发展。文中评述了间接对比方法的各个步骤的要点及直接对比方法的可能途径。     

复杂地层构造的地震模拟

复杂地层构造的地震模拟,是地震勘探理论的重要问题。本文提出了一个模拟系统:用有限元素表示复杂地层构造,在此基础上应用射线轨迹法计算二维合成记录,从而得到复杂地层构造的理论剖面。这个系统可以用来模拟构造和岩性,为解释人员提供理论依据。 作者研制了适用于任意复杂地层构造计算的通用程序,并用以计算了某测线的一批模型,所得合成剖面由电视加以显示,最后对结果作了简单的分析讨论。     

辽宁省丹东市新城区活动断裂探测及地震危险性评价

在分析区域地震地质资料的基础上,通过遥感信息解译,现场地震地质调查,控制性钻探和综合地球物理探测等项工作,查明了丹东新城区及附近的地震地质条件,认为新城区的地震危险性主要来自鸭绿江断裂,因此,需要考虑抗震设防。     

地震断裂系统的分形和在地震预报中的应用

本文基于地震断层具有自相似结构,从分形概念出发,依据地震资料提出了适合于前震和余震的分形模型,在此基础上讨论了前震和余震的b值特征,提出了地震矩不均衡度的概念。本文还探讨了利用分数维和地震矩不均衡度进行地震预报的可能性。研究发现,正常情况下,地震矩不均衡度呈现低值,大震发生前出现高值异常,震后降低;在空间上,震前2—4年其异常范围大,震前一年异常区逐渐缩小,大震发生在异常集中区内或其边缘。     

地震测深数据文件及其处理系统

本文在我国大规模开展地震测深工作的基础上,提出了地震测深数据文件系统,其目的在于在测深资料解释中引进数字处理技术,增加地震记录可用的信息量,建立地震测深数据库,从而促进我国深部研究工作中资料的相互交流. 在IBM-PC/XT微型计算机上,与本文的数据文件系统相应的数据采集、文件组织、绘制地震记录截面图以及对数字记录进行预处理等方面软件设计已经完成,从而建立了较完善的地震测深数据处理软件包.这个软件包在野外的使用表明,它具有效率高、使用方便的优点.     

ANALYSIS OF DIFFRACTION PATTERN FILTERING OF VARIABLY DENSITY AND VARIABLY AREA RECORDINGS

When a small, transparent replica of a seismic section is illuminated by a homogeneous beam of coherent, monochromatic parallel light a diffraction pattern is created that is representation of the double Fourier spectrum of the recorded seismic waves, i.e. their spectrum in terms of frequency, f, and apparent wave number, k. Masking selected parts of this diffraction pattern causes the spectrum to be filtered: the recomposition of the filtered spectrum then provides a filtered section. The ideal seismic section for this purpose would be a continuous variable density section obtained from recording made at every point of the seismic line. The light transmission coefficient (in terms of light amplitude) at each spot of the replica should be linearly related to the instantaneous seismic signal strength at the spot on the seismic section to which it refers. Unfortunately we cannot make recordings at every point of a seismic line and in our practically realisable recordings we have to be content with sampling in the direction of the location coordinate x. This means that with variable density recordings aliasing will be present and evident in the spectrum obtained in the direction of k; furthermore, the aliased spectrum is also multiplied by a sine function of k because the recording obtained at a given station is not presented on the seismic section as a single line along the time axis, but occupies the entire width of the trace. The diffraction patterns created by variable density recordings of dipping sine waves, including clipped recordings, and of the effect of dip filtering in such sections are discussed. The efficiency of dip rejection is shown to decrease with increasing dip. The diffraction pattern of a variable density recording is found to be characterised by a relationship between point pairs in the pattern. No such simple relationship has been found for the diffraction pattern of a variable area section; the spectra of such VAR sections belong to a very special class, because the amplitude transmission coefficient has only two values, viz. 0 and 1. Consequently, selective masking of the diffraction pattern of a VAR section may give rise to a filtered profile that does not look like a VAR section at all. General statements about the diffraction pattern of VAR sections are hard to give, because the transmission coefficient at a given point in the replica is not proportional to the signal level in the seismic section at the relevant point. In the case of VAR presentation of harmonic waves it was found that, as well as the aliasing effect in the k direction, higher harmonics of the frequency are also introduced. Some synthetic examples are given that show dip filtering to be less effective with VAR than with variable density recordings. Some arguments are advanced in favour of the opinion that high-pass filtering of VAR sections will have less success than low-pass filtering. This is demonstrated by two synthetic examples.     

全波形反演在缝洞型储层速度建模中的应用

速度是地震偏移成像准确与否的关键所在.全波形反演综合利用地震波场运动学和动力学信息,能够得到相比传统速度建模方法更高频的成分.全波形反演的理论比较成熟,但实际应用成功的例子相对较少,特别是对于陆上地震资料.塔里木盆地地震地质条件复杂,为了实现缝洞型储层的准确成像,本文开展了针对目标靶区的全波形反演精细速度建场研究.采用一种时间域分层多尺度全波形反演流程:首先通过层析成像建立初始速度模型;其次利用折射波反演浅层速度模型;最后利用反射波反演中深层速度模型.偏移成像结果表明基于全波形反演的速度建模技术能有效改善火成岩下伏构造的成像精度,显示了全波形反演在常规陆上采集资料的应用潜力.     

叠前时间偏移技术在东营凹陷北部陡坡带砂砾岩体成像中的应用

东营凹陷北部陡坡带是砂砾岩体发育地区.通过多年的研究与勘探实践,相继发现并开发了许多具有陡坡带沉积特色的砂砾岩扇体油藏,砂砾岩扇体油藏已成为胜利油田增储上产的重要勘探目标.为了使砂砾岩体的成像效果更加理想,在该区进行了叠前时间偏移应用研究.在充分认识砂砾岩体地震反射特征的基础上,对东营北带叠前时间偏移的适应性进行了分析,讨论了叠前时间偏移提高陡坡带砂砾岩体成像质量的技术优势,重点介绍了叠前时间偏移处理中的关键技术.处理结果表明,地质认识指导下的叠前时间偏移处理,大大地改善了东营北带砂砾岩体的成像质量,为提高砂砾岩扇体油气藏的勘探精度提供了可靠的基础资料.     

Velocity Structure Beneath Active Faults and Seismic Tectonic Characteristics in the Beijing Area

Velocity structure beneath active faults in the Beijing area has been discussed,based on the digital crustal model of velocity from deep seismic sounding profiles and tomography imaging of P waves. We found that there exists nearly vertical deep faults beneath the Shunyi Depression and the Machiko Depression,which are very likely to be seismogenic faults in the future. In the superficial crust,the north segment of the Shunyi-QianmenLiangxiang fault,the north segment of the Huangzhuang-Gaoliying fault and th...     

分时窗提取地震子波及在合成地震记录中的应用

提出了利用地震和测井资料精确提取井旁地震子波的分时窗提取地震子波方法，将此方法用于合成地震记录的制作，提高了合成地震记录与地震剖面的吻合度和分辨率，文中详细介绍了该方法的具体实现步骤，并给出了模型处理分析和实例分析。     

三维叠前深度偏移速度分析及蒙特卡洛自动层速度拾取

是否能够正确地建立深度域三维速度模型是三维叠前深度偏移成败的关键 .本文根据Deregowski循环 ,利用叠前深度域地震成像对速度模型变化的敏感性 ,采用偏移迭代逐次逼近最佳成像速度 ,研究开发了一套快捷有效的三维叠前深度偏移深度域速度模型建立技术 .借鉴时间域CDP(共深度点 )道集上常规叠加速度分析的策略 ,在深度域CRP(共反射点 )道集上 ,提出剩余慢度平方谱的概念并建立相应的实现技术 .导出深度域中均方根速度与层速度之间的关系 ;按照串级偏移原理确定偏移循环过程中初始速度、剩余速度及修改后速度之间的关系 ;采用蒙特卡洛非线性优化算法实现从剩余慢度平方谱中自动拾取层速度 ,讨论了其地质速度约束条件和蒙特卡洛非线性优化的收敛准则 ,使得所拾取的层速度模型具有合理的地质意义并获得最佳偏移成像效果 .SEG EAGE理论模型数值试算验证了方法的有效性 ,在海拉尔盆地霍多莫尔工区 ,5 8km2 三维资料的速度模型建立并获得满意的三维叠前深度偏移成像 .     

庐枞盆地浅表地壳速度成像与隐伏矿靶区预测

利用反射地震初至波可以精确反演地壳速度结构和构造信息.对庐枞盆地采集的高分辨地震数据,运用初至波层析成像方法,反演得到了该区1200 m以上的浅表地壳速度结构信息.找矿信息总是同地球化学异常联系紧密,对采集于炮孔深处的岩屑、泥砂样品进行了地球化学分析,发现了多处显著的金属元素异常.对五条测线的速度成像结果及其映射的地下岩性结构和炮孔地球化学分析结果进行了相关性对比分析.结果显示,庐枞盆地的沉积岩、火山岩和侵入岩具有不同的速度范围,侵入岩具有高波速特征,地球化学信息高异常往往出现在高速侵入岩体的上方;精细的速度信息蕴含着丰富的浅表地壳结构变化特征,与庐枞矿集区已知的地质、岩体、构造和矿体分布存在着良好的对应关系.速度成像结果可以准确刻画地下隐伏侵入岩体的空间分布形态,结合地球化学分析异常和重磁探测等信息,预测隐伏矿床,提供深部找矿靶区.     

测井技术在地震目标追踪应用中的方法探讨

目前常见测井-地震综合应用技术主要是测井约束地震反演技术,该技术存在两方面的研究精度问题.一是初始模型等因素影响反演的精度和分辨率,二是波阻抗的多解性难免造成反演的多解性.由于存在以上原因,有必要探索其他方式的"井震结合"研究思路,根据测井、地震信息具有地质背景成因一致的特点,提出寻找二者具有的、可相互辨识的共同地质属性的分析方法,研究认为,测井、地震共同地质属性的辨识可表现在两方面,一是测井、地震信息在地质界面上下地质成因相同.深入研究两者可相互辨识的地质属性信息,可提供地质界面的识别与追踪依据;二是对于一个具体的地质事件或地质体,也理应找到具有同一成因基础的、可相互辨识的信息响应依据.利用该方法对羊二庄某区的实践表明,在测井、地震信息成因一致性分析手段的基础上,开展测井-地质共同地质属性的辨识、追踪研究,对复杂油气区有利目标的追踪具有有效性,该研究成果被钻井结果证实.     

INTEGRATED INTERPRETATION, 3D MAP MIGRATION AND VSP MODELLING PROJECT,NORTHERN U.K. SOUTHERN GAS BASIN1

Depth conversion in the northern part of the U.K. Southern Gas Basin is complicated by the presence of Zechstein (Permian) salt swells and diapirs. In addition, the post-Zechstein (post-Permian) section displays large lateral velocity variations. The primary agents which control the velocity of this stratigraphic section are: (1) depth of burial, (2) lithological variation within individual formations, and (3) the effects of subsequent tectonic inversion. An integrated approach which combines well velocity, seismic velocity and seismic interpretation is required for accurate depth estimation. In 1988 Mobil and partners drilled an exploratory well in the northern part of the U.K. Southern Gas Basin. This well was located near the crest of a Zechstein salt diapir. Over 2000 m of Zechstein was encountered in the well. The Permian Rotliegendes objective was penetrated at a depth of over 3700 m. The initial delineation of the objective structure was based on the results of 3D map migration of the seismic time interpretation. Spatially-variant interval velocity functions were used to depth convert through five of the six mapped horizons. Both well and model-based seismic interval velocity analysis information was used to construct these functions. A moving-source well seismic survey was conducted. The survey was run in two critical directions. In conjunction with presurvey modelling, it was possible to confirm immediately the structural configuration as mapped to a distance of 7 km from the well. Post-survey 3D map migration and modelling was employed to further refine the structural interpretation. Although the question of stratigraphic anisotropy was considered in the evaluation of the long offset modelling, no evidence was found in the field data to support a significant effect. Finally, comparisons were made of: curved-ray versus straight-ray migration/modelling, midpoint-depth velocity versus (depth-dependant instantaneous velocity functions, and Hubral- versus Fermat-based map depth migration algorithms. Significant differences in the results were observed for structural dips exceeding 15^o and/or offsets exceeding 6 km. Map depth migration algorithms which employed both curved rays and spatially-variant instantaneous velocity functions were found to best approximate the ‘true’ geological velocity field in the study area.