Marquardt算法磁测三维反演方法策略   总被引：3，自引：0，他引：3  

李军  李才明  王邦华  侯智勇 《物探化探计算技术》2004,26(4):320-324

在磁法三维反演中,反演对象数目往往巨大,受限于当今的计算机技术与计算方法,在超大规模方程组的求解时遇到极大的困难。所以,在解决实际问题时能力非常有限,对于Marquardt算法磁测三维反演问题尤为突出。为此,在这里提出了两个较为有效的方法策略,即将模型组体化和将组体分类,并在理论应用中得到了较好的效果。  相似文献   

基于"3S"技术的地质灾害易损性面评价方法研究   总被引：4，自引：0，他引：4  

杨华  赵纯勇  张永兴 《中国地质灾害与防治学报》2004,15(3):96-101,125

易损性评价是地质灾害灾情评估体系中一个非常重要的内容。它揭示了地质灾害基本构成元素之一的承灾体的社会经济属性,对整个地灾灾情评估体系完整性的构建及评价结果的科学性、明确性、统一性表达具有举足轻重的作用。文章以等面积法作为基本手段,根据评价区实际情况算得的离散等积网格作为评价模板。结合“3S”(GIS、RS、GPS)技术,探讨地质灾害易损性面评价的实现方式,旨在提出“3S”框架内的地质灾害灾情评估实现模式。在具体操作层面上,文章提出了为实现易损性面评估而确立的评价对象之层次分类方法、评价模板构造计算方法、空间数据库建设方法、基于空间数据库的空间分析方法以及空间信息集成模型。重点讨论在空间数据完整性基础上的多源数据融合处理模式及以GIS空间分析算法模式为蓝本的空间数据处理和集成手段。在理念层面上,文章将空间分析实现方法分为现实评价区域、概念评价区域、数字评价区域、逻辑映射虚拟区域等4个阶段,认为在GIS环境下的地灾灾情评价原始空间数据、派生空间数据以及经专用模型集成后的结果空间数据的可视化表达方式可以分为颜色序列、数字高程模型(DTM)及等值线3种表达方式。最后提出实现地质灾害易损性面评价的一般操作流程。  相似文献   

A technique for modeling the Earth’s gravity field on the basis of satellite accelerations   总被引：3，自引：2，他引：1  

P. Ditmar  A. A. van Eck van der Sluijs 《Journal of Geodesy》2004,78(1-2):12-33

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.  相似文献   

Building and Testing an <Emphasis Type="Italic">a priori</Emphasis> Geophysical Model for Western Eurasia and North Africa     

Michael E.?PasyanosEmail author  William R.?Walter  Megan P.?Flanagan  Peter?Goldstein  Joydeep?Bhattacharyya 《Pure and Applied Geophysics》2004,161(2):235-281

We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q _P and Q _S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle _n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.  相似文献   

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea   总被引：2，自引：0，他引：2  

A. M.?DaviesEmail author  P.?Hall  M. J.?Howarth  P. J.?Knight  R. J.?Player 《Ocean Dynamics》2004,54(2):108-125

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke  相似文献   

Stable carbon isotope ratios in lake and swamp sediments as a proxy for prehistoric forest clearance and crop cultivation in the Neotropics   总被引：2，自引：0，他引：2  

Chad S. Lane  Sally P. Horn  Claudia I. Mora 《Journal of Paleolimnology》2004,32(4):375-381

Close correspondence between stable carbon isotope ratios ( ¹³ C), pollen, and charcoal profiles in sediment cores from Laguna Zoncho and Machita swamp, Costa Rica, shows that prehistoric forest clearance and crop cultivation can be detected in the stable carbon isotope ratios of total organic carbon ( ¹3C _TOC ). Analyses of δ 13C _TOC complement evidence from pollen, charcoal, and phytoliths and provide a proxy that is sensitive to the intensity and/or proximity to core sites of prehistoric forest clearance and agriculture in watersheds. Stable carbon isotope analyses are particularly useful in situations in which other evidence of forest clearance and agriculture is limited.  相似文献   

Effects of the number of stations and time resolution on derivation     

Odim Mendes  Jr.  Aracy Mendes da Costa  Fernando Celso Perin Bertoni 《Journal of Atmospheric and Solar》2006,68(18):2127-2137

For at least 30 years now it has been well known that the Dst index can be modelled using the solar wind as input. Since then, many attempts have been made to improve the predictability of Dst using different approaches. These attempts are useful, for instance, to understand which features of the solar wind–magnetosphere interactions are most important in producing magnetospheric activity and how the Dst index would improve the space weather forecast. The Dst index is by far the most reliable and simple indication that a magnetic storm is in progress. In this work, the effect of using more than four magnetic stations and shorter time intervals than the hourly averages used in Sugiura's procedures is evaluated. The discussion is based on the results presented by Burton in 1975 and Feldstein in 1984 considering 4 or 12 magnetic stations and time averages of 2.5 min for a magnetic disturbed period that occurred from February 7–28, 1967, including two geomagnetic storms. The analysis has shown that the global representation of a magnetic storm by the standard Dst (Sugiura) is well preserved either using 4, 6, 12 magnetic stations or using 1 h, 2.5 min 1 min averages. A brief review of the current understanding of Dst has been included to support the discussions. The analysis performed has shown that a more refined Dst index (time and number of stations>4) would be useful to investigate the intrinsic processes and the different current systems involved in the ring current development during magnetic storms; the standard Dst, as it is conceived, is quite adequate to monitor geomagnetic storms and identify their overall features; concerning the magnetic stations normally considered, the inclusion of higher magnetic latitude stations (>35) may underestimate the observed Dst.  相似文献   

Unstable kinetic Alfvén wave in partially ionized plasma     

J. Vranjes  D. Petrovic  M. Kono 《Planetary and Space Science》2006,54(6):641-644

The stability of kinetic Alfven waves is discussed for a partially ionized plasma with a flux of ionizing electrons which balance the plasma particle losses. Accidental electromagnetic perturbations are shown to be unstable due to the energy change of ionizing electrons.  相似文献   

面向地上下无缝集成建模的新一代三维地理信息系统   总被引：5，自引：0，他引：5  

吴立新  车德福  郭甲腾 《测绘工程》2006,15(5):1-6

地上下无缝集成三维建模是新一代3DGIS的主要标志,地上下无缝集成建模已成为当务之急。在三维空间建模技术现状与存在问题、目标层次与功能需求及三维空间模型现状与趋势分析的基础上,介绍了地上下集成建模的两个层次和无缝集成的基本原理;阐述以CD-TIN为纽带、以BRep-TIN-GTP为核心、以三层混合模型为成份的地上下集成空间数据模型的概念结构与典型逻辑关系。并结合城市与矿山应用,介绍该集成模型在GeoMo3D系统中初步实现后的可视化效果。同时,指出了与地上下集成建模相关的近年3DGIS的主要攻关方向和重点理论难题。  相似文献   

基于3D+-TPR-tree的点目标全时段移动索引设计   总被引：1，自引：1，他引：0  

郭晶  刘广军  郭磊  董绪荣 《测绘学报》2006,35(3):267-272

在经典3D R-tree基础上提出新的3D R-tree索引,通过改变待索引数据项的结构并重新设计查询处理算法,减少包容矩形死区,提高查询效率;为了满足全时段查询要求,设计一种称为3D -TPR-tree的联合索引结构,并对其中TPR-tree的参数包容矩形的调整算法进行优化。通过测试,证明3D R-tree的查询效率明显高于普通3D R-tree;此外,测试结果也表明经过优化的参数包容矩形的调整算法也部分提升了TPR-tree的查询性能。  相似文献