基于三维构造张量的位场边界识别滤波器

地质体构造边界位置的确定是位场数据解释中的一项重要工作,现有很多基于位场梯度张量数据的边界检测滤波器,但存在识别边界位置模糊且无法均衡深浅地质体异常的缺点.本文定义了位场数据的三维构造张量,并提出基于位场构造张量的边界滤波器.为了同时显示不同振幅异常的边界位置,对新定义的滤波器进行归一化处理.在高阶均衡滤波器的计算中需要计算位场的垂向高阶导数,本文引入一种计算的稳定算法,基于拉普拉斯方程利用位场水平导数求解垂向导数,可减小垂向导数计算中产生的误差.将定义的滤波器应用到合成的重磁数据中证明了新方法相比传统的滤波器能更加清晰、准确地圈定边界位置,而且针对同时含有正负异常的地质情况,可避免产生额外的错误边界.最后将新的滤波器应用到实测的重磁数据的解释中,结果显示基于构造张量的滤波器可更准确清晰地划分出断裂的边界位置,发现更多的构造细节.     

基于卫-卫跟踪观测技术利用能量守恒法恢复地球重力场的数值模拟研究

根据卫－卫跟踪观测技术的测量原理,基于能量守恒法建立了一种新的双星相互跟踪和三星相互跟踪的卫星观测方程. 通过数值模拟,采用预处理共轭梯度法恢复120阶地球重力场. 模拟结果表明：第一,双星相互跟踪恢复地球重力场的精度和美国喷气动力实验室公布的EIGEN GRACE02S的结果相符合;第二,三星相互跟踪恢复地球重力场的精度较双星提高约2倍.     

Retrieval algorithm for microwave surface emissivities based on multi-source, remote-sensing data: An assessment on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau plays a very important role in studying severe weather in China and around the globe because of its unique characteristics. Moreover, the surface emissivities of the Qinghai-Tibet Plateau are also important for retrieving surface and atmospheric parameters. In the current study, a retrieval algorithm was developed to retrieve the surface emissivities of the Qinghai-Tibet Plateau. The developed algorithm was derived from the radiative transfer model and was first validated using simulated data from a one-dimensional microwave simulator. The simulated results show good precision. Then, the surface emissivities of the Qinghai-Tibet Plateau were retrieved using brightness temperatures from the advanced microwave-scanning radiometer and atmospheric profile data from the moderate resolution imaging spectroradiometer. Finally, the features of the time and space distribution of the retrieved results were analyzed. In terms of spatial characteristics, a spatial distribution con- sistency was found between the retrieved results and surface coverage types of the Qinghai-Tibet Plateau. In terms of time characteristics, the changes in emissivity, which were within 0.01 for every day, were not evident within a one-month time scale. In addition, surface emissivities are sensitive to rainfall. The reasonability of the retrieved results indicates that the algorithm is feasible. A time-series surface emissivity database on the Qinghai-Tibet Plateau can be built using the developed algorithm, and then other surface or atmospheric parameters would have high retrieval precision to support related geological re- search on the Qinghai-Tibet Plateau.     

瞬变电磁波场转换算法的改进

地球物理电磁场数据与虚拟地震波场数据之间存在数学上的等效转换关系,通过这种等效转换,可有效提高地球物理电磁法对地下目标体分界面的辨识度.但是这种转换在数学上属于不适定问题,可采用奇异值分解法处理.由于大奇异值控制计算矩阵的主要信息,小的奇异值控制计算矩阵的次要信息,传统的截断奇异值分解法只保留大奇异值,而忽略小的奇异值,导致数值解不够精确.本文提出一种新的修正方案——改进截断奇异值法,采用岭估计方法计算由小的奇异值引起的虚拟波场.模型计算结果表明：改进截断奇异值法比传统的奇异值分解法得到的波场转换结果更好,对某煤矿采空区探测数据进行了处理,成功分辨出采空区分界面.

     

Impact of Monex-79 data on the objective analysis of the wind field over the Indian region

The monsoon sub-programme of FGGE deployed a variety of observing systems such as research aircraft, research vessels, floating ballooms, and a geostationary satellite, etc. during the special observational period of Monex-79, with a view to making the best possible data set available for studies on various aspects of the monsoon circulation. Whether the data obtained from the various observational platorms improve the representation of the monsson systems, flow patterns, etc. in the analysis is the basis of this study. For this, the objective analyses of the wind field were made with different data sets, first using only the data from conventional observation platforms, subsequently including the data from aircraft and ships and then finally including the data from the satellite. These analyses were compared with carefully made subjective analyses.It is inferred that the addition of data has a positive impact on the objective analysis and the improvement in the analysis is more marked in the data-void region like the Arabian Sea compared to the region where observational stations are distributed fairly stisfactorily. It is also inferred that the impact of the aircraft/ship data on the analysis is more consistent than the data from the satellite.     

A parallel domain decomposition algorithm for coastal ocean circulation models based on integer linear programming

This paper presents a new parallel domain decomposition algorithm based on integer linear programming (ILP), a mathematical optimization method. To minimize the computation time of coastal ocean circulation models, the ILP decomposition algorithm divides the global domain in local domains with balanced work load according to the number of processors and avoids computations over as many as land grid cells as possible. In addition, it maintains the use of logically rectangular local domains and achieves the exact same results as traditional domain decomposition algorithms (such as Cartesian decomposition). However, the ILP decomposition algorithm may not converge to an exact solution for relatively large domains. To overcome this problem, we developed two ILP decomposition formulations. The first one (complete formulation) has no additional restriction, although it is impractical for large global domains. The second one (feasible) imposes local domains with the same dimensions and looks for the feasibility of such decomposition, which allows much larger global domains. Parallel performance of both ILP formulations is compared to a base Cartesian decomposition by simulating two cases with the newly created parallel version of the Stevens Institute of Technology’s Estuarine and Coastal Ocean Model (sECOM). Simulations with the ILP formulations run always faster than the ones with the base decomposition, and the complete formulation is better than the feasible one when it is applicable. In addition, parallel efficiency with the ILP decomposition may be greater than one.     

First arrival time picking for microseismic data based on DWSW algorithm

The first arrival time picking is a crucial step in microseismic data processing. When the signal-to-noise ratio (SNR) is low, however, it is difficult to get the first arrival time accurately with traditional methods. In this paper, we propose the double-sliding-window SW (DWSW) method based on the Shapiro-Wilk (SW) test. The DWSW method is used to detect the first arrival time by making full use of the differences between background noise and effective signals in the statistical properties. Specifically speaking, we obtain the moment corresponding to the maximum as the first arrival time of microseismic data when the statistic of our method reaches its maximum. Hence, in our method, there is no need to select the threshold, which makes the algorithm more facile when the SNR of microseismic data is low. To verify the reliability of the proposed method, a series of experiments is performed on both synthetic and field microseismic data. Our method is compared with the traditional short-time and long-time average (STA/LTA) method, the Akaike information criterion, and the kurtosis method. Analysis results indicate that the accuracy rate of the proposed method is superior to that of the other three methods when the SNR is as low as ??10 dB.     

Methods and algorithms for reconstructing three-dimensional distributions of electric conductivity and polarization in the medium by finite-element 3D modeling using the data of electromagnetic sounding

The methods and algorithms for automated three-dimensional (3D) inversion of the induced polarization and time-domain electromagnetic data by 3D finite-element modeling are considered. The work of the algorithms of automated 3D inversion is illustrated by the examples based on the synthetic data provided by 3D modeling.     

Assimilation of SLA and SST data into an OGCM for the Indian Ocean

 Remotely sensed observations of sea-level anomaly and sea-surface temperature have been assimilated into an implementation of the Miami Isopycnic Coordinate Ocean Model (MICOM) for the Indian Ocean using the Ensemble Kalman Filter (EnKF). The system has been applied in a hindcast validation experiment to examine the properties of the assimilation scheme when used with a full ocean general circulation model and real observations. This work is considered as a first step towards an operational ocean monitoring and forecasting system for the Indian Ocean. The assimilation of real data has demonstrated that the sequential EnKF can efficiently control the model evolution in time. The use of data assimilation requires a significant amount of additional processing and computational resources. However, we have tried to justify the cost of using a sophisticated assimilation scheme by demonstrating strong regional and temporal dependencies of the covariance statistics, which include highly anisotropic and flow-dependent correlation functions. In particular, we observed a marked difference between error statistics in the equatorial region and at off-equatorial latitudes. We have also demonstrated how the assimilation of SLA and SST improves the model fields with respect to real observations. Independent in situ temperature profiles have been used to examine the impact of assimilating the remotely sensed observations. These intercomparisons have shown that the model temperature and salinity fields better resemble in situ observations in the assimilation experiment than in a model free-run case. On the other hand, it is also expected that assimilation of in situ profiles is needed to properly control the deep ocean circulation. Received: 8 January 2002 / Accepted: 8 April 2002     

An improved grey wolf optimizer algorithm for the inversion of geoelectrical data

The grey wolf optimizer (GWO) is a novel bionics algorithm inspired by the social rank and prey-seeking behaviors of grey wolves. The GWO algorithm is easy to implement because of its basic concept, simple formula, and small number of parameters. This paper develops a GWO algorithm with a nonlinear convergence factor and an adaptive location updating strategy and applies this improved grey wolf optimizer (improved grey wolf optimizer, IGWO) algorithm to geophysical inversion problems using magnetotelluric (MT), DC resistivity and induced polarization (IP) methods. Numerical tests in MATLAB 2010b for the forward modeling data and the observed data show that the IGWO algorithm can find the global minimum and rarely sinks to the local minima. For further study, inverted results using the IGWO are contrasted with particle swarm optimization (PSO) and the simulated annealing (SA) algorithm. The outcomes of the comparison reveal that the IGWO and PSO similarly perform better in counterpoising exploration and exploitation with a given number of iterations than the SA.     

基于二次场方法的并行三维大地电磁正反演研究

快速且高精度的三维大地电磁法正反演是目前研究的热点.由于大地电磁法场源的平面波特性,以往的正演方法大多采用直接求解总场的方法,在边界强加二维边界条件.本文提出了一种基于二次场方法的三维大地电磁法正演算法,将平面波在层状背景模型中的响应作为场源项,得到二次场满足的偏微分方程,并利用交错网格有限差分法求取二次场.与其他学者的基于总场方法的结果的对比证明了本文采用方法的正确性.在基于二次场的正演算法基础上,实现了基于L-BFGS的三维反演方法,并对公开的数据集进行了反演.另外,针对大地电磁法的多频率观测特性,采用了基于MPI的分频并行策略对程序进行并行化,可达到接近线性的加速比.

     

基于BFGS法融合InSAR和GPS技术监测地表三维形变

虽然InSAR技术具有高精度、大范围和高空间分辨率等优点,但只能监测雷达视线方向上的一维地表形变;而GPS技术虽可以监测地表的三维形变,但其空间分辨率很低.本文针对融合InSAR和GPS技术监测地表高空间分辨率三维形变展开研究.首先证明了简单的局部最优化迭代算法就能求得综合InSAR和GPS监测地表形变速率的能量函数模型的全局最优估值.随后提出了利用BFGS局部最优算法反演最优的地表三维形变速率.该方法既能避免全局最优化算法计算复杂且难以收敛的问题,又能克服传统的解析法中数值计算不稳定的缺点.最后,通过模拟实验和美国南加州真实数据实验表明,该方法能够得到高精度的地表三维形变速率场.而且当观测或插值误差导致解析法误差较大时,BFGS方法仍能得到高精度、稳定的全局最优解.     

An equivalent layer magnetization model for Australia based on Magsat data

An equivalent layer magnetization model for Australia and adjacent oceanic areas is presented. The model is obtained by linear inversion of Magsat anomaly data measured in the altitude range 325–550 km. The anomaly data set has been isolated from the raw data set by use of models of the core field and very long wavelength external fields, and is internally consistent. Certain major structural features of the Australian continent are geographically associated with magnetization anomalies. A first-order difference is seen between the Tasman Zone and the Precambrian cratonic areas: magnetization anomalies are much more subdued in the former, possibly reflecting a shallowing of the Curie isotherm within the crust. A profile of the vertical integral of magnetization is presented for a crustal section extending from the Gawler Block to the southeast coast. It is shown that the magnetization variations are probably due partly, but not wholly, to depth to Curie isotherm variations; gross magnetization variations among at least three distinct crustal units must be involved.     

Variations of electron density and temperature in ionosphere based on the DEMETER ISL data

Observations of the Langmuir Probe Instrument(ISL,Instrument Sonde de Langmuir) onboard the DEMETER satellite during four years from 2006 to 2009 were used to analyze the tempo-spatial variations of electron density(Ne) and temperature(Te) in the ionosphere.Twenty four research bins with each covering an area with 10° in longitude and 2° in latitude were selected to study the spatial distributions of Ne and Te.The results indicate that both Ne and Te have strong annual variations in the topside ionosphere at 660 km altitude.The semiannual anomaly and equinoctial asymmetry which are usually well known as the features of F-layer also exist in the topside ionosphere at low-and mid-latitudes.The yearly variation of Ne is opposite to the peak electron density of the F2-layer(NmF2) at higher latitudes in daytime and both are similar in nighttime.Also the yearly variations of Te at low-latitude are contrary to that at 600 km in daytime and similar in nighttime.An interesting feature of nighttime Te at low-latitude is an obvious annual variation in the northern hemisphere and semiannual variation in the southern hemisphere.The yearly variations of Te in daytime have negative and positive correlation with Ne at mid-and high-latitudes,respectively.Both Ne and Te in the neighborhood bins at the same latitude have a high correlation.In ionospheric events analyzing,this information may help to understand the characteristics of the variation and to distinguish the reliable abnormality from the normal background map.     

基于自适应遗传算法的MRS-TEM联合反演方法研究

地面磁共振法(MRS)因具有定性、定量分析地下水能力,而备受关注.传统磁共振地层含水量反演多采用均匀半空间模型,忽略电阻率分布信息对结果的影响.针对这一问题,本文基于多层电介质中磁共振响应理论,提出MRS与瞬变电磁(TEM)联合反演方法,通过电阻率分布信息对含水量反演过程的实时修正,提高了解释结果的准确度.反演算法采用自适应遗传算法(AGA)进行,基于繁殖规则,动态调整交叉概率和变异概率,解决了标准遗传算法易未成熟收敛而难以得到全局最优解问题.模型数据表明,含噪10%情况下,联合反演仍能较准确地反映地下含水单元模型结构,对比MRS单独反演优势明显.同时,内蒙古白旗野外观测数据联合反演结果与钻井资料基本一致,充分验证了AGA反演算法的实用性及MRS-TEM联合反演的实际意义.     

混凝土结构缺陷检测的探地雷达资料波场反演方法

本文针对工程建筑中的混凝土无损探测问题,提出了一种新的反演成像方法.它避免了传统方法观察偏移探测数据的情况.这种方法以探地雷达为探测手段,以电磁场原理为理论依据,其结果更具定量化和可视性.具体算法是将阻尼高斯牛顿法和GCV方法(Generalized Cross Validation)结合起来,其中阻尼牛顿法可以避免传统牛顿法计算Hessian矩阵的困难,GCV方法可以自适应地选择正则化参数,既减小了计算量又克服了噪声及不适定性的影响.数值算例验证了算法的有效性.     

An improved TES algorithm based on the corrected ALPHA difference spectrum

Different from visible signals, thermal infrared radiances depend on both temperature and emissivity. It is a key problem for us to separate temperature and emissivity in thermal infrared remote sensing re- search. Another difficulty encountered in the retrieval of surface temperature is the correction of downwelling sky irradiance, because it is closely related to surface emissivity. When emissivity is un- known, the downwelling sky irradiance is difficult to be removed. In this paper, we introduce a correc- tion term of downwelling sky irradiance developed by Li and Becker into Wien’s approximation, to de- rive an improved ALPHA difference spectrum which is independent of temperature, and furthermore develop a correction term to remove the error of Wien’s approximation. Under the support of the above work, attractive features of Alpha derived emissivity method and ASTER TES algorithm are combined together to acquire a new Improved TES algorithm based on Corrected ALPHA Difference Spectrum (ICADS TES). Because a multi-band inversion technique is applied, and the operations of band ratios and differences are included in the algorithm, it can partly remove the influence of atmosphere and noises. Numerical simulation experiments show that for various combinations of atmosphere, land covers and surface temperatures, the algorithm is applicable and stable. Its accuracy for temperature is 0―1.5 K, and that for emissivity is 0―0.015. Compared with current TES algorithms, our method has clear physical meaning, is easy to be implemented, and is applicable for a wide temperature range and surface types. The results are not influenced by the directional characteristic of emissivity. Because ICADS TES does not need the support of a priori information of surface types, it is also not influenced by the accuracy of classification and the problem of mixture pixels. Compared with our former TES algorithm based on corrected Alpha difference spectra (CADS TES), the new algorithm takes the effect of downwelling atmospheric radiation into account. When the quantity of atmosphere radiation can be estimated precisely, the performance of ICADS TES is much better.     

An automated algorithm for river ice monitoring over the Susquehanna River using the MODIS data

Reliable and prompt information on river ice condition and extent is needed to make accurate hydrological forecasts to predict ice jams breakups and issue timely flood warnings. This study presents a technique to detect and monitor river ice using observations from the MODIS instrument onboard the Terra satellite. The technique incorporates a threshold‐based decision tree image classification algorithm to process MODIS data and to determine the extent of ice. To differentiate between ice‐covered and ice‐free pixels within the riverbed, the algorithm combines observations in the visible and near‐infrared spectral bands. The developed technique presents the core of the MODIS‐based river ice mapping system, which has been developed to support National Oceanic and Atmospheric Administration NWS's operations. The system has been tested over the Susquehanna River in northeastern USA, where ice jam events leading to spring floods are a frequent occurrence. The automated algorithm generates three products: daily ice maps, weekly composite ice maps and running cloud‐free composite ice maps. The performance of the system was evaluated over nine winter seasons. The analysis of the derived products has revealed their good agreement with the aerial photography and with in situ observations‐based ice charts. The probability of ice detection determined from the comparison of the product with the high‐resolution Landsat imagery was equal to 91%. A consistent inverse relationship was found between the river discharge and the ice extent. The correlation between the discharge and the ice extent as determined from the weekly composite product reached 0.75. The developed CREST River Ice Observation System has been implemented at National Oceanic and Atmospheric Administration–Cooperative Remote Sensing Science and Technology Center as an operational Web tool allowing end users and forecasters to assess ice conditions on the river. Copyright © 2012 John Wiley & Sons, Ltd.