A data-driven approach to local gravity field modelling using spherical radial basis functions   总被引：3，自引：0，他引：3  

R. Klees  R. Tenzer  I. Prutkin  T. Wittwer 《Journal of Geodesy》2008,82(8):457-471

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.  相似文献   

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,19(2):16-22

2006年7月16日娃娃沟流域暴发的大规模泥石流,给下游3个电站造成巨大经济损失,是大渡河流域一次典型的灾害性泥石流。分析得出,娃娃沟泥石流重度高、搬运能力强,泥石流固体物质砂、石混杂,粗大砾石含量高;暴发频率低、规模大,流速及峰值流量分别高达10．78m／s及798．5m^3／s;在汇口处,泥石流堆积物堵塞河道是引起下游电站受灾的重要原因,高重度、粗颗粒、大流量的组合是此次泥石流堵江的重要原因。堵河判别计算结果显示在发生百年一遇泥石流时,该断面均有发生堵河的可能。娃娃沟泥石流表明：①在大渡河支流的泥石流沟周边的中小电站极有可能在泥石流暴发时受到破坏。因此,电站建设过程中应加强对周边泥石流沟的防灾减灾工作;②虽然娃娃沟流域植被良好,但仍然发生了大规模泥石流。表明植被不能完全避免泥石流的发生,对于此类泥石流沟不能疏忽大意。  相似文献   

浅谈煤田地震勘探中影响纵向分辨率的主要因素     

何黄生刘茂争  马玉生张玉良　沈骥千 《中国煤炭地质》2008,20(5):61-64

针对野外数据采集参数对分辨率的影响,对资料处理过程中静校正、反褶积、速度分析、Q补偿及谱白化等几个重要模块原理进行分析,指出在高速岩层中数据采集,应采用小药量激发,小道距,小排列,全三维,不组合、单点数字检波器接收。高密度多点采集,能够充分保护有效地质信息的高频成份,在更大范围内满足高分辨地震勘探的要求。资料处理过程中,在保证同相叠加和高信噪比的前提下,采用提高分辨率的模块以弥补地震波在传播过程中损失的高频成份。实践证明通过高保真、宽频带数据采集和高信噪比前提下的高分辨率数据处理,完全能够为矿井采掘提供高精度地震勘探资料。  相似文献   

南海海盆区莫霍面分布规律及其对深部钻探的意义     

秦绪文  张宝金  赵斌  路允乾  陈玺  王利杰  许振强  张如伟  耿明会  杨振  李建平  吕文超  尉建功 《地质学报》2022,96(8):2635-2646

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方，目前有研究推测南海是世界上莫霍面深度最浅的海域之一，但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态，是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析，结合前人的研究，形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张，岩浆供应充足，受扩张停止后岩浆活动影响较小，基底平坦，地质构造相对简单，同时洋壳地震速度结构不存在异常，且有较强的广角莫霍面反射波和可识别的地幔顶部折射波，具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰，西北次海盆次之，西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km，首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比，利用深反射多道地震计算的莫霍面深度细节更为丰富，并且可以在垂向上清晰刻画莫霍面的结构。整体上，南海海盆区莫霍面地震反射强烈和可信度高的区域中，深度较浅的区域之一是东部次海盆南部，最浅处仅约9. 5 km，其中水深4. 01 km，洋壳厚度仅5. 54 km。综合判断，东部次海盆南部是南海重要的莫霍面钻探备选区，这对南海莫霍面钻探选址具有重要意义。  相似文献   

含煤层地质环境下地震波场的数值模拟   总被引：3，自引：2，他引：1  

潘冬明  杨顶辉 《地球物理学进展》1997,12(1):74-83

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。  相似文献   

Scattering of elastic waves by a fracture zone containing randomly distributed cracks     

Jun Kawahara  Teruo Yamashita 《Pure and Applied Geophysics》1992,139(1):121-144

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.  相似文献   

Transmission and reflection of Rayleigh waves at a thin low velocity vertical layer: A laboratory and theoretical study     

P. F. Biagi  T. S. Gegechkory  E. N. Its  P. Manjgaladze  V. Sgrigna  D. Zilpimiani 《Pure and Applied Geophysics》1990,133(2):317-327

A three-dimensional ultrasonic modelling was used to study the Rayleigh waves reflected from a thin low velocity vertical layer and those transmitted through the same zone. The reflection and transmission coefficients of the Rayleigh waves have been calculated by an approximate method based on the Green function technique. Nonrigid contact boundary conditions were used for the low velocity layer modelling. A rather good correspondence between the theoretical and experimental absolute values of the reflection and the transmission coefficients was obtained. Some discrepancies between experimental and theoretical results, in particular with regard to phase shifts, may be explained by a certain inadequacy of the theoretical model to the low velocity layer investigated by the ultrasonic modelling.  相似文献   

Continuous air rate measurement in flotation cells: Some fundamental considerations     

J.A. Torrealba-Vargas  J.A. Finch 《International Journal of Mineral Processing》2006

A robust and reliable sensor to measure gas (air) superficial velocity (J_g) continuously in flotation systems is introduced. It is based on the sampling of bubbles by buoyancy with the accumulating air allowed to exit through an orifice. At steady state, pressure drop is measured and related to the J_g by prior calibration. The continuous device is readily automated and extended to a multi-unit set-up. The sensor and data collections are described. The governing gas flow equation and models of the dynamic response to air flow rate set point change and fluctuations in froth depth are derived. Model predictions are confirmed against plant data. By taking a moving average, the J_g remains valid in the face of typical plant disturbances.  相似文献   

Mais comment s'écoule donc un glacier ? Aperçu historique     

Frédérique Rémy  Laurent Testut 《Comptes Rendus Geoscience》2006,338(5):368-385

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).  相似文献