GPS技术在山区石油地震勘探测量中的应用

针对GPS技术在山区石油地震勘探测量中的工作方式,结合2007年度南盘江一桂中地区二维地震大剖面资料采集应用GPS的工程实例,介绍GPS技术应用的新进展。通过GPS静态测量和应用RTK技术放样物理点,从测前准备、数据处理、精度分析等方面,阐述提高作业精度、有效进行质量控制的方法,旨在最大限度地发挥GPS技术在山区施工中的作用。     

Research Advances and Exploration Significance of Large-area Accumulation of Low and Medium Abundance Lithologic Reservoirs

In recent years, a series of large low and medium abundance oil and gas fields are discovered through exploration activities onshore China, which are commonly characterized by low porosity-permeability reservoirs, low oil/gas column height, multiple thin hydrocarbon layers, and distribution in overlapping and connection, and so on. The advantageous conditions for large-area accumulation of low-medium abundance hydrocarbon reservoirs include： （1） large （fan） delta sandbodies are developed in the hinterland of large flow-uncontrolled lake basins and they are alternated with source rocks extensively in a structure like ＂sandwiches＂; （2） effective hydrocarbon source kitchens are extensively distributed, offering maximum contact chances with various sandbodies and hydrocarbon source rocks; （3） oil and gas columns are low in height, hydrocarbon layers are mainly of normal-low pressure, and requirements for seal rock are low; （4） reservoirs have strong inheterogeneity and gas reservoirs are badly connected; （5） the hydrocarbon desorption and expulsion under uplifting and unloading environments cause widely distributed hydrocarbon source rocks of coal measures to form large-area reservoirs; （6） deep basin areas and synclinal areas possess reservoir-forming dynamics. The areas with great exploration potential include the Paleozoic and Mesozoic in the Ordos Basin, the Xujiahe Formation in Dachuanzhong in the Sichuan basin, deep basin areas in the Songliao basin etc. The core techniques of improving exploration efficiency consist of the sweetspot prediction technique that focuses on fine characterization of reservoirs, the hydrocarbon layer protecting and high-speed drilling technique, and the rework technique for enhancing productivity.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

Laser probe Ar/Ar dating: History and development from a technical perspective

The ⁴⁰Ar/³⁹Ar method using a laser probe opened the door to microscale measurements and diffusion profiles frozen in samples. In the first decade since the initial application of a laser for ⁴⁰Ar/³⁹Ar dating in 1973, practical applications have been few. This is due not only to the fact that the laser and vacuum technologies were immature but that mass spectrometry was also in its infancy. In those days, the sensitivity of a mass spectrometer was generally insufficient to measure the small amount of argon degassed from a geological sample by a laser. These problems have subsequently been solved by new technologies. To understand their current status, a brief history of their development is outlined. This outline focuses on the required detection limit in micro scale measurement, practical approaches for accurate measurement are explained through examples in our laboratory specifically relating to the technical aspects of ⁴⁰Ar/³⁹Ar dating.     

U-Pb Age Determination for Seven Standard Zircons using Inductively Coupled Plasma–Mass Spectrometry Coupled with Frequency Quintupled Nd-YAG (λ = 213 nm) Laser Ablation System: Comparison with LA-ICP-MS Zircon Analyses with a NIST Glass Reference Material

This paper evaluates the analytical precision, accuracy and long‐term reliability of the U‐Pb age data obtained using inductively coupled plasma–mass spectrometry (ICP‐MS) with a frequency quintupled Nd‐YAG (λ = 213nm) laser ablation system. The U‐Pb age data for seven standard zircons of various ages, from 28 Ma to 2400 Ma (FCT, SL13, 91500, AS3, FC1, QGNG and PMA7) were obtained with an ablation pit size of 30 μm diameter. For ²⁰⁷Pb/²⁰⁶Pb ratio measurement, the mean isotopic ratio obtained on National Institute of Standards and Technology (NIST) SRM610 over 4 months was 0.9105 ± 0.0014 (n = 280, 95% confidence), which agrees well with the published value of 0.9096. The time‐profile of Pb/U ratios during single spot ablation showed no significant difference in shape from NIST SRM610 and 91500 zircon standards. These results encouraged the use of the glass standard as a calibration standard for the Pb/U ratio determination for zircons with shorter wavelength (λ = 213 nm) laser ablation. But ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ages obtained by this method for seven zircon standards are systematically younger than the published U‐Pb ages obtained by both isotope dilution–thermal ionization mass spectrometry (ID‐TIMS) and sensitive high‐resolution ion‐microprobe (SHRIMP). Greater discrepancies (3–4% younger ages) were found for the ²⁰⁶Pb/²³⁸U ages for SL13, AS3 and 91500 zircons. The origin of the differences could be heterogeneity in Pb/U ratio on SRM610 between the different disks, but a matrix effect accuracy either in the ICP ion source or in the ablation‐transport processes of the sample aerosols cannot be neglected. When the ²⁰⁶Pb/²³⁸U (= 0.2302) newly defined in the present study is used, the measured ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ages for the seven zircon standards are in good agreement with those from ID‐TIMS and SHRIMP within ±2%. This suggests that SRM610 glass standard is suitable for ICP‐MS with laser ablation sampling (LA‐ICP‐MS) zircon analysis, but it is necessary to determine the correction factor for ²⁰⁶Pb/²³⁸U by measuring several zircon standards in individual laboratories.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

抓“热点”、解“难点”，依法行政搞矿管

原达县地区是全省矿业活动最频繁的地区之一,也一度成了全省有名的“采煤混战”大区。近几年来,以《矿产资源法》为指针,抓住“重点”,调处“热点”,解决“难点”,以点带面,加强矿管行政执法工作,促进了全区矿业秩序的全面好转,矿业生产也相应得到健康发展。矿业产值占全区工业总产值的１５％左右,是全区的一大支柱产业,有力地促进了地方经济的发展。     

地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。