Tectonic core of a sedimentary drift: a potential ridge propagation feature beneath the Blake Outer Ridge

The Blake Outer Ridge is a 480–kilometer long linear sedimentary drift ridge striking perpendicular to the North American coastline. By modeling free-air gravity anomalies we tested for the presence of a crustal feature that may control the location and orientation of the Blake Outer Ridge. Most of our crustal density models that match observed gravity anomalies require an increase in oceanic crustal thickness of 1–3 km on the southwest side of the Blake Outer Ridge relative to the northeast side. Most of these models also require 1–4 km of crustal thinning in zone 20–30 km southwest of the crest of the Blake Outer Ridge. Although these features are consistent with the structure of oceanic fracture zones, the Blake Outer Ridge is not parallel to adjacent known fracture zones. Magnetic anomalies suggest that the ocean crust beneath this feature formed during a period of mid-ocean ridge reorganization, and that the Blake Outer Ridge may be built upon the bathymetric expression of an oblique extensional feature associated with ridge propagation. It is likely that the orientation of this trough acted as a catalyst for sediment deposition with the start of the Western Boundary Undercurrent in the mid-Oligocene.     

赤道西太平洋大气边界层的诊断分析

利用１９９３年ＥＮＳＯ事件爆发（４月）前酝酿时期“热带大洋与全球大气－海洋耦合响应试验”强化观测阶段“向阳红五号”科学考察船１５５°Ｅ，２°Ｓ定点海洋气象和高空大气探测资料，分析了赤道西太平洋大气边界层特征。结果表明：赤道西太平洋上空辐合对流区边界层内也有等温或逆温层存在，大气边界层物理参数变化与大尺度海－气变异有密切关系，赤道太平洋沃克环流加强、东移，边界层上部逆温层出现概率明显增大，边界层内高湿层湿度减小。另外，边界层内实测风速、风向廓线随高度的变化基本上不服从Ｅｋｍａｎ规律，但由其平均，分量合成得到的平均风随高度的变化却基本符合Ｅｋｍａｎ规律。     

Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data

A wide-angle seismic survey, combining ocean-bottom seismometers (OBS) and multi-channel seismic (MCS) profiling, was implemented in the southwestern Ryukyu subduction zone during August and September 1995. In this paper, we present the data analysis of eight OBSs and the corresponding MCS line along profile EW9509-1 from this experiment. Seismic data modeling includes identification of refracted and reflected arrivals, initial model building from velocity analysis of the MCS data, and simultaneous and layer-stripping inversions of the OBS and MCS arrivals. The velocity-interface structure constructed along profile EW9509-1 shows that the northward subduction of the Philippine Sea Plate has resulted in a northward thickening of the sediments of the Ryukyu Trench and the Yaeyama accretionary wedge north of the trench. The boundary between the subducting oceanic crust and the overriding continental crust (represented by a velocity contour of 6.75 km/s) and a sudden increase of the subducting angle (from 5 degrees to 25 degrees) are well imaged below the Nanao Basin. Furthermore, velocity undulation and interface variation are found within the upper crust of the Ryukyu Arc. Therefore, the strongest compression due to subduction and a break-off of the slab may have occurred and induced the high seismicity in the forearc region. This revised version was published online in November 2006 with corrections to the Cover Date.     

井中激电地-井方式井旁球体正反演

以激发极化法中的地-井方式中的体极化球形体作为研究对象,采用解析法对地下半空间的电位场进行求解,利用提取出的井中二次场异常电位差进行反演,得出极化球体的球心离井的距离和球心的埋深。通过计算机进行正反演模拟,验证了进行正反演计算的解析式的正确性,并从物理意义上对正反演过程中出现的现象进行分析。最后给出该方法的适用条件以及注意事项。研究对激发极化法地下勘探目标的定量分析具有明确意义。     

直流电测深法与中心回线瞬变电磁法联合反演

分析了直流电测深法(VES)与中心回线瞬变电磁法(TEM)对地层的响应特征,结果显示:TEM对低阻层反应灵敏,而VES对高阻层具有较高的分辨能力。基于TEM与VES对地层响应特征互补的特点,尝试将TEM与VES数据进行联合反演,模型试算结果表明,联合反演的效果明显比单一方法反演的效果好。     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

地面核磁共振模型约束反演含水层参数

导电模型的地面核磁共振感应电动势是含水量的非线性函数.引入模型约束的迭代反演方法求解该非线性问题的反演.在反演过程中,根据理论公式计算迭代过程中灵敏度矩阵,并采用平坦模型和光滑模型两种约束.对均匀半空间、层状导电模型和实际数据进行了反演模拟,结果表明,模型约束迭代反演方法能从地面核磁共振感应电动势获得含水层较为合理的含水量及分布,且该结果可以从作为初始模型的均匀含水量分布反演得到.对无噪音数据,平坦模型和光滑模型约束对反演的含水量分布影响不大;但当数据存在一定的噪音时,平坦模型约束将比光滑模型约束获得更为精确的含水层参数.     

小波边缘分析建模波阻抗反演方法在煤层解释中的应用

受限于地震波主波长的1／4,常规地震处理难以分辨较薄煤层。基于模型的波阻抗反演虽然克服了厚、薄煤层计算问题．但因难以得到准确的高频部分,很多情况下反演结果的精度和分辨率往往不能满足地质预测的要求。小波边缘分析建模波阻抗反演方法是从地震数据中提取地震特征参数,用于建模并参与迭代反演。该方法即可以弥补反演过程中高频成分的误差,减少对钻井数据的依赖程度,得到较为合理的初始模型,又可提高反演的精度和分辨率．使反演结果更好的反映实际地下地质情况。利用该方法,对大井-将军庙及新疆红沙泉2个勘探区的侏罗系中统西山窑组含煤地层进行波阻抗反演,反演结果与钻井资料进行比较,其中大井一将军庙0802及0809号孔单煤层厚度误差分别为3％、8％．新疆红沙泉多煤层厚度及煤层之间的夹层厚度预测结果除B6煤层误差较大外,其它误差均小于7％．可见该反演方法的计算结果比较准确。     

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!