蒙古中部剖面重力异常研究

在蒙古中部106°E附近布设一条南北走向的重力剖面,首次获得该地区精确的实测重力值,经过数据预处理和改正计算,得到剖面重力异常结果。分析表明,剖面自由空气异常和地形有很强的相关性|剖面布格重力异常变化总体较为平缓,只在戈壁-阿尔泰-曼达尔戈壁加里东区带(GAB)和哈拉加里东区带（HRB）附近发生相对较快的变化,推测这两个区域可能存在较明显的深部构造变化|在蒙古-鄂霍兹克造山带附近,并未发现与板块缝合带相对应的重力异常变化|剖面均衡重力异常以负值为主,多数地区地壳未达到重力均衡状态,剖面附近地震活动和均衡重力异常之间的关系与中国南北地震带北段相关研究结果类似。     

一种基于丹麦法的改进型双步M估计

针对抗差估计方法初值与临界值选取不当的问题,提出一种以L1-范估计的平差值作为初值、以改进的丹麦法权函数作为第一步抗差阶段权函数的双步M估计,该方法既具有较强的抗差性,又具有最小二乘的最优性。试验表明,该方法估计结果与无粗差时的LS估计结果基本一致,抗差效果明显优于传统抗差估计方法与其他两步抗差估计。     

陆地惯性重力测量静态试验研究

利用高精度捷联式惯性测量设备开展测定重力场的初步静态试验,并对静态惯性重力测量影响因素作初步探讨。静态测试表明,基于试验采用IMU内置的石英挠性加速度计,在稳定的观测环境下通过30 s的滑动平均滤波可以消除高频噪声,得到标准差约为1.5 mGal的重力观测值,但仪器自身的不稳定性和环境变化引起的漂移会使重力观测值呈线性增大的趋势。针对该变化,通过回归分析进行拟合,并结合初始重力已知点上观测得到的初始偏移量,在车载试验中对零速修正点上的静态重力测量误差进行补偿,使大多数点上的重力观测精度由12~20 mGal提升至优于2 mGal。     

A theoretical study on terrestrial gravimetric data refinement by Earth gravity models

The idea of this paper is to present estimators for combining terrestrial gravity data with Earth gravity models and produce a high‐quality source of the Earth's gravity field data through all wavelengths. To do so, integral and point‐wise estimators are mathematically developed, based on the spectral combination theory, in such a way that they combine terrestrial data with one and/or two Earth gravity models. The integral estimators are developed so that they become biased or unbiased to a priori information. For testing the quality of the estimators, their global mean square errors are generated using an Earth gravity model08 model and one of the recent products of the gravity field and steady‐state ocean circulation explorer mission. Numerical results show that the integral estimators have smaller global root mean square errors than the point‐wise ones but they are not efficient practically. The integral estimator of the biased type is the most suited due to its smallest global root mean square error comparing to the rest of the estimators. Due largely to the omission errors of Earth gravity models the point‐wise estimators are not sensitive to the Earth gravity model commission error; therefore, the use of high‐degree Earth gravity models is very influential for reduction of their root mean square errors. Also it is shown that the use of the ocean circulation explorer Earth gravity model does not significantly reduce the root mean square errors of the presented estimators in the presence of Earth gravity model08. All estimators are applied in the region of Fennoscandia and a cap size of 2° for numerical integration and a maximum degree of 2500 for generation of band‐limited kernels are found suitable for the integral estimators.     

Two-dimensional physical-based inversion of confined and unconfined aquifers under unknown boundary conditions

An inverse method is developed to simultaneously estimate multiple hydraulic conductivities, source/sink strengths, and boundary conditions, for two-dimensional confined and unconfined aquifers under non-pumping or pumping conditions. The method incorporates noisy observed data (hydraulic heads, groundwater fluxes, or well rates) at measurement locations. With a set of hybrid formulations, given sufficient measurement data, the method yields well-posed systems of equations that can be solved efficiently via nonlinear optimization. The solution is stable when measurement errors are increased. The method is successfully tested on problems with regular and irregular geometries, different heterogeneity patterns and variances (maximum K_max/K_min tested is 10,000), and error magnitudes. Under non-pumping conditions, when error-free observed data are used, the estimated conductivities and recharge rates are accurate within 8% of the true values. When data contain increasing errors, the estimated parameters become less accurate, as expected. For problems where the underlying parameter variation is unknown, equivalent conductivities and average recharge rates can be estimated. Under pumping (and/or injection) conditions, a hybrid formulation is developed to address these local source/sink effects, while different types of boundary conditions can also exert significant influences on drawdowns. Local grid refinement near wells is not needed to obtain accurate results, thus inversion is successful with coarse inverse grids, leading to high computation efficiency. Furthermore, flux measurements are not needed for the inversion to succeed; data requirement of the method is thus not much different from that of interpreting classic well tests. Finally, inversion accuracy is not sensitive to the degree of nonlinearity of the flow equations. Performance of the inverse method for confined and unconfined aquifer problems is similar in terms of the accuracy of the estimated parameters, the recovered head fields, and the solver speed.     

山东洞体应变观测质量分析

选取山东省3个洞体应变台2007-2011年观测数据,分析资料的完整性、内在质量及干扰因素.结果表明：①观测资料连续性和完整性较好;②仪器零漂满足设计指标要求;③M2波振幅因子相对误差和相对噪声水平达到Ⅰ类应变台标准;④仪器更新后观测质量提升;⑤洞体应变观测主要受仪器调修、气压、降雨和观测环境变化等因素影响.     

箱包CT检测系统弱信号情况下的CT成像实验研究

在箱包CT检测中,当放射源活度很低的情况时,对CT成像进行了实验研究。利用多次累加的方法降低统计涨落和噪声的影响,使CT图像质量得到提高。实验证明对于弱信号情况下,利用多次累加的方法,CT成像具有比较高的质量,能够实现物性判别的要求。另外对于A/D量化误差也进行了分析。     

Upper mantle seismic anisotropy beneath a convergent boundary: SKS waveform modeling in central Tibet

This study used SKS waveforms from the International Deep Profiling of Tibet and the Himalayas (INDEPTH) III dataset and a new 2D method for modeling seismic waves in anisotropic media to construct an image of anisotropic structures beneath central Tibet. A preferred model revealed three-segment anisotropic structures in the upper mantle beneath the study region. Waveform modeling demonstrated that the anisotropy was mainly generated by the lithosphere but not the asthenosphere, and that an anisotropic model with a flatter axis of symmetry provides a more consistent interpretation of the observations than models having steeply dipping symmetry axes. A relatively low velocity zone may underlie or intermingle with the anisotropic structures in the northern portion of the region. Synthetic tests also indicate that variations in the elastic constants and depth extent of the anisotropy assumed by the calculations do not affect the general conclusions, although trade-offs exist among certain model parameters. The modeling results suggest that the complex seismic structures in central Tibet were associated with underthrusting of the Indian lithosphere beneath the Asian lithosphere; the inferred flat symmetry axis of the anisotropy was likely generated during this collision process. If this were not the case, the inherited anisotropy would exhibit a steeply dipping axis of symmetry, parallel to the direction of underthrusting.     

黑龙江省三江平原垦区中心城市选择

黑龙江省三江平原垦区是中国耕地规模最大的连片农垦区和第三大区际商品粮供应基地。近年来,随着这一地区工业化和城镇化进程的加快,如何科学的选择和建设区域的中心城市成为该区完善区域城镇体系和推进城镇化进程的紧迫任务和现实课题。基于此,应用改进的引力模型来计算县级单位与农场间的空间联系强度,界定出各县级单位的腹地范围,明确三江平原垦区的中心城市;以各农垦分局为研究单元,研究各农垦分局的碎化程度;通过明确县级单位拥有农场腹地数目,确定各农垦分局的中心城市;结果发现建三江、虎林、富锦、友谊和密山等宜作为三江平原垦区的中心城市;提出通过壮大中心城市的综合实力、推动场县共建和深化对外合作等措施来推动三江平原垦区城镇化进程。     

中天山天湖东铁钼矿含矿片麻状花岗岩年代学、地球化学和锆石Hf同位素-对于中天山早古生代构造演化的启示

以中天山东段的天湖东铁钼矿含矿花岗岩为例,在LA-ICP-MS锆石U-Pb年代学测定其为早古生代花岗岩（（445.3±4.6） Ma）基础上,通过岩相学、地球化学及锆石原位Hf同位素组成等多方面研究,探讨该岩体的成岩作用及其构造背景。天湖东含矿片麻状花岗岩的主要矿物为斜长石、石英、钾长石,并含少量黑云母和角闪石等。全岩地球化学分析结果表明,该片麻状花岗岩高硅、弱富铝、富钙、富钠而贫钾,ASI值为0.68~0.82,属于准铝质钙碱性花岗岩,总体上富集大离子亲石元素（LILE）Rb、Ba等和轻稀土元素La、Ce、Nd等,而亏损高场强元素Nb、Ta、Ti、Yb等,轻重稀土分异明显,轻稀土分异较为明显,而重稀土分异不明显,表现出典型岛弧岩浆岩的地球化学特征。锆石的ε_Hf(445 Ma)值为-6.31~-1.77,二阶段Hf模式年龄（T_DM2）为1.538~1.825 Ga,表明该花岗岩的源区主要为壳源物质。综合分析上述资料,认为天湖东铁钼矿片麻状花岗岩是由俯冲过程中地壳物质重熔的产物。结合前人的研究和本课题组的新近研究成果认为,在早古生代时,中天山为岩浆弧构造环境,形成一系列的钙碱性岩浆岩,而该岩浆弧的形成可能是受到介于吐哈陆块和塔里木板块之间的古天山大洋在早古生代时期向南俯冲而形成的。