任意空间取向TI介质中速度随方位变化特征

本文基于任意空间取向TI介质中坐标变换的方法,扩展研究了任意强弱、具有任意空间取向对称轴的TI介质中体波相速度和群速度的方位变化;通过模型数值计算,获得了相速度与群速度方位变化图案,并比较了二者偏差.研究发现,各向异性越强群速度与相速度的偏差越大;并且,TI对称轴的空间取向和测线方位影响速度方位变化.相对于TI对称轴,速度方位变化图案不变;但是,随着TI对称轴空间取向和测线方位的改变,其速度方位变化呈现一定的规律性.研究结果可以精确地预测任意空间取向TI介质中速度方位变化,有利于地震各向异性数据处理与解释.     

半空间饱和介质内圆形洞室对平面P1波的散射

考虑土颗粒和孔隙流体的压缩性以及孔隙流体与土骨架之间的粘性耦合作用,采用修正的Biot模型,假定半空间表面不透水,得到了平面P1波(快压缩波)在半空间表面的反射P1波、P2波(慢压缩波)和SV波(剪切波)的幅值.采用大圆弧假定将半空间内圆形洞室的散射问题转化为大圆弧和圆形洞室的多重散射问题,运用波函数展开法将入射波、反射波以及半空间表面和洞室的散射波的势函数展开成Fourier-Bessel函数的无穷级数形式,由Graf加法定理得到同一坐标系下的势函数的表达式,根据半空间表面和洞室完全自由的边界条件得到了待定复系数的理论解.通过数值计算,着重分析了平面P1波垂直向上入射时无量纲入射频率和洞室埋藏深度等对洞室的动应力集中因子和半空间表面的归一化水平和竖向位移的影响.     

Exact traveltime computation in multi-layered transversely isotropic media with vertical symmetry axis

An approach to calculate the accurate ray paths and traveltimes in multi-layered VTI media (transversely isotropic media with a vertical symmetry axis) is proposed. The expressions of phase velocity, group velocity and Snell’s law used for computation are all explicit and exact. The calculation of ray paths and traveltimes for a given ele-mentary wave is equivalent to that of a transmission problem which is much easier to be treated with the formulae proposed. In the section of numerical examples, the proce...     

Imaging and image processing in porous media research

Three-dimensional imaging and image processing has become an important part for investigations of fluid distribution and flow in porous media. We describe two methods of computed tomography with different characteristics, namely X-ray- and neutron-based. We give an overview of image processing sequences and their methods. We investigated image enhancement with a focus on filters using partial differential equations, classification and structure identification that we used to prepare our images for quantitative evaluations. These methods are demonstrated on a partially saturated sand sample. Finally, we show an application with soil aggregates where investigations using synchrotron X-rays and thermal neutrons have led to new insights and refined fluid distribution and flow models.     

Drainage in heterogeneous sand columns with different geometric structures

This paper discusses multi-step drainage experiments in two heterogeneously packed sand columns (10 × 10 × 20 cm³). Different packing structures were generated using two different sand types. One purpose of the study was to test the influence of packing structures on the movement of water. The second purpose was to assess the quality of predictions for the outflow curves in both columns made with an upscaled model. The heterogeneous structures of the columns can be considered as two opposing extremes. The first column was packed with a random arrangement of two sand types that is not stochastically homogeneous and where a cluster running through the column exists for both materials. The second column was packed with a periodic pattern of coarse-sand inclusions in a fine-sand background and has a clearly defined unit cell. The depth-averaged (2D) spatial distribution of the water content in the columns was monitored during the whole multi-step outflow experiment using neutron radiography. The 3D water content was measured at the steady states by neutron tomography. The experimental results are compared with the model predictions of an upscaled model derived with the homogenization theory. The parameters for the upscaled model are calculated from the hydraulic parameters of the two sand types. These hydraulic parameters were first identified in independent measurements on samples of the two individual sand types, separately. Additionally, the hydraulic parameters of both sands were identified by fitting a numerical model to the measured outflow curves. The different column structures showed a significant effect on water retention and the effective retention function, as water was trapped in the coarse-sand inclusions of the periodic structure. We included this trapping effect in the effective retention function of the upscaled model with an apparent air entry pressure. Contrary to the retention, the different packing structures had no large effect on the dynamic behavior of the outflow. The effective conductivity of the columns is therefore not significantly influenced by the structure. The upscaled models predicted the movement of the averaged water content in the two columns well. This confirms the applicability of upscaled models even if the underlying requirements are not strictly met.     

Accumulation characteristics and the main controlling factors of Lunnan multilayer oil province,Tarim Basin,China

Lunnan region is a large-scale paleohigh with many coexisting oil and gas bearing series.At present, about 2 billions tons of proved,probable and possible oil and gas reverses have been proved there.Eight oil and gas bearing series have been found in the Ordovician,Carboniferous,Triassic and Jurassic of Lunnan region,they all bear the characteristics of large-scale multilayer oil-gas province.Ordovician is the main reservoir series where over 0.8 billion tons of oil geologic reserves were discovered,and a super large-scale marine carbonate oil and gas field has formed.Reservoir space of the carbonate reservoirs is mainly composed of dissolved hole,dissolved pore and fracture in Lunnan paleo-burial hill.Generally, dissolved holes are widely distributed among them.Reservoir developments are mainly controlled by karstification and tectonic disruption.Due to the similar geochemical characters,the Ordovician,Carboniferous,Triassic and Jurassic oil and gas reservoirs present the same oil source rock of Mid-Upper Ordovician,the latter except Ordovician are mostly of secondary oil and gas reservoirs migrated vertically by faults during the process of multiple phase tectonic movement,adjustment and reconstruction. Lunnan composite oil and gas accumulation region is situated in the vicinity of large-scale hydrocarbon generation depressions in three directions,ample oil and gas from hydrocarbon generation depressions supplied the adjacent oil and gas reservoirs once.Hereby,the succeed paleohigh is the long-term hydrocarbon accumulation region,which is favor for the formations of high quality reservoirs,fault systems and huge-scale composite oil and gas accumulation.     

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格（LG）进行TTI介质的正演模拟研究,避免了Virieux标准交错网格（SSG）算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性：（1）结合最小二乘思想推导出新的频散改进差分系数（DIC）,该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;（2）将分裂的多轴完全匹配层（M-PML）吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.     

改进的TTI介质纯P波方程正演模拟与逆时偏移

逆时偏移作为一种高精度偏移方法已成为复杂构造成像的重要技术,描述纵波独立传播的延拓方程是各向异性介质逆时偏移的一个关键问题.在对VTI介质几个经典相速度近似公式回顾的基础上,针对常用于描述纯P波的Harlan近似公式在各向异性参数ε较大情况下近似精度较低的问题,本文对Harlan公式中的非椭圆项进行了修正,在非椭圆项前添加了一个与各向异性参数ε有关的修正系数,得到了三种改进型Harlan公式,并以近似精度最高的改进式为基础,推导了TTI介质纯P波方程.针对该伪微分方程,本文利用伪谱法和有限差分法联合实现波场延拓,对于常密度二阶方程,基于中心网格实现;对于一阶应力-速度方程则基于旋转交错网格实现.通过数值试验分析了TTI介质纯P波一阶应力-速度方程的近似精度,并以一阶纯P波方程为基础进行了TTI介质逆时偏移数值模拟试验.结果表明,本文给出的方法能够较准确地描述TTI介质纯P波波场特征,可以应用至各向异性介质逆时偏移.     

利用2.5维数值模式匹配算法高效高精度建立垂直井眼中多分量阵列感应井眼校正库

本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器(TDIT)的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响.