Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

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.     

山西沁水盆地煤层气定向井钻井工艺

山西沁水盆地蕴藏着丰富的煤层气资源,采用定向井开采煤层气具有开采成本低、煤层气产量高的优点。本文结合山西省沁水县名水村的工程实例,介绍了定向井的钻井工艺,总结了施工经验,对类似工程具有借鉴意义。     

中国西部现今地壳形变与地震活动特征

以1999-2004年间中国西部地区近650个GPS站点的观测资料为基础,采用二维“高张力样条”函数内插算法获得了连续地壳形变场。结果表明,相对于稳定欧亚参考框架,中国西部现今地壳运动西强东弱,北向运动从西向东逐渐减弱,东向运动逐渐增强,青藏高原东缘及附近地区是东向运动的消减区带;主应变从西向东由压缩应变转变为拉伸应变;最大剪应变主要位于喜马拉雅地块及其东部地区;中国西部地区地震活动在空间分布上似乎集中发生在印度板块北东向挤压欧亚板块的两条共扼带内,与地壳形变场有一定对应关系,地壳形变的高应变率区为中强地震多发地带。     

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

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.     

兰州—白银多目标地球化学调查野外工作方法浅析

多目标区域地球化学调查从调查区表层、深层采样土壤元素分布特征入手建立土壤元素地球化学分配机制,以景观、地貌及地质特征为基础,选择合理的采样位置,正确的采样方法、采样层位,应用掌上机（Pocket PC）及数字填图系统（Rgmap）采集野外样品,本文以兰州—白银黄土覆盖区景观区为例,对表层、深层采样方法和技术要求进行探讨,对全省开展多目标地球化学调查,特别是野外工作方法提供借鉴。     

河西走廊疏勒河流域地质灾害危险性现状评估方法的应用——以甘肃省河西走廊(疏勒河)农业灌溉暨移民安置综合开发建设项目8个移民社区及配套工程建设为例

依据“建设用地地质灾害危险性评估技术要求”、利用野外踏勘及相关资料,对“甘肃省河西走廊（疏勒河）农业灌溉暨移民安置综合开发建设项目8个移民社区及配套工程”进行地质灾害危险性现状评估。结果认为：评估区发育和分布的地质灾害类型有土壤盐渍化、风蚀沙埋和洪水冲蚀。现状评估结果：土壤盐渍化和洪水冲蚀危险性均小;风蚀沙埋有1处危险性中等,其它危险性小。     

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

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

利用参考重力场模型基于能量法确定GRACE加速度计校准参数

利用多个参考重力场模型分别对GRACE一个月的实测加速度计观测数据进行检校.数值计算结果的比较分析表明了利用参考重力场模型确定加速度计校准参数是有效的.     

地下矿产开采对天然气管道工程的安全影响评价

西气东输天然气管道工程途经兖州市,穿越杨庄、杨村、兴隆庄等煤矿矿区,压覆煤层埋深变化较大。该文采用预测法对各煤矿在不同开采条件下分煤层对管道运营安全的危害程度进行了评价分析。通过计算表明,煤矿开采埋深越浅,煤层越厚,对管道危害越严重。兴隆庄煤矿3煤层、杨村煤矿煤16t、煤17开采对管道影响强烈,杨庄煤矿煤16t、煤17开采对管道影响中等,兴隆庄煤矿煤16t、煤17开采对管道运营影响小。     

广西佛子冲铅锌矿田成岩成矿时代研究

广西佛子冲铅锌矿田位于桂东南博白-岑溪多金属成矿带的NE端,区内发育燕山期花岗岩、花岗闪长岩、二长花岗斑岩和英安质凝灰熔岩等;其中花岗闪长岩和二长花岗斑岩与铅锌成矿关系密切,成矿岩体形成时代分别为156Ma和127Ma。