THE COLLOCATION APPROACH TO THE INVERSION OF GRAVITY DATA1

Gravity data inversion can provide valuable information on the structure of the underlying distribution of mass. The solution of the inversion of gravity data is an ill-posed problem, and many methods have been proposed for solving it using various systematic techniques. The method proposed here is a new approach based on the collocation principle, derived from the Wiener filtering and prediction theory. The natural multiplicity of the solution of the inverse gravimetric problem can be overcome only by assuming a substantially simplified model, in this case a two-layer model, i.e. with one separation surface and one density contrast only. The presence of gravity disturbance and/or outliers in the upper layer is also taken into account. The basic idea of the method is to propagate the covariance structure of the depth function of the separation surface to the covariance structure of the gravity field measured on a reference plane. This can be done since the gravity field produced by the layers is a functional (linearized) of the depth. Furthermore, in this approach, it is possible to obtain the variance of the estimation error which indicates the precision of the computed solution. The method has proved to be effective on simulated data, fulfilling the a priori hypotheses. In real cases which display the required statistical homogeneity, good preliminary solutions, useful for a further quantitative interpretation, have also been derived. A case study is discussed.     

水平电偶源频率测深中全区视电阻率定义的新方法

从电磁场的统一性出发，分析了水平电偶源频率测深中等效电阻率的算法原理和具体定义．分析表明，等效电阻率可明显地改善波区视电阻率在非波区场的畸变，较为形象、直观地反映出地电断面的电性变化．数值计算、三维模型计算和实际数据处理表明，等效电阻率具有良好的特性和应用前景．     

压扭性断裂的震前应力场与地应力变化异常

为了探讨浅源构造地震孕育和发生的力学过程,本文提出了一种力学模式,其中把发震断层的震前状态划分为几个阶段,即相对稳定的闭锁区内应力积累——滑动区逐渐扩展、闭锁区收缩并造成闭锁应力高度集中——剩余闭锁区缩小至一定程度而导致断层剧烈错动的几个阶段.计算并分析了震源力学状态这种变化在其周围产生的附加应力场,以及反映这一过程的地应力变化特征.其结论与多年来地应力观测实践中取得的某些经验性认识相符,并较好地解释了1975年海城地震前的若干实测资料.      

拜城—大柴旦剖面的上地壳Q值结构

利用横跨塔里木盆地、阿尔金造山带和柴达木盆地的拜城—大柴旦综合地球物理剖面的人工地震宽角反射/折射资料，对拜城—大柴旦剖面研究区上地壳的Q值结构进行了反演. 结果表明，研究区上地壳的Q值结构具有明显的分层性，Q值随深度的加大而增大；横向对应的不同地质构造单元的Q值具有不同的分布特点. 塔里木盆地的上地壳具有稳定的Q值结构，但盆地北缘（特别是库车坳陷）的Q值比南缘的低，表明北缘的介质比较破碎. 这可能与盆地北缘活跃的构造特征与巨大的沉积厚度有关，是引起盆地南北两侧地震活动性差异（北缘强，南缘弱）的重要原因之一. 阿尔金造山带上地壳的平均Q值较高（约500），这可能由于结晶基底埋藏较浅，基岩出露，因此地震波在此处传播能量衰减较慢所致. 柴达木盆地西半部分上地壳的Q值除了表层的（厚1.0～2.0km）较高外，其余各层的Q值与塔里木盆地中部的相当（平均约400），显示了与塔里木盆地相似的稳定构造.     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

北京以西至晋冀蒙交界地区的应力场

本文根据北京以西至晋冀蒙交界地区的地震地质构造条件和地震活动资料,用有限单元法与位错理论分别研究了该地区的构造应力场与地震应力场,结合应力场的分布特征和近期小地震的活动特点,推断该地区存在两个应力相对集中区,即(1)黑山寺断裂以北、牛心川——四回庄断裂以东及梁家庄断裂以南所围成的地区;(2)天镇北断裂与韭菜沟——黄土堌堆断裂和张皋断裂汇交区所围成的长轴呈北西向的区域.其中,第一个应力相对集中区内怀来以西的下花园断裂、浑井东断裂和桑干河断裂的汇而不交处及其附近地区,未来发生强震的      

地壳形变预报地震途径的探索

本文用系统科学理论河北省及邻区形变台网为一预报系统，通过对实用化攻关ＰＤＢＳ软件系统等方面处理的结果进行全面分析，得到表征系统的整体协同性状态的形变异常指标，从而寻求适合本区预报地震的途径。     

水平荷载情况下深海吸力锚稳定性研究

深海吸力锚基础的极限承载能力是海洋工程结构设计中的一个关键问题,达到极限平衡状态时,吸力锚的极限承载能力与其失稳模式密切相关。本文基于Cou lom b摩擦对原理,给出了一种精确模拟吸力锚承载能力的有限元模型。在该数值模型基础上,利用通用有限元分析软件ABAQU S,研究了系泊点位置、长径比对吸力锚承载能力的影响,并给出了深海吸力锚失稳模式。研究结果表明,系泊点位置极大地影响着吸力锚的极限承载力和稳定性,系泊点位置的变化,会导致吸力锚出现前倾转动、平移滑动和后仰转动等的失稳模式。     

渤海地震是一次地幔地震

１９６９年７月的１８日１３时２４分５０秒，渤海发生了一次７．４级地震，震源深度３５ｋｍ，从地壳厚度分布来看，震中区地壳厚度为２９ｋｍ。所以渤海地震很可能是一次地幔地震。本文从地震活动性、宏观特征分析开始，结合地壳厚度、居里等温面起伏以及地壳上地幔的电性结构和岩石实验材料，对这次地震的深部构造做了论述。     

唐山地震区的地壳结构及大陆地震成因的探讨

本文对唐山地震区的深地震测深资料进行了分析解释。结果表明,该区地壳结构的几个特点,与大陆地区的强震发生有明显关系。 在该震区,相应于不同记录剖面的不同地段,其地壳结构可分为三种类型:类型Ⅰ,类型Ⅱ和类型Ⅲ。每种类型都有它们独特的速度——深度分布。在两种类型的地壳结构边界,或者边界附近,通常存在莫霍界面深度的急剧变化,这种变化可视为某种类型的断层。上述三种类型地壳结构的差异,主要是由上地幔物质向地壳的垂直迁移程度不同造成的。震中区的地壳结构,显示出极高度的物质迁移。然而,应该指出,现今地壳结构是在漫长的地质年代里,经过多次变形积累而形成的。因此,这种结构特征,是表示过去的构造活动痕迹,还是表示现今的构造活动,尚难定论。如我们将不同类型地壳结构的边界,与剥去第三纪的地质构造图进行对比,可以发现,不同地壳结构类型的分界,恰与一定的古老地质构造边界相符合。象北京——天津——唐山这样距离海洋不远的地区,古老的不同种类的小地质体,可能在古代的全球构造运动中聚集起来,并粘结在一起。 1976年唐山主震发生在地壳结构类型Ⅱ内、靠近类型Ⅰ的边界处。因此,唐山地震可能是由于上地幔物质向上迁移的结果,也可能是由于板块的横向挤压、引起脆弱边界重新活化的结果。另有一点是,      

长周期地震观测中的环境干扰

文章根据长周期地震观测的理论和我国763长周期台网的观测实践,分析了长周期观测中的环境干扰。认为长周期观测对环境干扰(主要指温、湿度)敏感,但影响长周期仪器正常工作的关键是潮湿,尤其是记录室的潮湿。文章以763台网优秀台站的观测实践为例,认为只要掌握了环境干扰的规律,在防潮问题上采取有效措施,长周期地震观测是完全可以做到仪器工作稳定,管理得心应手。     

库木库里背斜活动性初步研究

库木库里盆地位于青藏高原北缘,与柴达木盆地一山之隔,是二者的过渡地带,也是高原主体部分向NE扩展的前缘地区;现今构造表现为被3条大型活动构造带(走滑的阿尔金断裂带、东昆仑断裂带和逆冲的祁漫塔格褶皱逆冲系)所夹持。因此,该盆地对于研究青藏高原北缘的构造活动性、活动历史,探讨高原的扩展模式具有十分重要的意义。虽然库木库里盆地南、北两侧均发育活动性很强的大型走滑断裂,但是在盆地中央发育1条大型背斜,走向NWW-SEE,与祁漫塔格褶皱逆冲系和柴达木盆地内的褶皱构造走向一致,说明盆地目前遭受NNE向的挤压。通过对盆地地形横、纵剖面和阶地展布形态的分析,得出背斜有自西向东扩展变形的特征;野外调查和测年结果显示,背斜东段冰川融水形成了大型冰水扇,形成年龄为(87.09±2.31)~(102.4±3.7)ka,进而获得背斜东段自晚更新世以来平均隆升速率的最大值为(2.78±0.28)~(3.28±0.28)mm/a。库木库里盆地整体的活动性很强,在构造上与其北边的柴达木盆地类似,都受控于阿尔金断裂南侧的NNE向的区域挤压作用。     

CRITICAL CROSS POWER SPECTRAL DENSITY FUNCTIONS AND THE HIGHEST RESPONSE OF MULTI-SUPPORTED STRUCTURES SUBJECTED TO MULTI-COMPONENT EARTHQUAKE EXCITATIONS

The highest response of multi-supported structures subjected to partially specified multi-component earthquake support motions is considered. The seismic inputs are modelled as incompletely specified vector Gaussian random processes with known autospectral density functions but unknown cross spectral densities and these unknown functions are determined such that the steady state response variance of a given linear system is maximized. The resulting cross power spectral density functions are shown to be dependent on the system properties, autospectra of excitation and the response variable chosen for maximization. It emerges that the highest system response is associated neither with fully correlated support motions, nor with independent motions, but, instead, specific forms of cross power spectral density functions are shown to exist which produce bounds on the response of a given structure. Application of the proposed results is demonstrated by examples on a ground based extended structure, namely, a 1578 m long, three span, suspension cable bridge and a secondary system, namely, an idealized piping structure of a nuclear power plant.     

用多层变密度直接反演法研究下扬子地区的岩石圈构造

本文用多层变密度直接反演方法研究下扬子地区的岩石圈结构,并讨论其大地构造涵义及演化问题。文中给出了三次样条密度-深度函数的重力正演公式,采用非线性最小二乘优化算法,得到岩石圈内G_5、G_6、G_7三个主要密度界面,它们分别相应于中地壳低密度层底面、Moho面和上地幔顶部低密度层上界面。通过岩石圈结构分析,得到以下结果:杭州湾-海盐-嘉兴-长兴-高淳-芜湖隐伏深断裂是一条北西向的一级构造边界,它对下扬子区的大地构造发展和地壳演化有重要影响。     

用多层变密度直接反演法研究下扬子地区的岩石圈构造

本文用多层变密度直接反演方法研究下扬子地区的岩石圈结构,并讨论其大地构造涵义及演化问题。文中给出了三次样条密度-深度函数的重力正演公式,采用非线性最小二乘优化算法,得到岩石圈内G_5、G_6、G_7三个主要密度界面,它们分别相应于中地壳低密度层底面、Moho面和上地幔顶部低密度层上界面。通过岩石圈结构分析,得到以下结果:杭州湾—海盐—嘉兴—长兴—高淳—芜湖隐伏深断裂是一条北西向的一级构造边界,它对下扬子区的大地构造发展和地壳演化有重要影响。     

一种反演不均匀介质大地电磁测深资料的方法

通过反演受近地表不均匀介质影响的大地电磁测深资料的新方法,可获得较准确的一维地电断面参数、不均匀体埋藏深度及其电性特征。本方法计算速度快,单测ρxy,ρyx可独立反演并相互验证结果的准确性     

普通民众认知与响应汶川M_S8.0地震特征分析——以陕西勉县为例

为进一步了解普通民众对汶川MS8.0地震灾害的认知与响应程度,并为政府防震减灾战略的制定和实施提供科学依据,根据地震灾害认知与响应结构体系设计了普通民众地震灾害认知与响应调查问卷,主要指标体系包括地震灾害知识、防震减灾技能、自救互救和震后信息传播等4个方面。基于权重和赋分,构建了指标体系各部分评价模型和地震灾害认知与响应总体评价模型,以此来计算普通民众地震灾害认知与响应的综合得分。测评结果表明:1)勉县灾区民众认知与响应地震灾害的综合水平与能力均处于不及格水平;2)在认知与响应地震灾害方面,民众的防震减灾技能最差,对地震灾害知识的了解水平次之,而震时及震后自救互救与震后信息传播2个层面稍好一些;3)普通民众对地震知识和观点的理解还停留在表面现象上,对复杂、抽象的地震知识的理解和掌握,特别是对地震知识掌握的准确性和深层次理解还有待于进一步提高;4)从震后民众获取地震灾害信息的主渠道来看,加大政府关于地震灾害知识和防灾技能的宣传力度以及使地震信息公开透明化是降低地震灾害风险的有效途径之一     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

某些电化学分析方法在地震预报中的应用

本文简要介绍了离子选择电极、氧化还原电位(Eh)等电化学分析方法的原理、测试技术。並以1980年甘肃闾井3.4级地震取得的现场观测资料说明,氟离子、硝酸根、Eh等项目在地震前有较明显的前兆反应。以离子选择电极为主的电化学分析方法具有所需仪器设备较少、测定手段简便、便于现场实地观测、易于实现连续自动化等优点。是水化多项观测的切实可行的方法,必将在地震预报的研究中发挥一定的作用。     

初论拱坳构造及其与地震的关系——以鄂西江汉地区为例

本文提出了一个大陆新构造基本类型—拱坳构造。对它的构造特征、性质、应力作用状况以及与地震活动的关系进行了初步论述。以鄂西江汉地区发育的拱坳构造为典型作了较详细的剖析。