基于倾角中值滤波法的浅层反射地震叠前信噪分离技术

浅层反射地震资料中面波一般比较发育,对有效信号形成干扰,现有压制面波等相千干扰的处理手段主要适用于中深部、深部地震勘探,但并不适用于浅层反射地震资料的处理。为此提出了倾角中值滤波法。该方法是基于叠前地震单炮道集上反射波与强线性干扰在￡吨域视速度上的差异,在多组视倾角范围内,求取一个最佳视倾角,将最佳视倾角的振幅序列的中值视为相干干扰,利用“减去法”,达到信噪分离。通过模型与实际资料的处理,汪明该方法能够有效分离强线性干扰,并保留了强线性干扰的特性,提高了资料的信噪比,保真度高,是浅层高分辨率地震勘探的一种理想的叠前信噪分离技术。     

准噶尔盆地东部石炭系火山岩气藏石油地质特征及有利勘探方向分析

准噶尔盆地火山岩油气勘探取得重大进展。笔者经综合研究,提出影响该区天然气藏形成的石油地质主要特征是：前陆坳陷控制了烃源岩区域展布,后碰撞火山岩是主要的储层,古生代晚期的构造优化了储层物性并形成圈闭。松喀尔苏组为前陆坳陷磨拉石建造的近端粗碎屑岩沉积,是克拉美丽造山事件后快速沉积产物,作为主要烃源岩的下石炭统滴水泉组,是与松喀尔苏组同期前陆坳陷不同位置的沉积层系,滴水泉组是前陆坳陷中央的细粒沉积层系。后碰撞火山岩覆盖在前陆坳陷沉积层系之上,分布广泛,分布范围不受古生代造山带和现今构造单元控制。二叠纪发育断陷构造,形成凹隆相间的构造格局,在构造高部位的隆起区石炭系火山岩风化壳和构造裂缝,优化改善了火山岩储层的储集性能,形成有利油气运移通道和储集空间。克拉美丽山南缘的前陆坳陷中央为有利烃源岩分布区,该范围内的构造高部位的火山岩发育区,将是准噶尔盆地东部下一步火山岩油气藏勘探的有利区域。     

A fast integral equation solver for 3D induction well logging in formations with large conductivity contrasts

Simulation of induction logging responses in formations with large conductivity contrasts is an important but challenging problem due to the singularity of a linear system caused by large contrasts. Also, three‐dimensional (3D) analysis of complex geophysical structures usually encounters high computational demands. In this paper, a pre‐corrected fast Fourier transform (pFFT)‐accelerated integral equation method is applied to overcome these difficulties. In the approach, the entire formation is included in the solution domain. The volume integral equation is set up in the region based on the fact that the total field is the summation of the excitation field and the secondary field. The emitted field by the transmitter coil (treated as a magnetic dipole) is regarded as the excitation of the system. Then the method of moments (MoM) is used to solve the integral equation. To reduce the high computational requirements of the MoM, the pFFT method is used to speed up the solution of the matrix equation and reduce the memory requirement as well. The resultant method is capable of computing induction logging problems involving large and complex formations. For problems with high conductivity contrasts, the solution of the matrix equation usually converges very slow or even fails to converge due to the large condition number of the coefficient matrix. To overcome this difficulty, an incomplete LU pre‐conditioner is used to significantly speed up the convergence of the matrix equation, thus further reducing the computation time. Numerical results show that the present method is efficient and flexible for 3D simulation of induction logging and is specifically superior for problems with high conductivity contrasts.     

Comparison of infiltration models with regard to design of rectangular infiltration trenches

ABSTRACT

Guidelines for the design of infiltration trenches of rectangular cross-section are compared. Four concepts of infiltration into a native soil from linear soakaways of rectangular cross-section are discussed. The results of calculations using analytical equations are compared with numerical two-dimensional (2D) simulations using HYDRUS computer code and field measurements. Satisfactory agreement is achieved for Kozeny’s concept of infiltration and the free-surface approach, applied to a trench in a highly saturated sandy loam. Theoretical solutions, neglecting matric potential, suggest that the mean infiltration rate along the wetted perimeter of a rectangular trench varies in the range 1.0 to 1.5 of the value of saturated hydraulic conductivity.     

3D electrical resistivity inversion using prior spatial shape constraints

To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical data regarding the extent of a geological anomaly are important prior information. We propose the use of shape constraints in 3D electrical resistivity inversion, Three weighted orthogonal vectors （a normal and two tangent vectors） were used to control the resistivity differences at the boundaries of the anomaly. The spatial shape of the anomaly and the constraints on the boundaries of the anomaly are thus established. We incorporated the spatial shape constraints in the objective function of the 3D resistivity inversion and constructed the 3D resistivity inversion equation with spatial shape constraints. Subsequently, we used numerical modeling based on prior spatial shape data to constrain the direction vectors and weights of the 3D resistivity inversion. We established a reasonable range between the direction vectors and weights, and verified the feasibility and effectiveness of using spatial shape prior constraints in reducing excessive structures and the number of solutions. We applied the prior spatially shape-constrained inversion method to locate the aquifer at the Guangzhou subway. The spatial shape constraints were taken from ground penetrating radar data. The inversion results for the location and shape of the aquifer agree well with drilling data, and the number of inversion solutions is significantly reduced.     

东乌旗中铁陨石硅酸盐地球化学特征及对成因的启示

东乌旗中铁陨石为1995年在内蒙发现的、为数不多的降落型中铁陨石,由金属—硫化物和硅酸盐相紧密混合而成。硅酸盐相由角砾部分（主要为角砾状橄榄石和辉长岩岩屑）和非角砾部分组成。岩相学特征、主量元素成分及REE地球化学特征表明角砾状橄榄石、非角砾橄榄石和辉长质硅酸盐（辉长岩岩屑和除橄榄石外的非角砾部分）三者可能分别结晶自不同的岩浆源区。进一步的分析表明Mg含量非常高的角砾状橄榄石可能来自一个分异小行星的橄榄岩质地幔,非角砾橄榄石和辉长质硅酸盐可能和HED陨石具有相同的母体Vesta小行星。冲击碰撞导致两个小行星的物质发生了混合。冲击作用产生的高温（约1200℃—1450℃）使得被撞击小行星表壳玄武质和辉长质岩石发生混合重熔,熔体结晶后形成中铁陨石的辉长质硅酸盐。该温度下非角砾橄榄石和角砾状橄榄石均未发生重熔。该模型较好地解释了角砾状橄榄石、非角砾橄榄石和辉长质硅酸盐三者之间的不平衡现象。