GPS授时地震仪走时误差校正

受制造成本和工艺的限制,GPS授时地震仪中的守时晶振不可避免地会有频率误差,继而在数据采集过程中,随着采集时间的增加,各台站之间会出现采集不同步现象。介绍了Akima分段多项式插值重采样方法及其在GPS授时地震仪时间校正中的应用。通过对理论数据与实际数据的Akima插值重采样试验表明,Akima插值重采样方法可应用于GPS授时地震仪走时误差校正,是一种简便、高效、保形性好的插值重采样方法。      

Occurrence of persistent organic pollutants in marine fish from the Natuna Island,South China Sea

Five marine fish species were collected from the Natuna Island, South China Sea to investigate the occurrence of polybrominated diphenyl ethers (PBDEs), polychlorinated biphenyls (PCBs), dichlorodiphenyltrichloroethane and its metabolites (DDTs). Concentrations of PBDEs, PCBs, and DDTs in marine fish ranged from 2.85 to 7.82, 14.3 to 48.1, and 7.99 to 40.3 ng/g lipid weight, respectively. Higher concentrations of PBDEs, PCBs, and DDTs were observed in Snakefish (Trachinocephalus myops), which might be attributed to their different feeding and living habits. PCBs were the predominant POPs in all marine fish, followed by DDTs and PBDEs. BDE 47 and PCB 153 were the predominant congener of PBDEs and PCBs, respectively. Compositional distribution of DDTs indicated the possible presence of fresh input sources around the Natuna Island. The ratios of o,p′-DDT/p,p′-DDT being less than 1 in fish samples suggested that DDT contributions from dicofol seemed considerably low. New input sources of DDT in South China Sea are worth further research.     

Regional gravity field modeling using free-positioned point masses

A two-step free-positioned point mass method is used for regional gravity field modeling together with the remove-compute-restore (RCR) technique. The Quasi-Newton algorithm (L-BFGS-B) is implemented to solve the nonlinear problem with bound constraints in the first step, while in the second step the magnitudes of the point masses are re-adjusted with known positions in the least-squares sense. In order to reach a good representation of the gravity field, a number of parameter sets have to be defined carefully before the computations. The effects of four important parameter sets (depth limits, number of point masses, original/reduced basis functions and optimization directions) are investigated for regional gravity field modeling based on two numerical test cases with synthetic and real data. The results show that the selection of the initial depth and depth limits is of most importance. The number of point masses for obtaining a good fit is affected by the data distribution, while a dependency on the data variability (signal variation) is negligible. Long-wavelength errors in the predicted height anomalies can be reduced significantly by using reduced basis functions, and the radial-direction optimization proves to be stable and reliable for regular and irregular data scenarios. If the parameter sets are defined properly, the solutions are similar to the ones computed by least-squares collocation (LSC), but require fewer unknowns than LSC.     

卫星高光谱大气CO2遥感反演精度地基验证研究

利用TCCON网站提供的北半球7个地面观测站CO₂干空气混合比（XCO₂）数据,对3种卫星反演的XCO₂产品进行了验证,包括SCIAMACHY产品、NIES-GOSAT产品和ACOS-GOSAT产品。结果表明：卫星CO₂遥感反演产品与地基遥感资料具有较一致的季节性周期变化,一年中月平均浓度最高值均出现在4月和5月,最低值均出现在8月和9月;相对于地面观测,3种卫星产品均低估了XCO₂;ACOS-GOSAT产品与NIES-GOSAT产品的精度大体相当,误差标准差分别为2.26×10^-6和2.27×10^-6;SCIAMACHY产品的精度略差,误差标准差为2.91×10^-6。     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.     

间断有限元方法的数值频散分析及其波场模拟

数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题,从理论推导和数值分析的角度进行了深入研究,并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method,简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid,简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法,特别是空间精度为3阶的RKDG方法,即使当空间步长取波长的一半,即一个波长内取2个网格点时,最大的频散误差也不超过1.67%.同时,我们也通过波场模拟对比研究了不同数值方法的数值频散问题,进一步直观地验证了数值频散的理论分析结果.     

伊洛瓦底盆地热-沉降史模拟及构造-热演化特征

本文首先运用EASY% Ro反演法对伊洛瓦底盆地由北向南进行了热史的恢复,北部钦敦凹陷的平均古地温梯度为13.0~15.0 ℃/km,中部沙林凹陷的平均古地温梯度为18.0~22.0 ℃/km,南部三角洲凹陷的平均古地温梯度为33.0~37.0 ℃/km.从模拟结果可以看出,盆地由北向南地温梯度逐渐升高,生烃门限的深度由深变浅.然后模拟了盆地的构造沉降史.模拟结果表明,盆地具有幕式构造沉降特征,这反映了伊洛瓦底盆地可能处于弧间或弧后的构造背景.伊洛瓦底盆地北部和南部具有不同幕次的构造沉降史,北部在早始新世时期（53~51 Ma）经历了一幕拉伸过程,然后进入了热沉降期,并伴随局部的快速隆升;南部则经历了两幕拉张过程,分别是在早始新世时期（53~51 Ma）和中中新世时期（21~13 Ma）.盆地的这种南北构造沉降的差异很可能是造成盆地地温梯度北低南高的原因.     

黑龙江及邻区地壳应力场初步研究

采用综合震源机制方法,利用5 002条P波极性数据反演黑龙江地区二维应力结构,求得该区0.5°×0.5°综合震源机制解.初步研究表明：21世纪以来,构造应力场总特征基本一致,具有明显的继承性,主压应力和主张应力的优势方向明显,与板块运动关系具有成因上的联系.     

图像ROI选择及其应用研究

在辐射成像应用中,DR及CT检测会产生大量的图像数据,对这些图像数据的处理涉及图像处理技术的方方面面。这些处理技术包括重建、可视化、图像增强、图像压缩和图像缺陷分析等。通常人们需要对图像中的特定区域（ROI）进行处理,并对图像上的ROI进行选择。ROI的确定,编辑是辐射成像图像处理的常用功能。本文介绍常用各种ROI,实现了这些ROI的勾画、编辑、选择等功能,并介绍了ROI与孔隙度计算及缺陷识别等应用的结合。     

Rupture pattern of the Oct 23, 2011 Van-Merkez,Eastern Turkey earthquake

High-frequency rupture process of the Oct 23, 2011 Van-Merkez earthquake is imaged by back-projection method using high-quality teleseismic P wave data from the US Array, and prestack Kirchhoff migration using P wave data from a subarray of global seismic networks. The rupture model with two asperities is confirmed by previous two methods. In low-frequency imaging, a large asperity derived from the migration method corresponds to the second one from the high-frequency P waves. The consistency of the locations of asperities from datasets with different frequency bands indicates that there is possible insignificance of the frequency-dependent feature for the earthquake. The resultant images illustrate the spatial and temporal evolution of the rupture, which mainly propagated WSW over a length of 33 km during the first 18 s, accompanying with bursts of two asperities at 3 and 11–13 s. The rupture direction is confirmed by the S wave corner frequency variations of strong ground accelerations. The rupture fronts are mainly located at the updip of the causative fault. Based on polarities of the P waveforms and focal mechanisms of the mainshock and aftershocks, the failure of these two asperities is determined to have occurred on a reverse fault with a dip angle of 47°. Hence, the rupture pattern of the 2011 Van-Merkez earthquake was dominated by a unilateral rupture toward the west-southwest direction.