利用ICESat数据确定格陵兰冰盖高程和体积变化   总被引：1，自引：0，他引：1       下载免费PDF全文

陈国栋  张胜军 《地球物理学报》2019,62(7):2417-2428

两极冰盖消融是造成海平面上升的重要原因，作为世界第二大冰盖，格陵兰冰盖消融速度在进入21世纪以后明显加快，引起了广泛关注.本文利用ICESat卫星激光测高数据，探讨了坡度改正的方法，通过改进平差模型解决了病态问题，并采用重复轨道方法计算了2003年9月至2009年10月间格陵兰冰盖的体积和高程变化趋势，对格陵兰冰盖各冰川流域系统的变化情况进行了详细分析.结果表明，格陵兰冰盖在这6年间平均高程变化趋势为-16.79±0.84 cm·a^-1，体积变化速率为-301.37±15.16 km³·a^-1，体积流失主要发生在冰盖边缘，其中DS1、DS8等流域的体积损失正在加剧，而高程在2000 m以上的冰盖内陆地区表现出高程积聚的状态，但增长速度明显减缓.与现有研究成果的对比表明，算法优化后的本文结果更具可靠性.  相似文献   

一种卫星重力梯度数据的ARMA最优滤波模型构建方法          下载免费PDF全文

朱广彬  常晓涛  付兴科  瞿庆亮 《地球物理学报》2018,61(12):4729-4736

由于卫星重力梯度观测的有色噪声特性和海量观测特征，在利用直接法进行重力场模型的最小二乘求解时，观测值的协方差阵为超大型的非对角阵，这给数值求解带来了极大困难.本文提出了一种基于先验误差功率谱密度的最优ARMA滤波模型构建方法，结合法方程的分块求解策略，可实现对卫星重力梯度观测值的高效滤波处理.数值仿真结果表明，利用最优ARMA滤波器进行时域滤波后，法方程的态性得到了明显改善，重力梯度观测值中的有色噪声得到了有效的"白化"处理，大地水准面精度得到了显著提升.  相似文献   

Scattering pattern analysis and true-amplitude generalized Radon transform migration for acoustic transversely isotropic media with a vertical axis of symmetry     

Quan Liang  Wei Ouyang  Weijian Mao  Shijun Cheng 《Geophysical Prospecting》2020,68(4):1211-1227

In multi-parameter ray-based anisotropic migration/inversion, it is essential that we have an understanding of the scattering mechanism corresponding to parameter perturbations. Because the complex nonlinearity in the anisotropic inversion problem is intractable, the construction of true-amplitude linearized migration/inversion procedures is needed and important. By using the acoustic medium assumption for transversely isotropic media with a vertical axis of symmetry and representing the anisotropy with P-wave normal moveout velocity, Thomsen parameter δ and anelliptic parameter η, we formalize the linearized inverse scattering problem for three-dimensional pseudo-acoustic equations. Deploying the single-scattering approximation and an elliptically anisotropic background introduces a new linear integral operator that connects the discontinuous perturbation parameters with the multi-shot/multi-offset P-wave scattered data. We further apply the high-frequency asymptotic Green's function and its derivatives to the integral operator, and then the scattering pattern of each perturbation parameter can be explicitly presented. By naturally establishing a connection to generalized Radon transform, the pseudo-inverse of the integral operator can be solved by the generalized Radon transform inversion. In consideration of the structure of this pseudo-inverse operator, the computational implementation is done pointwise by shooting a fan of rays from the target imaging area towards the acquisition system. Results from two-dimensional numerical tests show amplitude-preserving images with high quality.  相似文献   

Separation of magnetic anomalies into induced and remanent magnetization contributions     

Jamaledin Baniamerian  Shuang Liu  Xiangyun Hu  Maurizio Fedi  Mahak Singh Chauhan  Mahmoud Ahmed Abbas 《Geophysical Prospecting》2020,68(7):2320-2342

Inversion of magnetic data is complicated by the presence of remanent magnetization, and it provides limited information about the magnetic source because of the insufficiency of data and constraint information. We propose a Fourier domain transformation allowing the separation of magnetic anomalies into the components caused by induced and remanent magnetizations. The approach is based on the hypothesis that each isolated source is homogeneous with a uniform and specific Koenigsberger ratio. The distributions of susceptibility and remanent magnetization are subsequently recovered from the separated anomalies. Anomaly components, susceptibility distribution and distribution of the remanent and total magnetization vectors (direction and intensity) can be achieved through the processing of the anomaly components. The proposed method therefore provides a procedure to test the hypotheses about target source and magnetic field, by verifying these models based on available information or a priori information from geology. We test our methods using synthetic and real data acquired over the Zhangfushan iron-ore deposit and the Yeshan polymetallic deposit in eastern China. All the tests yield favourable results and the obtained models are helpful for the geological interpretation.  相似文献   

Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation     

Raphael Di Carlo Silva dos Santos  João Carlos Ribeiro Cruz  Manuel de Jesus dos Santos Costa 《Geophysical Prospecting》2020,68(7):2046-2063

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.  相似文献   

Inter-comparison of remote sensing platforms for height estimation of mango and avocado tree crowns     

《International Journal of Applied Earth Observation and Geoinformation》2020

To support the adoption of precision agricultural practices in horticultural tree crops, prior research has investigated the relationship between crop vigour (height, canopy density, health) as measured by remote sensing technologies, to fruit quality, yield and pruning requirements. However, few studies have compared the accuracy of different remote sensing technologies for the estimation of tree height. In this study, we evaluated the accuracy, flexibility, aerial coverage and limitations of five techniques to measure the height of two types of horticultural tree crops, mango and avocado trees. Canopy height estimates from Terrestrial Laser Scanning (TLS) were used as a reference dataset against height estimates from Airborne Laser Scanning (ALS) data, WorldView-3 (WV-3) stereo imagery, Unmanned Aerial Vehicle (UAV) based RGB and multi-spectral imagery, and field measurements. Overall, imagery obtained from the UAV platform were found to provide tree height measurement comparable to that from the TLS (R² = 0.89, RMSE = 0.19 m and rRMSE = 5.37 % for mango trees; R² = 0.81, RMSE = 0.42 m and rRMSE = 4.75 % for avocado trees), although coverage area is limited to 1–10 km² due to battery life and line-of-sight flight regulations. The ALS data also achieved reasonable accuracy for both mango and avocado trees (R² = 0.67, RMSE = 0.24 m and rRMSE = 7.39 % for mango trees; R² = 0.63, RMSE = 0.43 m and rRMSE = 5.04 % for avocado trees), providing both optimal point density and flight altitude, and therefore offers an effective platform for large areas (10 km²–100 km²). However, cost and availability of ALS data is a consideration. WV-3 stereo imagery produced the lowest accuracies for both tree crops (R² = 0.50, RMSE = 0.84 m and rRMSE = 32.64 % for mango trees; R² = 0.45, RMSE = 0.74 m and rRMSE = 8.51 % for avocado trees) when compared to other remote sensing platforms, but may still present a viable option due to cost and commercial availability when large area coverage is required. This research provides industries and growers with valuable information on how to select the most appropriate approach and the optimal parameters for each remote sensing platform to assess canopy height for mango and avocado trees.  相似文献   

利用海洋重力场变化分析洋底板块运动          下载免费PDF全文

李倩倩  鲍李峰  C. K. Shum 《地球物理学报》2020,63(7):2506-2515

洋底板块运动是地球动力学和全球变化研究的重要内容.本文根据质量迁移与地球外部重力场变化的对应关系,利用不同时期测高资料推算的1995—2019全球海洋重力场变化结果,反演分析全球洋底板块运动特征.研究表明,板块汇聚边界、板块内无震海岭、海山群、断裂带等区域重力异常变化显著,而在板块离散边界无明显变化趋势;西南印度洋中脊、大西洋中脊、中印度洋中脊等地区重力异常垂直梯度变化显著,且在西太平洋俯冲带、部分海岭区域也存在明显变化,其空间分布与地形基本吻合.海洋重力场变化整体上准确反映了全球洋底板块构造运动.相较于重力异常变化反演结果,重力垂直梯度的变化能够更为准确地反映洋底板块运动特征,特别是在洋中脊区域,扩张速率越小,垂直重力梯度变化越显著.此外,详细讨论了测高海洋重力场不确定因素对洋底板块运动分析结果的影响,海面坡度改正是主要因素之一.  相似文献   

Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement     

João B.C. Silva  Darcicléa F. Santos  Daniele P. Monteiro 《Geophysical Prospecting》2020,68(3):1057-1071

We advance a principle directed to assigning numerical values to free parameters usually present in inversion methods. It may be formulated as: ‘Optimum estimates of free parameters in an inversion procedure must lead, in tests using synthetic data, to solutions whose geometrical expression reflects the main qualitative or semiquantitative geological characteristic of the study area.’ To this end, the interpreter should (i) specify a typical anomalous source geometry which incorporates the most relevant geological information for the study area, (ii) compute the corresponding gravity anomaly and (iii) invert the anomaly for the source geometry finding the numerical values of the free parameters that lead to a solution closest to the typical source. Application of the above methodology to synthetic and real data from the basement relief of a rift basin has asserted its efficacy.  相似文献   

Deep neural network application for 4D seismic inversion to changes in pressure and saturation: Optimizing the use of synthetic training datasets     

Gustavo Côrte  Jesper Dramsch  Hamed Amini  Colin MacBeth 《Geophysical Prospecting》2020,68(7):2164-2185

In this work, we tackle the challenge of quantitative estimation of reservoir dynamic property variations during a period of production, directly from four-dimensional seismic data in the amplitude domain. We employ a deep neural network to invert four-dimensional seismic amplitude maps to the simultaneous changes in pressure, water and gas saturations. The method is applied to a real field data case, where, as is common in such applications, the data measured at the wells are insufficient for properly training deep neural networks, thus, the network is trained on synthetic data. Training on synthetic data offers much freedom in designing a training dataset, therefore, it is important to understand the impact of the data distribution on the inversion results. To define the best way to construct a synthetic training dataset, we perform a study on four different approaches to populating the training set making remarks on data sizes, network generality and the impact of physics-based constraints. Using the results of a reservoir simulation model to populate our training datasets, we demonstrate the benefits of restricting training samples to fluid flow consistent combinations in the dynamic reservoir property domain. With this the network learns the physical correlations present in the training set, incorporating this information into the inference process, which allows it to make inferences on properties to which the seismic data are most uncertain. Additionally, we demonstrate the importance of applying regularization techniques such as adding noise to the synthetic data for training and show a possibility of estimating uncertainties in the inversion results by training multiple networks.  相似文献   

Inversion-based deblending using migration operators     

Amr Ibrahim  Daniel Trad 《Geophysical Prospecting》2020,68(8):2459-2470

In this paper, we compare the denoising- and inversion-based deblending methods using Stolt migration operators. We use Stolt operator as a kernel to efficiently compute apex-shifted hyperbolic Radon transform. Sparsity promoting transforms, such as Radon transform, can focus seismic data into a sparse model to separate signals, remove noise or interpolate missing traces. Therefore, Radon transforms are a suitable tool for either the denoising- or the inversion-based deblending methods. The denoising-based deblending treats blending interferences as random noise by sorting the data into new gathers, such as common receiver gather. In these gathers, blending interferences exhibit random structures due to the randomization of the source firing times. Alternatively, the inversion-based deblending treats blending interferences as a signal, and the transform models this signal by incorporating the blending operator to formulate an inversion problem. We compare both methods using a robust inversion algorithm with sparse regularization. Results of synthetic and field data examples show that the inversion-based deblending can produce more accurate signal separation for highly blended data.  相似文献