起伏地形上规则二度体复重磁场正演和直接反演 连载3

上接 119页9　反演方法推导中的有关公式为了对复场展开式作简化表示 ,将会多次用到2个递推公式 ,这里预先导出。设 :x1=∑ni=1σi,　x2 =∑ni,j=1σiσj,(i 相似文献   

地球化学样品分析质量的计算机监控程序

本文提供的计算机监控程序可以用于 1 :5万和 1 :2 0万化探样品分析质量监控和覆盖区多目标地球化学样品分析质量监控。1　程序的编制和设计　　主体程序如下INPUT″请输入起始样号″,BR=B- 1 (注 :实际工作中要求以 50件样品为一批 )DIM T( R+ 55) ,S( R+ 50 ) ,y ( 1 0 0 0 )FOR I=R+ 1 TO R+ 50 :READ T ( I) :NEXT IFOR n=1 TO 5PRINT″请输入第″;n;″组数据 ,每输入一个数据后按回车″2 0 FOR I=R+ ( n- 1 ) * 1 0 + 1 TO R+ ( n- 1 ) * 1 0 + 1 0 :INPUT S( I)NEXT IPRINT″请检查数据是否正确 ,正确输…     

起伏地形上规则二度体复重磁场正演和直接反演

上接 38页5　复场直接积分法和板状体复磁场公式5 .1　薄板状体复磁场公式5 .1.1　水平薄板设 :水平薄板厚为Ｄ、宽为 2ｂ ,左端点为σ1,右端点为σ2 =σ1+2ｂ ,板内变点σ =σ1+ｌ,0≤ｌ≤2ｂ。该板复场Ｔ(Ｓ)相当小水平圆柱体的复场沿水图 6　水平薄板断面示意平方向积分 (图 6)。在 ( 1 4 7)式中 ,取Ｍｍ =ＤｄｌＭｓ,　σ0 =σ1+ｌ ,( 5 1)则水平薄板状体复磁场为Ｔ(Ｓ) =∫2ｂ0μ04π2ＤＭｓｄｌ{Ｓ- (σ1+ｌ) } 2 =μ04π2Ｍｍ1σ1-Ｓ- 1σ2 -Ｓ ,( 5 2 )Ｍｍ=ＤＭｓ。 ( 5 3)5 .1.2　有限延深倾斜薄板设板宽为 2ｂ( Ｉｍσ1)、…     

重力勘探任意点解释法初探

在矿体的xoh剖面图中,我们将x分解成:x_2=x_1+△x;x_3=x_2+△x=x_1+2△x;x_n=x_(n-)。~1+△x=x_1+(n-1)△x。式中,△x—两相邻点之点距。然后利用g的导数立微分方程:g~'_1=(g_1-g_2)/△x和联立微分方程。再化为代数方程组。解联立方程组,便可求出矿体的埋深h、坐标x_1及其他物理量。     

厌氧环境中硫代As(V)在水-矿物界面的吸附研究

<正>水体中砷的迁移转化行为受其赋存形态影响。一般有氧水体和缺氧水体中,常见的无机砷形态为砷酸根(HAs O42-,H2As O4-和As O43-)或亚砷酸(H3As O3);在一些含硫的缺氧水体中,元素硫参与砷的地球化学循环,与之形成含砷硫化物沉淀或者可溶的硫代砷(thioarsenic species)。硫代砷形态十分复杂,按照砷的价态可分为硫代As(V)(thioarsenate,H3AsVOx S4-x,x=0~3)和硫代As(III)(thioarsenite,H3AsIIIOx S3-x,x     

Study on the Crystal of Daomanite

It was not possible to carry out a complete analyses of crystal, as the experiment by Ding and Shi et al.. It's analysis precision R=0.25 or more big than this, which value are not satisfied for single crystal study, but we through many test and found the best: [R(int)=14.5%]. The final fullmatix least-squares refinement on F2 converged to R1=0.0791 and wR2=0.1864 for 704 observed reflections [I 3 2s(I)]. Daomanite is orthorhombic system, space group Cmc21, a=3.7520(8))?, b=15.844(4) ?, c=5.8516(12) ?, α=β=γ=90°. V=347.86(14)?3, Z=4. Daomanite chemical formula is Cu Pt AsS 2. Idealized composition Me+M2+M2+S2=CuS ·PtA s S. There is no other similar mineral in the world.     

阜平群下部太古代变质岩石的REE、Rb-Sr和Sm-Nd年龄及其意义

对采自阜平群下部16个变质岩石样品进行了Sm-Nd年龄测定,部分样品还进行了Rb-Sr年龄测定和REE分析。Sm-Nd同位素结果形成两条等时线。一条等时线具有年龄为:2790±171(2σ)Ma,I_(Nd)=0.50917±16(2σ),ε_(Nd)(t)=±1.6;另一条等时线,t=1949±84(2σ)Ma,I_(Nd)=0.50986±5(2σ),ε_(Nd)(t)=-5.0±1.0。Rb-Sr同位素数据分散,参考等时线年龄为2070Ma。     

TI-59计算器磁法程序介绍

、无限延深倾斜厚板磁法正演1.计算公式Z一ZJ sin。〔音5‘。(。一‘),n(x一b)2 hZ行千句面干玉豆十eos(日一i)tan一’ ZbhxZ一bZ hZ上式可改写为: ..几 一 g 儿L 、、刀/za=命,S‘n日(Za以Y为单位)。05甲一i心tg(。一i)Ln(卜 4bx(x b)“ li“ ZbhxZ一bZ hZ〕 其中:J:视倾角(度),垂直厚板走向磁化强度分量(10一6CGSM);i:视磁倾角(度);日:厚板h:厚板顶面深度(米);b:厚板半宽度(米)2.操作说明步骤}1{2… 一3:…6{2几dLRNCP000 00LRN AR/S 一三R/SR/S结果值清除程序存贮器进入学习态输入计算Za程序骨退出学习态开始输入数据输…     

红河石Ca18(□,Ca)2Fe2+Al4(Fe3+,Mg,Al)8(□,B)4BSi18O69(O,OH)9——符山石族新矿物

作为符山石族矿物的新成员,红河石(Hongheite,IMA 2017-027新矿物),Ca18(,Ca)2Fe2+Al4(Fe3+,Mg,Al)8(,B)4BSi18O69(O,OH)9发现于个旧世界级Sn-多金属矿田东北缘、与马拉格Sn矿床毗邻的北沙冲花岗岩(77.43Ma)内矽卡岩中。红河石常呈横径达4~25mm的放射状针-柱状集合体产出。当位于晶洞中时,红河石则呈发育良好的自形柱状晶体(0.5~4.0mm长,0.3~1.0mm宽)产出。与红河石共生的矿物见有赛黄晶、萤石、斧石-(Fe)、硅硼钙石、枪晶石、硼锡钙石、石英和羟鱼眼石-(K)等。红河石为墨绿色,条痕浅灰绿色,玻璃光泽,性脆,断口不规则。主要的晶面是:{100}、{110}、{101}和{001}。红河石的显微硬度:988.3N/mm2,相当于摩氏硬度6~7。其实测密度与计算密度分别是3.446g/cm3和3.423g/cm3。红河石一轴正晶,No=1.720(2),Ne=1.725(2);多色性弱。红河石的化学成分:SiO235.85%;TiO20.01%;Al2O311.00%;Fe2O37.92%;FeO2.14%;CaO 33.57%;MnO 0.42%;MgO 3.48%;B2O32.82%;Cr2O30.01%;Na2O 0.01%;F 0.40%(F≡O-0.17);Cl 0.14%(Cl≡O-0.03);H2O 0.75%,总量98.32%。依据晶体结构精测和Si在单位分子式中的原子数(即Si=18 apfu),计算和书写的红河石简化晶体化学式:Ca18(,Ca)2Fe2+Al4(Fe3+,Mg,Al)8(,B)4BSi18O69(O,OH)9。其三条最强粉晶线[d(?)(I/I0)(hkl)]为:2.9289(47)(004),2.7661(100)(342)和2.6079(68)(243)。红河石属四方晶系,空间群为P4/nnc,晶胞参数:a=15.667(3)?,c=11.725(1)?,V=2878(1)?3,Z=2。红河石晶体结构精测的R因子为0.063。红河石殊异于为已知的符山石族矿物种,在于其X(4)位以空位()为主、Y(3)位以Fe3+居优和T(2)位被B所占。顺便对符山石族矿物晶体-化学式的计算与书写予以讨论并提出建议。     

实验变差函数估计中的某些问题综述

变差函数(Variogram)又称半变差函数(Semivariogram),是研究区域化变量的基本工具,记Z(x)为所研究的区域化变量(如矿石品位),γ_z为其变差函数,则变差函数的数学定义为γ_z(x,h)=1/2E{[Z(x+h)-Z(x)]~2}(1)(1)式又称为理论变差函数.按一般的统计学观点,应由下列统计基来估计γ_z,即γ_2(x,h)=1/2N sum from i=1 to N′(i/1)[Z_i(x+h)-Z_i(x)]~2(2)但是在地质问题的应用中,通常是不可能由(2)式来估计γ_z的.因为,在确定的空间位置点X_0上只能有Z(X)的一次观测值,即Z(X_0).因而,只得利用下式来估计γ_z(但这样做仍是有一定条件的).     

局部重磁场源全方位成像理论概要(续)

作者近年来首创了"复杂条件下局部重磁场源全方位成像"理论体系。在本文中,概要阐述该理论体系中下延有限二度体复重磁场级数正演通式、复场模值全方位延拓和全方位反演等。     

重磁数据的线性滤波与反演失真问题

文中对国内研制的4种主要线性滤波方法进行了对比分析。选出一种具代表性的滤波方法,用计算机模拟观测场,探讨了在复杂干扰下重磁异常因滤波引起的失真,并就重磁异常在因滤波而引起失真的情况下如何反演解释,提出了可行的办法和建议。     

局部重磁场源全方位成像理论概要

本文概要阐述了复杂条件下局部重磁场源全方位成像理论体系的实用价值与学术意义、球坐标系内三度体重磁场球谐级数正演通式、三度体重磁场全方位延拓和全方位反演等。     

复小波变换反演重力异常

对国外最新研究位场连续复小波变换理论的一系列文章进行了系统归纳、推导与总结,推出了板状体重力场的小波变换系数解析式.对齐次场源形成的重力场进行柯西小波变换,等效于重力场向上延拓和求导两种运算相结合,并服从小波变换双尺度定律.利用小波变换系数模沿极值轴线的变化特征同时反演重力场源的几何形状、埋深、倾角.独立编写了所有程序,并对不同理论模型进行了大量计算,对所有结果进行了分析与讨论.此法为重力场反演提供了一条新的途径,实用于重磁资料自动化处理、反演与解释.     

勘查油气田化探新方法——土壤相态烃测量法

烃气测量是油气化探最主要的方法。目前应用的烃气测量方法有两类:一类是测量土壤中游离烃气如K—V指纹法、吸附丝法、物上气法及直接取土壤气测量法等。此类方法由于需特殊设备及壤气测量方法本身固有的缺陷,目前还未得到普遍应用,另一类为土壤吸附烃气法即国内外目前普遍使用的酸解烃气法、此法在多数油气田上取得好效果,但也有相当多的地区如黄土、荒漠及水网发育等特殊景观区效果不好或不显著。本文介绍一种新的烃气测量法──土壤相态烃测量法。这是用一种特殊的技术,只选择性提取由油气藏渗漏过程中保存在运移通道围岩及近地表土壤各种次生矿物结构中以"穴居"态存在的轻烃组分。"穴居"态烃为累积性的指标,不受外界温度、修水等的影响,保存较牢固,测定结果稳定。相态烃法避免了酸解烃法存在的生物成因烃、矿物同生烃及碳酸盐的干扰影响,提高了测量指标的信噪比,突出了应有的异常。本方法已在黄土厚覆盖区两个已知油气田和一个风成沙覆盖的断裂发育油田上试验获得成功,结果表明相态烃法优于酸解烃法。     

重力异常三维反演——视密度成像方法技术的应用

将起伏地形上测量的重力异常,延拓至平均水平面;对平面上的重力场进行不同深度层源的切割分离;再将各层的重力场向下延拓至相应的深度,并反演出该深度层的密度结构.以此得到的密度,反映该深度层密度的近似分布,故称为视密度成像.在重力场的曲面向下延拓过程中,采用了新的延拓方法——迭代法,延拓深度大,延拓过程稳定;且不需解线性代数方程组,避开了制约三维反演实用化的计算时间的瓶颈问题.在主频1.73 GHz的微机上,完成71×81数据点反演的计算时间为10 s.以新疆一个地区的实例说明这种方法的反演效果.     

The Rapid Inversion of 3-D Potential Field and Program Design

The application of three-dimensional inversion of gravity and magnetic fields is very important not only in geophysical researches, but also in the study of geological structures. A formula of potential field in frequency-domain, developed by Parker in 1973, can be used as a rapid and effective algorithm in gravity and magnetic inversion. The technique has been improved then by Oldenburg, Sprenke, Feng and others.In addition to a brief introduction of Parker's algorithm and its applications, this paper includes the following five parts: basic computational techniques, inversion of single layer, convergence and constraints, simultaneous inversion for density and topography as well as inversion of multilayers. The authors present relevant practical iterative formulas and its varieties when density distribution varies with depth in linear or exponential relation. In order to maintain computation stability and speed up iteration convergence, some approaches are taken in the program design, for instance shifting lower interface of the studied layer, inverting corrections of topography, reducing grid boundary effects and utilizing low-pass filter. With the consideration of the nonuniqueness of the inversion, a method of using seismic data to constrain the range ofpossible models is discussed. It is pointed out that the density variation generates less effects than those of topography on the spectrum of gravity anomaly in second order. Therefore density contrast and topography can be inverted simultaneously by an alternative weighting iteration. By analogy, the inversion of multilayer model can be done in the above procedure. An approach of model decomposition is useful in the computation of multilayer model. The techniques discussed in the present paper for gravitational field are also valid for magnetic field.     

重磁异常的人机联作校正－迭代反演

本文提出了三度体重磁异常的人机联作校正－迭代反演方法。该方法用二度半组合多边形棱柱体来逼近三度体，从而把三度体重磁异常反演问题转化为二度半体的反演问题；为了消除组合体迭加场的影响，该方法采用了一种校正－迭代技术。理论模型反演计算表明，该方法实际可行。     

重磁异常的人机联作校正－迭代反演

本文提出了三度体重磁异常的人机联作校正－迭代反演方法。该方法用二度半组合多边形棱柱体来逼近三度体，从而把三度体重磁异常反演问题转化为二度半体的反演问题；为了消除组合体迭加场的影响，该方法采用了一种校正－迭代技术。理论模型反演计算表明，该方法实际可行。