利用Jason-1卫星确定全球平均海面高

介绍了卫星测高数据处理的基本方法,分析研究了卫星测高的主要误差,提出了Jason-1的数据编辑准则,并用其测高数据确定了全球平均海面高,与CLS01模型进行了比较分析.     

矿井涌水量混沌时间序列分析与预测

矿井地下水系统是个非线性系统,由于地下水系统的复杂性,很难完全确定影响矿井涌水量的因素,因此传统的矿井涌水预测表现出较大的局限性。而混沌时间序列分析能够利用单变量处理方法,分析矿井地下水系统中,由于多因素耦合作用所产生的矿井涌水量时间序列。通过对某矿区1985年9月至2005年2月矿井涌水量时间序列资料的分析,显示矿井涌水量时间序列具有混沌特征;通过对比预测值和实测值,表明预测精度高。因此,混沌时间序列分析方法用于矿井涌水量预测是可行的。     

GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals

This paper presents a general modeling strategy for ambiguity resolution (AR) and position estimation (PE) using three or more phase-based ranging signals from a global navigation satellite system (GNSS). The proposed strategy will identify three best “virtual” signals to allow for more reliable AR under certain observational conditions characterized by ionospheric and tropospheric delay variability, level of phase noise and orbit accuracy. The selected virtual signals suffer from minimal or relatively low ionospheric effects, and thus are known as ionosphere-reduced virtual signals. As a result, the ionospheric parameters in the geometry-based observational models can be eliminated for long baselines, typically those of length tens to hundreds of kilometres. The proposed modeling comprises three major steps. Step 1 is the geometry-free determination of the extra-widelane (EWL) formed between the two closest L-band carrier measurements, directly from the two corresponding code measurements. Step 2 forms the second EWL signal and resolves the integer ambiguity with a geometry-based estimator alone or together with the first EWL. This is followed by a procedure to correct for the first-order ionospheric delay using the two ambiguity-fixed widelane (WL) signals derived from the integer-fixed EWL signals. Step 3 finds an independent narrow-lane (NL) signal, which is used together with a refined WL to resolve NL ambiguity with geometry-based integer estimation and search algorithms. As a result, the above two AR processes performed with WL/NL and EWL/WL signals respectively, either in sequence or in parallel, can support real time kinematic (RTK) positioning over baselines of tens to hundreds of kilometres, thus enabling centimetre-to-decimentre positioning at the local, regional and even global scales in the future.     

动态最小二乘支持向量机学习算法

针对大样本集的训练问题和动态训练样本的模型更新问题,提出了动态最小二乘支持向量机学习算法.该算法充分利用已建好的模型,逐渐加入新样本,并可删除位于任何位置的非支持向量,避免了矩阵求逆运算,保证了算法的高效率.大坝变形及电离层延迟两个时间序列的预报实例表明,该算法具有计算时间短、预报精度高的特点.     

葛店矿双庙扩大区三维地震勘探

葛店煤矿双庙扩大区位于豫、皖两省交界处。勘探区为一轴向北北东向的不对称向斜构造,地层倾角26°～30°,测区深部有岩浆岩侵入体,严重影响煤层赋存形态;地表条件复杂,村庄密集、河流较多,测区东侧塌陷区水塘大面积分布．三维地震勘探野外施工难度大。本区地震勘探观测系统设计为12线18炮制,为避免偏移归位及叠加次数降低等问题,在向斜轴部两侧、村子周边及河岸处适当加密炮点,并在塌陷区利用钢钎插置检波器以保证CDP点完整性。数据处理采取了叠前反褶积、叠前时间偏移等技术,数据处理效果显著。遵循多波组对比、多个数据体综合解释的原则,以纵向、横向和任意时间剖面相结合,充分利用三维可视化技术进行地震资料解释。本次勘探成果查明了主采煤层三2^2、二2煤层的构造形态及埋藏深度;解释落差5m以上断层34条,其中修改原有断层5条;另外新发现小于5m的异常断点22个。经巷道验证煤层底板标高误差0．43％～1．96％,巷道已揭露断层7条,与解释成果基本一致。     

第二次土地调查若干问题探讨

第二次土地调查的目的是全面查清全国土地利用状况,掌握真实的土地基础数据,如何保证调查的质量?文章以江阴市为例进行探讨。     

Group Time Delay in the Reflection of P-Polarized Electromagnetic Plane Waves at a Stratigraphic Interface

The group time delay (GTD) formulae of quasi-total reflection of inhomogeneous P-polarized electromagnetic waves (P-PEW, whose electric field parallels to the plane of incidence), caused by the Goos–Hänchen effect, is derived by the phase shift of the wave. A numerical example where the frequency equals to 1 GHz is given. The curves of calculating results show that there are three discontinuous points at the critical angle of phase shift vector, the critical angle of attenuation vector and 90^○. When the angle of incidence is equal to one of these three angles, the GTD will become infinite, and the electromagnetic wave will propagate along the interface. The GTD will be very large, when the angle of incidence is close to one of these three angles. The results indicate that there is an important relationship between the two critical angles and the conductivity and permittivity of the two strata, and that the GTD has an important relationship with the angle of incidence. These results can be used to determine the lithology of the strata and to delimit the interfaces of strata more effectively. It is suggested that this research may prove useful in electromagnetic logging analysis and, perhaps, in the design of logging instruments.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

随钻测录井地质导向二维分解实时绘图方法

根据随钻测录井实时地质导向和大斜度、水平井评价成图技术需求，针对传统绘图方法存在的弊端，提出了将测录井信息、井眼轨迹和地质模型进行二维分解的实时绘图方法。针对二维分解绘图方法绘图时空复杂度较高的问题，给出了不同事件驱动下的局部实时计算和拷屏重绘算法，控制了对CPU和内存的消耗、提高了绘图效率，消除了实时绘图的闪烁和卡顿现象。应用实例表明，二维分解实时绘图方法能够实现大尺度随钻测录井地质导向图形的流畅、无卡顿实时绘图，可提高大斜度、水平井储层模型评价的刻画精度和时效。     

中国大陆第四纪地壳运动时程

本文从中国大陆有代表性的青藏高原、黄土高原、三峡地区着手,进行了新构造运动剖析。并结合全国火山活动期次,总结出中国大陆第四纪地壳运动的时程变化。从而为进一步探讨新构造运动特点和规律,为国民经济建设及人类生存环境提供了一定的地质背景资料。