非规则三角网格有限差分法地震正演模拟

针对有限差分方法在地震正演模拟中的广泛应用,提出了1种基于三角形网格的地震正演模拟差分算法。相对于矩形网格差分法而言,该方法可以更为精确地描述弯曲的速度界面,可直接用于起伏地表模型的正演模拟计算而无需对模型做特殊处理。由于空间步长可随速度的变化而改变,加之在递推计算中采用了显式差分算法,该方法可视为对传统的规则网格有限差分法的有效改进。数值模拟的结果表明该方法是稳定、精确的且具有较高的计算效率。     

地震正演模拟技术在准噶尔盆地西北缘七区地震解释中的应用

地震正演模拟技术是用数学或物理的方法研究地震波在某种具体的简化模型中传播的特点,用来模拟真实地质结构条件下的地震波场.随着油气勘探的难度越来越大,该技术得到了越来越广泛的应用.以准噶尔盆地西北缘七区三维地震解释中的应用为例,正演模拟技术在地质模式的识别、断裂的精细刻画、速度异常的分析、小断裂的识别等方面都有较好的效果.     

基于多重网格的电火花枪阵震源波场正演模拟

基于电火花震源的超高分辨率地震勘探在我国的海洋区域调查以及工程地球物理勘探中发挥着重要的作用。实际勘探时,电火花震源间隔通常为厘米级,为电火花枪阵波场的数值模拟带来困难。针对这一问题,在常规变网格有限差分算法的基础上,在电火花枪阵波场的数值模拟中提出了多重网格策略,实现了基于多重网格算法的电火花枪阵波场高精度三维数值模拟。数值模拟实验表明,相比于常规的交错网格有限差分数值模拟,基于多重网格有限差分数值模拟算法能够显著提高计算效率,降低内存消耗,同时还可有效地压制常规变网格算法中的虚假反射现象,实现电火花枪阵震源波场的高精度数值模拟。     

基于纵横波保幅分离的粘弹介质弹性波正演模拟

地下介质往往表现为粘弹特性,研究基于粘弹假设的数值模拟方法对于正确认识地震波的传播规律和提高地震勘探精度具有重要意义。波动方程正演过程中纵横波的保幅解耦是研究准确的粘弹介质中地震波传播机理的前提,基于散度和旋度算子的纵横波解耦方法会使波场的相位和振幅产生畸变,且解耦后的波场在极性反转位置上无法与分离前混合波场各分量对应。在散度和旋度算子上再做一次梯度和旋度运算的波场分离方法虽然能够克服上述缺陷,但存在保幅性差等问题。本文从粘弹介质中的一阶速度-应力方程出发,推导了矢量纵横波分离的波数域表达式,结合有限差分思路给出了其在空间域的求解方法。本文方法利用纵横波的传播速度对现有的矢量波场分离方法进行振幅校正,并将校正结果分别作为纵波与横波对时间的二阶偏导,实现了粘弹介质中的纵横波分离。模型试算结果表明,本文方法能够克服现有方法的缺陷,获得更具保幅性的波场分离结果。     

一种斜井VSP的正演模型

在垂直入射斜井ＶＳＰ的条件下，假定地下为倾斜层状平面结构，在各层中地震射线为直线，界面的反射和透射服从Ｓｎｅｌｌ定律时，本文提出一种射线追踪法正演模型。射线追踪是采用迭代法实现的。用理论模型试算可以得到较为满意的结果。     

复杂地形结构下瞬变电磁三维正演模拟

如何确定海底多金属硫化物的空间分布并对其进行资源量评价至关重要。瞬变电磁法是陆地金属矿床勘探的重要手段,但海底的崎岖地形、硫化物的复杂结构以及不稳定的近底观测条件为瞬变电磁在海洋中的应用提出了挑战。为验证瞬变电磁方法对热液区硫化物的勘探应用潜力,本文采用了有限元方法,结合大西洋洋中脊TAG热液区实测地形数据和硫化物深部剖面进行三维正演模拟,通过对比仪器在不同位置、姿态、离底高度的正演模拟结果发现:重叠回线探测装置在离底高度小于60 m时,可以有效探测到TAG丘体的高电导率异常。海底复杂地形以及仪器拖曳方式都会对二次场早期响应产生干扰,同时仪器姿态变化也会改变所探测到的响应,这表明,应结合研究区域的海底地形、仪器海底定位以及姿态数据才能更好的对实测瞬变电磁数据进行合理解释。     

琼东南盆地南部早中新世碳酸盐岩台地的地震解释

The southern uplift of the Qiongdongnan Basin is a deepwater area in which no wells have beens drilled. The Miocene-Quaternary strata in the Xisha Islands, which are located 40–100 km to the south, are composed of carbonate reef formations. Paleotectonic and paleogeographic analyses of the basin suggest that the southern uplift experienced favorable geological conditions for the development of carbonate reefs during the Miocene.The high-impedance carbonates have high amplitudes and low frequencies on seismic profiles. The reefs are distributed on paleotectonic highs and are thicker than the contemporaneous formations. A forward model of the variation in carbonate thickness based on lithological and velocity information from wells in nearby regions can simulate the seismic response of carbonates with different thicknesses. We identified several important controlling points for determining the thickness of carbonates from seismic profiles, including the pinchout point,the λ/4 thickness point, and the λ/2 thickness point. We depict a carbonate thickness map in the deepwater area of the southern Qiongdongnan Basin based on this model. The carbonate thickness map, the paleotectonic and paleogeographic background, and the seismic response characteristics of reefs suggest that the carbonates that developed on the southern uplift of the Qiongdongnan Basin during the Miocene were mainly an isolated carbonate platform peninsula and ramp deposits. It consisted of gentle ramp platform, steep slope platform,platform depression, gravity flow, and reef bank facies.     

苏北朱家墩气田成藏特征对南黄海南部盆地勘探的意义

南黄海南部盆地是下扬子地台沿北东方向向海域的自然延伸,盆地的区域构造及油气地质特征与陆上苏北盆地相似。苏北盆地朱家墩气田的成功开发为南黄海南部盆地的勘探提供了借鉴,在分析朱家墩气田成藏特征的基础上,将南黄海南部盆地的成藏条件与朱家墩气田成藏条件进行对比,发现两者具有很强的可比性,得出南黄海南部盆地具有良好的勘探前景,其勘探方向为寻找"古生新储"的气藏,其发育部位为邻近大断裂一侧的深凹带。     

复杂地形结构下瞬变电磁三维正演模拟以大西洋洋中脊TAG热液区为例

如何确定海底多金属硫化物的空间分布并对其进行资源量评价至关重要。瞬变电磁法是陆地金属矿床勘探的重要手段,但海底的崎岖地形、硫化物的复杂结构以及不稳定的近底观测条件为瞬变电磁在海洋中的应用提出了挑战。为验证瞬变电磁方法对热液区硫化物的勘探应用潜力,本文采用了有限元方法,结合大西洋洋中脊TAG热液区实测地形数据和硫化物深部剖面进行三维正演模拟,通过对比仪器在不同位置、姿态、离底高度的正演模拟结果发现：重叠回线探测装置在离底高度小于60 m时,可以有效探测到TAG丘体的高电导率异常。海底复杂地形以及仪器拖曳方式都会对二次场早期响应产生干扰,同时仪器姿态变化也会改变所探测到的响应,这表明,应结合研究区域的海底地形、仪器海底定位以及姿态数据才能更好的对实测瞬变电磁数据进行合理解释。     

珠江口盆地白云深水区含气储层AVO模板的建立及应用

白云深水区多口油气钻探揭示砂质储层含气和含水都有可能表现为Ⅲ类或Ⅳ类的AVO及亮点异常,单纯依靠亮点+Ⅲ类AVO异常进行烃类气体预测具有多解性,是钻前储层的流体预测面临的新挑战和难点。本文利用珠江口盆地白云深水区测井数据,根据储层特征与含气性差异性优选了多口测井数据齐全井,利用测井数据对不同井位储层段的AVO特征进行分析。通过Aki-Richards公式计算了不同井的截距和梯度属性,建立了白云深水区P和G属性图版。对比分析发现不同含气储层,随着岩性、含气性差异及其岩性组合不同,呈现不同AVO异常特征,岩心分析表明该异常不仅与储层含气性有关,也与盖层岩性密切相关,尤其是当盖层或者储层含有灰岩时,对AVO异常影响较大。利用M矿区目标储层的叠前反演弹性参数,再结合本文建立的P和G属性图版和流体识别因子方法,对目标储层含气性进行了预测,发现了该储层在横向上含气性不同,钻探结果证实预测结果有效与可靠,表明该AVO模板具有较好的实用性。     

天然气水合物赋存地层的地震波场模拟及对比

通过分析天然气水合物在海洋中的6种主要赋存状态类型,总结了每种赋存状态之间的相互转化关系及其物性参数计算方法,并应用到地震波场的正演模拟中.对比研究了声波模型、弹性波模型和双相介质模型对各种水合物赋存地层的响应特征,结果表明:1)当地层中存在孔隙充填型水合物且下伏地层不含游离气时,双相介质模拟的含水合物层底界表现负极性...     

Application of AVO attribute inversion technology to gas hydrate identification in the Shenhu Area,South China Sea

AVO (Amplitude Versus Offset) is a seismic exploration technology applied to recognize lithology and detect oil and gas through analyzing the feature of amplitude variation versus offset. Gas hydrate and free gas can cause obvious AVO anomaly. To find geophysical evidence of gas hydrate and free gas in Shenhu Area, South China Sea, AVO attribute inversion method is applied. By using the method, the multiple seismic attribute profiles and AVO intercept versus gradient (I-G) cross plot are obtained. Bottom-simulating reflector (BSR) is observed beneath the seafloor, and the AVO abnormal responses reveal various seismic indicators of gas hydrate and free gas. The final AVO analysis results indicate the existence of gas hydrate and free gas in the upper and lower layers of BSR in the study area.     

计算域的选取对风暴潮数值模拟的影响

作者以 72 0 3号台风过程中引起的黄海沿岸水文测站的风暴潮过程为例 ,选取不同计算域 ,模拟计算的同一测站的风暴潮增水值有很大差异。只有选取整个黄渤海海域 ,才能得到较精确的风暴潮增水结果。从而揭示在封闭或半封闭海域中由风暴过程激发的区域自由振荡 ,是风暴潮增水中不容忽视的量。说明该海域中的风暴潮过程是海域整体效应的响应 ,因此认为数值模拟中计算域应选取整个封闭或半封闭海域     

Numerical Simulation of Wave Field at the Kemena River Mouth in Kalimantan Island

-By use of the parabolic equation of numerical simulation of wave which is suitable forlarge-angle propagation and Crank-Nicolson differential method,the wave field at the Kemema Rivermouth has been studied for analysis of sediment movement in the area.In order to reflect wave energy lossaccurately,the Bretchneider-Reid formula is quoted and the friction coefficient in the formula is discussedin this paper.The calculation results indicate that the wave becomes a little damped at the mouth ofKemena River influenced by the topography and bottom friction,and the wave at the east beach is higherthan that at the west beach,because the east beach extends out.     

吕宋海峡流场季节变化的数值模拟

基于海洋模式HYCOM(Hybrid Coordinate Ocean Model),利用大小区嵌套技术,分别对全球海洋和西北太平洋进行了网格嵌套数值模拟,研究了吕宋海峡海域环流场的季节性变化。考虑全球海洋环流影响的西北太平洋模式,成功地刻画了黑潮的流结构及季节变化。吕宋海峡海洋环流流场在不同深度处差异较大,存在着明显的季节变化。黑潮入侵南海主要发生在500m深度以上,冬季最明显,夏秋两季不明显。在500m层常年存在一支南海暖流流入西北太平洋,在800m层南海暖流消失。一年四季黑潮主要通过吕宋海峡的南部和中部进入南海。1 000m层流场表明,黑潮主要通过吕宋海峡的中部入侵南海。在800~1 000m处主要是黑潮水流入南海。     

水声通信半物理仿真平台设计与应用

信道自身的复杂性和现场级试验的高投入、高风险限制了水声通信技术的发展。针对高效率水声通信技术研究的需求,提出1种水声通信半物理仿真平台水声通信半物理仿真平台(A Hardware-in-Loop Simulation Platform for Underwater Acoustic Communication)。水声通信半物理仿真平台以可配置水声Modem为核心,包括现场测试采集系统、算法仿真评估系统、算法实现装载系统、水声波形存储播放系统。基于建立的水声通信半物理仿真平台,进行水声直接序列扩频通信系统的半物理仿真研究,结果表明水声通信半物理仿真平台可有效提高水声通信技术的研究效率。     

Velocity and AVO analysis for the investigation of gas hydrate along a profile in the western continental margin of India

The occurrence of gas hydrate has been inferred from the presence of Bottom-Simulating Reflectors (BSRs) along the western continental margin of India. In this paper, we assess the spatial and vertical distribution of gas hydrates by analyzing the interval velocities and Amplitude Versus Offset (AVO) responses obtained from multi-channel seismics (MCSs). The hydrate cements the grains of the host sediment, thereby increasing its velocity, whereas the free gas below the base of hydrate stability zone decreases the interval velocity. Conventionally, velocities are obtained from the semblance analysis on the Common Mid-Point (CMP) gathers. Here, we used wave-equation datuming to remove the effect of the water column before the velocity analysis. We show that the interval velocities obtained in this fashion are more stable than those computed from the conventional semblance analysis. The initial velocity model thus obtained is updated using the tomographic velocity analysis to account for lateral heterogeneity. The resultant interval velocity model shows large lateral velocity variations in the hydrate layer and some low velocity zones associated with free gas at the location of structural traps. The reflection from the base of the gas layer is also visible in the stacked seismic data. Vertical variation in hydrate distribution is assessed by analyzing the AVO response at selected locations. AVO analysis is carried out after applying true amplitude processing. The average amplitudes of BSRs are almost constant with offset, suggesting a fluid expulsion model for hydrate formation. In such a model, the hydrate concentrations are gradational with maxima occurring at the base of hydrate stability zone.     

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model,a quasi-3D numerical model is establishedfor coastal regions of shallow water.The characteristics of this model are that the velocity profiles can be ob-tained at the same time when the equations of the value of difference between the horizontal current velocityand its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent withthe results of the 2D model.The circulating flow in the rectangular area induced by wind is simulated and ap-plied to the tidal flow field of the radial sandbanks in the South Yellow Sea.The computational results fromthis quasi-3D model are in good agreement with analytical results and observed data.The solution of the finitedifference equations has been found to be stable,and the model is simple,effective and practical.     

Field Test and Numerical Simulation of Tidal Bore Pressures on Sheet-Pile Groin in Qiantang River

A field test of bore pressures on a sheet-pile groin in the Qiantang River was carried out to investigate the characteristics of tidal bore pressures in which temporal and spatial distributions of bore pressures on the sheet-pile groin were obtained through the test. The standard turbulence model was used to simulate the propagation of the tidal bore over a sheet-pile groin with the free surface tracked by the volume of fluid (VOF) method. With the simulated bore pressure over the sheet-pile groin, a finite difference method was used to simulate the dynamic response of the sheet-pile groin subjected to bore pressures. Good agreement between the field test and the simulated results validated the present numerical method. It was found that the bore pressures on the sheet-pile groin varied linearly in the vertical direction and obtained maximum value at the still water level. Bore pressures below the still water level decreased slightly with the depth and were characterized by a trapezoidal distribution.