断层垂向封闭性定量研究方法及其在准噶尔盆地白家海凸起东道海子断裂带应用   总被引：2，自引：0，他引：2  

梁全胜  刘震  何小胡  王德杰  王晓东 《现代地质》2008,22(5)

断层的垂向封闭能力主要取决于断面的紧闭程度和断裂带内填充物的岩性。在定量计算断面正压力和断层两侧地层泥地比的基础上,提出了断层垂向封闭因子(Fvs)的概念,并将其定义为断面正压力与断层两侧地层泥地比的乘积。通过该参数可以定量评价断层的垂向封闭能力。以准噶尔盆地东道海子断裂为例,在建立断层垂向封油、气能力评价标准的基础上,应用断层垂向封闭因子(Fvs)定量评价了该断裂带的成藏期垂向封闭性及现今垂向封闭性,分析了该断裂带对油气成藏的控制作用,得到较好效果。  相似文献   

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C   总被引：2，自引：1，他引：1  

Christoph Förste  Roland Schmidt  Richard Stubenvoll  Frank Flechtner  Ulrich Meyer  Rolf König  Hans Neumayer  Richard Biancale  Jean-Michel Lemoine  Sean Bruinsma  Sylvain Loyer  Franz Barthelmes  Saskia Esselborn 《Journal of Geodesy》2008,82(6):331-346

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

Methodology and use of tensor invariants for satellite gravity gradiometry   总被引：2，自引：1，他引：1  

Oliver Baur  Nico Sneeuw  Erik W. Grafarend 《Journal of Geodesy》2008,82(4-5):279-293

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.  相似文献   

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition   总被引：3，自引：0，他引：3  

Oliver Baur  Gerrit Austen  Jürgen Kusche 《Journal of Geodesy》2008,82(4-5):207-221

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.  相似文献   

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry     

Wenbin Shen  Jin Li  Jiancheng Li  Zhengtao Wang  Jinsheng Ning  Dingbo Chao 《地球空间信息科学学报》2008,11(4):273-278

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.  相似文献   

斜坡应力分布的系统测试与反演分析     

沈军辉  崔建凯  徐进  贾留杰 《工程地质学报》2008,16(5):616-619

采用简易、快速、经济的洞壁应力恢复法对某工程边坡勘探平硐洞壁应力场作了系统测试, 并采用FLAC2D软件对测试成果作了斜坡初始应力的反演分析。结果表明, 斜坡应力场随深度变化具有明显分带性, 应力量值随硐深曲线具有典型的驼峰式特征, 与理论分析结果相符, 表明该方法能较好地应用于斜坡应力场的系统测试和研究。成果对认识斜坡岩体的变形破裂及卸荷特征等具有重要意义。  相似文献   

Time-dependent stress change during failure of rocks     

V.T. Belikov  A.F. Shestakov 《Russian Geology and Geophysics》2008,49(5):350-356

The energy balance of a solid subject to fracture has been explored using heat and mass transfer equations with regard to the volumetric and superficial components. In the suggested model, brittle fracture of a cracked solid considered as a heterogeneous two-phase medium is described by an equation analogous to the Griffith’s criterion for propagation of a single crack. The derived equation is used, together with estimates of relative change in specific interface area, to study the respective change of free strain energy and pressure in rocks associated with failure.  相似文献   

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

西藏札达盆地控盆断裂有限元数值模拟   总被引：1，自引：0，他引：1  

段岩  孟宪刚  邵兆刚  朱大岗 《地质通报》2008,27(10)

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

On　the　physical　model　of　earthquake　precursor　fields　and　the　mechanism　of　precursors＇timespace　distribution──origin　and　evidences　of　the　strong　body　earthquake──generating　model     

梅世蓉 《地震学报(英文版)》1995,(3)

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...  相似文献