GPS单频精密单点定位方法与实践

讨论了单频精密单点定位的解算方法与思路,编写了单频精密动态单点定位的程序,实现了事后动态单频精密单点定位.     

制图综合中移位算法概述分析

本文首先对移位算法的发展史进行简要介绍,然后,在其基础上,根据各移位算法的思想特征,对各移位算法进行归纳,将其划分为七类,并对每一类算法的核心思想进行详细论述,最后,对各类移位算法进行评价,分析其优缺点,总结各类算法中的可借鉴思想,为更好地实现移位奠定基础。     

垦71区块井间地震资料处理研究与应用

由于避开了地表低降速带和风化层的影响,井间地震能够采集到高频地震信号,经过资料处理可以得到较高分辨率的速度和上行反射波成像剖面。2004～2005年,胜利油田在垦71 区块实施了8对井间地震资料采集。针对原始资料特点,研究了以层析成像、管波衰减、VSP-CDP技术和有限角度叠加成像等关键处理技术,较好地完成了该区块的井间地震资料处理。初步分析解释表明,井间地震可以识别低序级断层、小幅度构造和薄储层。     

小波变换在遥感图像上矿区空间定位方面的分析与应用——以鄂东南铁矿矿区为例

随着社会的发展,人们对自然地物的影响越来越显著,大量地改变了地表的状况,在遥感图像上则表现为灰度值变化剧烈,其中的纹理信息发生改变.在频率域上说,人工改变的区域部分的高频成分比其他区域部分要丰富很多,小波变换能够有效地将其中的高频信息提取出来,用假彩色分割图将其中的异常值明显地表示出来.湖北大冶的铁山矿区分布着六大露天矿区,地表开挖范围广,影响程度大,以该地区为例,分析小波变换在遥感图像上矿区空间定位方面的应用.使用本文介绍的方法处理效果比较明显.矿区位置的确定还需要参考断层分布信息,这意味着使用小波变换进行处理的方法具有一定的局限性.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

Multi-Scale Texture Modeling

Patterns of crystallographic preferred orientation are referred to as texture. The specific subject of texture analysis is the experimental determination and interpretation of the statistical distribution of orientations of crystals within a specimen of polycrystalline material, which could be metals or rocks. The objective is to relate an observed pattern of preferred orientation to its generating processes and vice versa. In geosciences, texture of minerals in rocks is used to infer constraints on their tectono-metamorphic history. Since most physical properties of crystals, such as elastic moduli, the coefficients of thermal expansion, or chemical resistance to etching depends on crystal symmetry and orientation, the presence of texture imparts directional properties to the polycrystalline material. A major issue of mathematical texture analysis is the resolution of the inverse problem to determine a reasonable orientation density function on SO(3) from measured pole intensities on , which relates to the inverse of the totally geodesic Radon transform. This communication introduces a wavelet approach into mathematical texture analysis. Wavelets on the two-dimensional sphere and on the rotational group SO(3) are discussed, and an algorithms for a wavelet decomposition on both domains following the ideas of Ta-Hsin Li is given. The relationship of these wavelets on both domains with respect to the totally geodesic Radon transform is investigated. In particular, it is shown that the Radon transform of these wavelets on SO(3) are again wavelets on . A novel algorithm for the inversion of experimental pole intensities to an orientation density function based on this relationship is developed.     

Latitudinal variation of measured O3 photolysis frequencies J(O1D) and primary OH production rates over the Atlantic Ocean between 50° N and 30° S

The latitudinal variation of the photolysis frequency of ozone to O(¹D) atoms, J(O¹D), was measured using a filter radiometer during the cruise ANT VII/1 of the research vessel Polarstern in September/October 1988. The J(O¹D) noon values exhibited a maximum of 3.6×10^-5 s^-1 (2 sr) at the equator and decreased strongly towards higher latitudes. J(O¹D) reached highest values for clean marine background air with low aerosol load and almost cloudless sky. The J(O¹D) data, measured under these conditions and a temperature of 295 K, can be expressed by: % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOsaiaacI% cacaqGpbWaaWbaaSqabeaaiiaacqWF8baFaaGccaqGebGaaeykaiaa% bccacqWF9aqpcaqGGaGaaeyzaiaabIhacaqGWbGaaeiiaiaabUhacq% GHsislcaaI4aGaaiOlaiaaicdacaaIYaGaeyOeI0IaaGioaiaac6ca% caaI4aGaaiiEaiaaigdacaaIWaWaaWbaaSqabeaacqGHsislcaaIZa% aaaOGaaeiiaiaabIhacaqGGaGaam4uaiabgUcaRiaaiodacaGGUaGa% aGinaiaacIhacaaIXaGaaGimamaaCaaaleqabaGaeyOeI0IaaGOnaa% aakiaadofadaahaaWcbeqaaiaaikdaaaGccaGG9bGaaeikaiaaboha% daahaaWcbeqaaiabgkHiTiaaigdaaaGccaGGPaaaaa!5EE9!\[J({\text{O}}^| {\text{D) }} = {\text{ exp \{ }} - 8.02 - 8.8x10^{ - 3} {\text{ x }}S + 3.4x10^{ - 6} S^2 \} {\text{(s}}^{ - 1} )\] where S represents the product of the overhead ozone column (DU) and the secant of the solar zenith angle. The meridional profile of the primary OH radical production rate P(OH) was calculated from the J(O¹D) measurements and simultaneously recorded O₃ and H₂O mixing ratios. While the latitudinal distribution of J(O¹D) and water vapour was nearly symmetric to the equator, high tropospheric ozone levels up to 40 ppb were observed in the Southern Hemisphere, SH, resulting in higher P(OH) in the SH.     

弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.