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.  相似文献   

The Creation of a National Agricultural Land Use Dataset: Combining Pycnophylactic Interpolation with Dasymetric Mapping Techniques     

Alexis Comber  Chris Proctor  Steve Anthony 《Transactions in GIS》2008,12(6):775-791

Agricultural Census data is summarised over spatially coarse reporting units for reasons of farm confidentiality. This is problematic for research at a local level. This article describes an approach combining dasymetric and volume preserving techniques to create a national land use dataset at 1 km² resolution. The results for an English county are compared with contemporaneous aggregated habitat data. The results show that the accurate estimates of local agricultural land use (Arable and Grass) patterns can be estimated when individual 1 km squares are combined into blocks of > 9 squares, thereby providing local estimates of agricultural land use. This in turn allows more detailed modelling of land uses related to specific livestock and cropping activities. The dataset created by this work has been subject to extensive external validation through its incorporation into a number of other national models: nitrate leaching (e.g. MAGPIE, NEAP‐N), waste, and pathogen modelling related to agricultural activity.  相似文献   

Numerical simulation of a stratigraphic model     

Véronique Gervais  Roland Masson 《Computational Geosciences》2008,12(2):163-179

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.  相似文献   

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.  相似文献   

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula   总被引：7，自引：3，他引：4  

Vicente Medina  Marcel Hürlimann  Allen Bateman 《Landslides》2008,5(1):127-142

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.  相似文献   

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献     

孙升林 《中国煤炭地质》2008,(6):1-5

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!  相似文献   

三维地震小面元采集技术在晋城矿区的应用     

张胤彬 《中国煤炭地质》2008,(6):70-72

依托“西部煤炭资源高精度三维地震勘探技术”工程,对晋城矿区进行了旨在提高小断层,小陷落柱探测能力的高密度三维地震勘探。根据面元选择因素及该区地质任务,采用5m×5m网格进行野外数据采集;考虑炮检距、方位角、覆盖次数、排列片横纵比及煤层埋深（350～500m）等因素,采用中点放炮、60道接收,24次覆盖（横向4次,纵向6次）的8线16炮束状观测系统,基岩中激发。原始资料经同一处理流程后,获得5m×5m×1ms、5m×10m×1ms、10m×10m×1ms及2.5m×2.5m×1ms不同单元的三维数据体多个,通过对比可以发现小断层,小陷落柱在其小面元叠加时间剖面、顺层切片及相干切片都有清晰的反映。实例说明,小面元采集技术可以提高对小构造的纵、横向分辨能力,满足山区对三维地震精确勘探的要求。  相似文献   

Hydrometallurgical processing and recovery of molybdenum trioxide from spent catalyst   总被引：5，自引：0，他引：5  

Kyung Ho Park  B. Ramachandra Reddy  D. Mohapatra  Chul-Woo Nam 《International Journal of Mineral Processing》2006,80(2-4):261-265

Hydrometallurgical processing of spent hydrodesulphurisation (HDS) catalyst for the recovery of molybdenum using sodium carbonate and hydrogen peroxide mixtures was investigated. The results indicated that the recovery of molybdenum was largely dependent on the concentrations of Na₂CO₃ and H₂O₂ in the reaction medium, which controls the acidity of the leach liquor and carry over of impurities such as Al, Ni, P, Si and V. Leaching process was exothermic and leaching efficiency of molybdenum decreased with increasing solid to liquid ratio. Large scale leaching of spent catalyst, under optimum conditions: 20% pulp density, 85 g/L Na₂CO₃, 10 vol.% H₂O₂ and 1 h reaction, resulted a leaching efficiency of 84% Mo. The obtained leach liquor contained (g/L): Mo — 22.0, Ni — 0.015 and Al — 0.82, P — 1.1, Si — 0.094 and minor quantities of V — 8 mg/L, As and Co — < 1 mg/L. Recovery of Mo from leach solution as MoO₃ of 97.30% purity was achieved by ammonium molybdate precipitation method.  相似文献   

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D     

Peter Frolkovič  Jozef Kačur 《Computational Geosciences》2006,10(3):279-290

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).  相似文献   

Statistical analysis for the hydrogeological evaluation of the fracture networks in hard rocks     

J. Hofrichter  G. Winkler 《Environmental Geology》2006,49(6):821-827

The hydrogeological effectiveness of fracture sets is determined and evaluated by the fuzzy c-mean and hierarchical clustering. These cluster analyses combine the geological spatial attributes and the hydraulic relevant attributes of fractures. Based on the results of the clustering the fracture set volumes are estimated.  相似文献