浅层SH反射波法和多道瞬态面波法 进行联合工程勘察的讨论

针对目前单一工程物探方法调查浅部地质结构的局限性,提出采用弹性波法勘察中研究浅部地质结构的SH波反射勘察和多道瞬态面波法勘察这2种有效手段进行联合勘探,提高其可靠性和精度,并给出了3个实例,指出了在数据采集和处理时应注意的问题。     

Benthic nutrient fluxes in the intertidal flat within the Changjiang （Yangtze River） Estuary

In an annual cycle from March 2005 to February 2006, benthic nutrient fluxes were measured monthly in the Dongtan intertidal flat within the Changjiang （Yangtze River） Estuary. Except for NH4^＋, there always showed high fluxes from overlying water into sediment for other four nutrients. Sediments in the high and middle marshes, covered with halophyte and consisting of macrofauna, demonstrated more capabilities of assimilating nutrients from overlying water than the low marsh. Sampling seasons and nutrient concentrations in the overlying water could both exert significant effects on these fluxes. Additionally, according to the model provided by previous study, denitrification rates, that utilizing NO3- transported from overlying water （Dw） in Dongtan sediments, were estimated to be from -16 to 193 μmol·h^-1·m^-2 with an average value of 63 μmol·h^-1·m^-2 （n=18）. These estimated values are still underestimates of the in-situ rates owing to the lack of consideration of DN, i.e., denitrification supported by the local NO3^- production via nitrification.     

Environmental Investigation and Evaluation of Land Subsidence in the Datong Coalfield Based on InSAR Technology

Heavy mining of Jurassic and Carboniferous horizontal coal seams in the Datong coalfield has seriously affected the local geological environment, which is mainly manifested by such geohazards as soil avalanches, landslides, mudflows, surface subsidence, surface cracks, surface solid waste accumulation and surface deformation. More seriously, coal mining causes groundwater to leak. Overpumping of groundwater has resulted in substantial land subsidence of the urban area in Datong City. Based on the previous geo-environmental investigation in the work area, the authors used radar remote sensing techniques such as InSAR （synthetic aperture radar interferometry） and D-InSAR （differential synthetic aperture radar interferometry）, supplemented by the optical remote sensing method, for geo-environmental investigation to ascertain the geo-environmental background of the Datong Jurassic and Carboniferous-Permian coalfield and evaluate the effects of the geohazards, thus providing a basis for the geo-environmental protection, geohazard control and prevention, land improvement and optimization of the human environment. In this study 8 cog-nominal ERS-1/2 SAR data frames during 1992 to 2003 were used, but the following processing was made： （1） the multitemporal SAR magnitude images were used to interpret the geological structure, vegetation, microgeomorphology and drainage system; （2） the multi-temporal InSAR coherent images were used to make a classification of surface features and evaluate the coherence change due to coal mining; and （3） the multi-temporal cog-nominal SAR images were used to complete D-InSAR processing to remove the information of differential deformation areas （sites）. In the end, a ten-year time series of differential interferograms were obtained using the multi-temporal cog-nominal SAR images. In the tests, 84 deformed areas （sites） were obtained, belonging to those in 1993-1996, 1996-1997, 1997-1998, 1998-2001, 1998-2002 and 2001-2003 respectively. Of the 84 areas, 44 are m     

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.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

Prediction of UT1–UTC,LOD and AAM χ<Subscript>3</Subscript> by combination of least-squares and multivariate stochastic methods

This article presents the application of a multivariate prediction technique for predicting universal time (UT1–UTC), length of day (LOD) and the axial component of atmospheric angular momentum (AAM χ ₃). The multivariate predictions of LOD and UT1–UTC are generated by means of the combination of (1) least-squares (LS) extrapolation of models for annual, semiannual, 18.6-year, 9.3-year oscillations and for the linear trend, and (2) multivariate autoregressive (MAR) stochastic prediction of LS residuals (LS + MAR). The MAR technique enables the use of the AAM χ ₃ time-series as the explanatory variable for the computation of LOD or UT1–UTC predictions. In order to evaluate the performance of this approach, two other prediction schemes are also applied: (1) LS extrapolation, (2) combination of LS extrapolation and univariate autoregressive (AR) prediction of LS residuals (LS + AR). The multivariate predictions of AAM χ ₃ data, however, are computed as a combination of the extrapolation of the LS model for annual and semiannual oscillations and the LS + MAR. The AAM χ ₃ predictions are also compared with LS extrapolation and LS + AR prediction. It is shown that the predictions of LOD and UT1–UTC based on LS + MAR taking into account the axial component of AAM are more accurate than the predictions of LOD and UT1–UTC based on LS extrapolation or on LS + AR. In particular, the UT1–UTC predictions based on LS + MAR during El Niño/La Niña events exhibit considerably smaller prediction errors than those calculated by means of LS or LS + AR. The AAM χ ₃ time-series is predicted using LS + MAR with higher accuracy than applying LS extrapolation itself in the case of medium-term predictions (up to 100 days in the future). However, the predictions of AAM χ ₃ reveal the best accuracy for LS + AR.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

多道瑞利波在矿井独头巷道超前探中的应用

瑞利波超前探测作为一种新的矿井地球物理探测方法,具有简便、快速、经济、分辨率高、适用场地小、应用范围广等优点。介绍了瑞利波探测技术原理,并将其用于煤矿井下独头巷道大距离超前探测中。实例表明,探测距离达80 m以上,条件允许时可达100 m以上。     

基于AHP和BP神经网络的深部地热水可持续开发能力评价

采用资料完整性、开采潜力、回灌量、平均水压下降速率、地面沉降速率、水温、水质、地热井布局8项指标构建天津地热可持续开发能力评价指标体系;运用层次分析法确定了各项指标的权重,建立起评价因素集和评语集,给出了归一化数值;建立了天津地热可持续开发能力的BP神经网络模型,以层次分析法得出的结果作为样本,对BP网络进行了训练和测试,实例评价结果表明了AHP和BP神经网络方法的可行性,为地热资源的可持续开发能力评价提供了一种新的评价方法。