基于Laplace方程的垂线偏差法反演全球海域重力异常

海洋重力场模型反演的质量主要依赖于采用测高数据的精度、空间分辨率和数据分布密集程度。本文联合Geosat GM/ERM、ERS-1 GM/ERM、TOPEX/Poseidon、Envisat、Cryosat-2、Jason-1 ERM/GM和SARAL/AltiKa等多种测高观测数据集,深入比较了多种波形重跟踪算法的效果,回波数据重跟踪处理后的沿轨海面高标准差。统计表明,Sandwell算法优于MLE-4算法、Davis阈值法、改进阈值法和β参数拟合法;基于不同测高数据波形重采样的结果给出了沿轨海面梯度计算中低通滤波的参数选择方法,并采用Sandwell提出的垂线偏差法,反演了全球海域1′×1′的重力场模型。检核表明,反演结果与DTU13和SIO V23.1模型检核的差值均方根分别为3.4、1.8 mGal,与NGDC船测数据的检核精度为4~8 mGal,且本文模型在部分典型海区内精度更优。     

一种多前缘多阈值的波形重构算法

为了改善近海测高精度,本文对近海测高波形的特征进行了研究,并在现有算法的基础上提出基于波前缘识别算法的多前缘多阈值波形重构方法.然后将该方法用于中国近海的卫星测高波形数据的处理,得到精度与分辨率更好的海面高数据.作为验证,将此海面高数据换算为重力异常并与船测重力进行比较,结果显示该方法可以明显提高卫星测高精度.     

Identification of tsunami deposits considering the tsunami waveform: An example of subaqueous tsunami deposits in Holocene shallow bay on southern Boso Peninsula, Central Japan

This study proposes a tsunami depositional model based on observations of emerged Holocene tsunami deposits in outcrops located in eastern Japan. The model is also applicable to the identification of other deposits, such as those laid down by storms. The tsunami deposits described were formed in a small bay of 10–20-m water depth, and are mainly composed of sand and gravel. They show various sedimentary structures, including hummocky cross-stratification (HCS) and inverse and normal grading. Although, individually, the sedimentary structures are similar to those commonly found in storm deposits, the combination of vertical stacking in the tsunami deposits makes a unique pattern. This vertical stacking of internal structures is due to the waveform of the source tsunamis, reflecting: 1) extremely long wavelengths and wave period, and 2) temporal changes of wave sizes from the beginning to end of the tsunamis.

The tsunami deposits display many sub-layers with scoured and graded structures. Each sub-layer, especially in sandy facies, is characterized by HCS and inverse and normal grading that are the result of deposition from prolonged high-energy sediment flows. The vertical stack of sub-layers shows incremental deposition from the repeated sediment flows. Mud drapes cover the sub-layers and indicate the existence of flow-velocity stagnant stages between each sediment flow. Current reversals within the sub-layers indicate the repeated occurrence of the up- and return-flows.

The tsunami deposits are vertically divided into four depositional units, Tna to Tnd in ascending order, reflecting the temporal change of wave sizes in the tsunami wave trains. Unit Tna is relatively fine-grained and indicative of small tsunami waves during the early stage of the tsunami. Unit Tnb is a protruding coarse-grained and thickest-stratified division and is the result of a relatively large wave group during the middle stage of the tsunami. Unit Tnc is a fine alternation of thin sand sheets and mud drapes, deposited from waning waves during the later stage of the tsunami. Unit Tnd is deposited during the final stage of the tsunami and is composed mainly of suspension fallout. Cyclic build up of these sub-layers and depositional units cannot be explained by storm waves with short wave periods of several to ten seconds common in small bays.     

CryoSat-2卫星海冰区域波形识别及海冰干舷高确定

利用40%阈值法对CryoSat-2卫星波形数据进行重跟踪,将波形特征参数和海冰浓度相结合,对海冰和Lead(浮冰之间的开阔水域)进行有效识别。利用沿轨前后搜索算法计算海冰干舷高,并引用AWI结果,绘制2011~2013年北冰洋多年冰区域和一年冰区域平均海冰干舷高变化趋势图。比较本文结果与AWI结果的各年同期数据,验证本文结果的可靠性。     

小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

基于波形拟合的中国东海地区410 km间断面附近速度结构研究

410 km间断面是地幔转换带的顶界面,对其速度结构和起伏形态开展地震学探测有助于认识地球内部物质组成和相关的地球动力学过程.本文选取了由中国数字地震台网记录到的位于琉球俯冲区的一个中源地震P波宽频带波形资料,利用三重震相波形拟合研究了中国东海地区410 km间断面附近的精细速度结构.结果表明：中国东海地区下方410 km间断面整体表现为一尖锐的速度界面且有8~15 km的小幅抬升;该间断面之上存在52~62 km厚的低速层,其P波速度降低0.5%~1.6%;440 km深度以下存在1.0%~3.0%的P波高速异常.结合前人在该地区的层析成像结果,我们推测该高速异常体与西太平洋俯冲板片在中国东海地区地幔转换带内的滞留有关;板片内水相E分解使得转换带内水含量增加,这引发了410 km间断面的抬升;410 km间断面之上的低速层应与含水矿物脱水导致的地幔橄榄岩部分熔融有关.     

考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟

2008年5月12日四川汶川地区发生M_W7.9地震,震中位置103.4°E,31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡,且因为高山等地形复杂区域抢险救灾的艰巨性,为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响,建立了包含地形起伏影响及去除地形影响的两类模型.同时,依据震源破裂过程运动学反演结果,建立了包含障碍体破裂过程的震源滑动模型,实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法,通过对地震波传播过程的数值计算和后处理分析,模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用,距离断层破裂带越近,地形效应越明显,反之,距离越远,则地形效应越微弱;（2）因为地形高差与障碍体的影响,地震造成的峰值可能出现在震中区域之外;（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近;而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段,处清平、安县附近;对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义.     

Improving gravity anomalies over China marginal sea from retracked Geosat and ERS-1 data

The quality of altimeter data and ocean tide model is critical to the recovery of coastal gravity anomalies. In this contribution, three retracking methods (threshold, improved threshold and Beta-5) are investigated with the aim of improving the altimeter data over a shallow water area. Comparison indicates that the improved threshold is the best retracking method over China Sea. Two ocean tide models, NAO99b and CSR4.0, are analyzed. Results show that different tide models used in the processing of altimeter data may result in differences more than 10 mGal in recovered coastal gravity anomalies. Also, NAO99b is more suitable than CSR4.0 over the shallow water area of China Sea. Finally, gravity anomalies over China Sea are calculated from retracked Geosat/GM and ERS-1/GM data by least squares collocation. Comparison with shipborne gravimetry data demonstrates that gravity anomalies from retracked data are significantly superior to those from non-retracked data. Our results have the same order as the other two altimeter-derived gravity models: Sandwell&Smith(V16) and DNSC08.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

2011年11月6日美国俄克拉何马5.6级地震震源参数测定及成因初步研究

2011年11月6日3时,美国俄克拉何马州发生了5.6级地震。此次地震是该地区100多年以来的最大地震,因该区地震活动性低,故需要对此次地震成因开展研究。本文采用CAP(Cut And Paste)方法对该地震分别进行近震、远震以及近远震联合反演,得到了震源机制解与震源深度。其中联合反演得到的机制解参数为:节面Ⅰ:234°/81°/-170°;节面Ⅱ:143°/81°/-9°;震级为MW=5.6。为了验证近震结果的可靠性,采用不同模型对近震数据进行反演测试,发现震源断层面解结果有3°的偏差,震源深度结果有0.5km偏差。CAP结果显示,该地震的质心深度为5.0km左右,属于浅源地震。同时,根据地震震源参数的标度分析,得到破裂区宽度约为7.0km,此次地震破裂区上缘可达1.4km深,下缘深度达8.4km。破裂上缘在该区域井的注水深度范围之内,由此推测该地震可能与当地页岩气开采过程中的注水活动有关。