岩石复电阻率Dias模型及其参数求取方法

Dias模型是众多描述岩石复电阻率的频散特性模型中的一种,它不但参数少,且可由这些参数导出具有特定物理意义的参数。针对Dias模型参数的求取及反演解的不确定性评价,这里提出了一种模拟退火阻尼最小二乘联合反演方法。该方法首先将参数求解空间离散化,利用模拟退火找到全局最优解区间,然后利用阻尼最小二乘法在解空间进行求解。该方法有效地克服了原来两种方法的反演速度慢、反演结果依赖于初值的缺点,并实现了对反演结果多解性的评价。研究结果表明,对于岩石复电阻率数据的处理解释,Dias模型的使用比传统的Cole模型及其变形更具有潜在的优势。     

不等电极距情况下高密度电阻率法的反演结果分析

通过高密度电阻率法的模型实验,得出不等电极距情况对直立高阻薄板和低阻球体反演结果的影响:缺失电极对它正下方的ρ_s影响最大,而异常形态主要在水平方向有明显的变化,异常体随电极缺失的位置有趋势变化,电极缺失的数目越多,异常形态变化越大,但是从整体来看反演剖面图基本相似。     

考虑关断时间效应的瞬变电磁一维反演

为克服因关断时间校正而带来的误差,尝试直接对包含关断时间效应的瞬变电磁资料进行反演。正演采用将斜阶跃波离散为多个阶跃波之和的方法,反演采用经典的马奎特法。对模型试算的结果表明,采用该方法直接对斜阶跃响应进行反演,结果比较理想,计算速度也可以接受。     

水气两相流系统K-S-P模型参数反演的最优估计

应用渗流力学和参数反演辨识最优化理论,建立了水气两相流控制方程,给出了控制两相流系统K-S-P模型参数之间的本构关系,并对其中的相关参数进行数值反演辨识,得出待求参数的最优估计值。对不同质地土壤的实验测试值进行对比分析,其结果吻合较好,验证了模型的适定性和可靠性,解决了实验数据处理以及结果对初值的敏感性问题,并且该数值方法可通过多组初值参数进行数值实验,结合实验结果来满足实际工程的需求。这不仅对非饱和带内污染物传输、气藏工程试井计算、填埋气体的释放以及挥发性有机污染物的去污分析中两相流数值计算提供基础,而且对K-S-P本构关系之间参数的确定提供了可靠的理论依据。     

Electrical aerosol spectrometer of Tartu University

The electrical aerosol spectrometer (EAS) of the parallel measuring principle at Tartu University is an efficient instrument for rapid measurement of the unstable size spectrum of aerosol particles. The measuring range from 10 nm to 10 μm is achieved by simultaneously using a pair of differential mobility analyzers with two different particle chargers. The particle spectrum is calculated and measurement errors are estimated in real time by using a least-squares method. Experimental calibration ensures reliability of measurement. The instrument is well suited for continuous monitoring of atmospheric aerosol.     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Structure and tectonic evolution of the Southern Eurasia Basin, Arctic Ocean

Multichannel seismic reflection data acquired by Marine Arctic Geological Expedition (MAGE) of Murmansk, Russia in 1990 provide the first view of the geological structure of the Arctic region between 77–80°N and 115–133°E, where the Eurasia Basin of the Arctic Ocean adjoins the passive-transform continental margin of the Laptev Sea. South of 80°N, the oceanic basement of the Eurasia Basin and continental basement of the Laptev Sea outer margin are covered by 1.5 to 8 km of sediments. Two structural sequences are distinguished in the sedimentary cover within the Laptev Sea outer margin and at the continent/ocean crust transition: the lower rift sequence, including mostly Upper Cretaceous to Lower Paleocene deposits, and the upper post-rift sequence, consisting of Cenozoic sediments. In the adjoining Eurasia Basin of the Arctic Ocean, the Cenozoic post-rift sequence consists of a few sedimentary successions deposited by several submarine fans. Based on the multichannel seismic reflection data, the structural pattern was determined and an isopach map of the sedimentary cover and tectonic zoning map were constructed. A location of the continent/ocean crust transition is tentatively defined. A buried continuation of the mid-ocean Gakkel Ridge is also detected. This study suggests that south of 78.5°N there was the cessation in the tectonic activity of the Gakkel Ridge Rift from 33–30 until 3–1 Ma and there was no sea-floor spreading in the southernmost part of the Eurasia Basin during the last 30–33 m.y. South of 78.5°N all oceanic crust of the Eurasia Basin near the continental margin of the Laptev Sea was formed from 56 to 33–30 Ma.     

Interpretation of seismic and potential field data from western New York State and Lake Ontario

Lithoprobe and industry seismic profiles have furnished evidence of major zones of easterly dipping Grenville deformed crust extending southwest from exposed Grenville rocks north of Lake Ontario. Additional constraints on subsurface structure limited to the postulated Clarendon–Linden fault system south of Lake Ontario are provided by five east–west reflection lines recorded in 1976. Spatial correlations between seismic structure and magnetic anomalies are described from both Lake Ontario and the newly reprocessed New York lines.In the Paleozoic to Precambrian upper crust, the New York seismic sections show: (1) An easterly thickening wedge of subhorizontal Paleozoic strata unconformably overlying a Precambrian basement whose surface has an apparent regional easterly dip of 1–2°. Minor apparent normal offsets, possibly on the order of tens of meters, occur within the Paleozoic section. The generally poorly reflective unconformity may be locally characterized by topographic relief on the order of 100 m; (2) Apparent local displacement on the order of 90 m at the level of the Black River Group diminishes upward to little or no apparent offset of Queenston Shale; (3) Within the limited seismic sections, there appears to be no evidence that the complete upper crustal section is vertically or subvertically offset; (4) Dipping structure in the Paleozoic strata (15° to 35°) resembles some underlying Precambrian basement elements; (5) The surface continuity of inferred faults constituting the Clarendon–Linden system is not strongly supported by the seismic data.Beneath the Paleozoic strata, the seismic sections show both linear and arcuate reflector geometry with easterly apparent dips of 15° to 35° similar to the deep structures imaged on seismic lines from nearby Lake Ontario and on Lithoprobe lines to the north. The similarity supports an extension of easterly dipping Central Metasedimentary Belt structures of the Grenville orogen from southern Ontario to beneath western New York State.From a comparison of the magnetic and gravity fields with the New York seismic sections, we suggest: (1) The largely nonmagnetic Paleozoic strata appear to contribute negligibly to magnetic anomalies. Seismically imaged fractures in the New York Paleozoic strata appear to lie mainly west of a positive gravity anomaly. The relationship between magnetic and gravity anomalies and the changes in the geometry of interpreted Precambrian structures remains enigmatic; (2) North to northeast trending curvilinear magnetic and gravity anomalies parallel, but are not restricted to the principal trend of the postulated Clarendon–Linden fault system. Paleozoic fractures of the Clarendon–Linden system may partly overlie a southward extension of the Composite Arc Belt boundary zone.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.