Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden's update and 1D derivatives     

Anders V. Christiansen  Esben Auken   《Journal of Applied Geophysics》2004,56(4):247

Dense profile-oriented resistivity data allows for 2D and 3D inversions. However, huge amounts of data make it practically impossible to do full 2D or 3D inversions on a routine basis. Therefore, a number of approximations have been suggested over the years to speed up computations. We suggest using a combination of Broyden's update on the Jacobian matrix with derivatives calculated using a 1D formulation on a parameterized 2D model of locally 1D layered models. The approximations bring down the effective number of 2D forward responses to a minimum, which again gives us the ability to invert very large sections. Broyden's update is not as useful with a parameterized problem as is the case with a smooth minimum structure problem that has been the usual application. 1D derivatives, however, seem to be very effective when initiating a full 2D solution with Broyden's update. We compare the different methods using two different kinds of data on two synthetic models and on two field examples. The most effective and reliable optimization combines 1D derivatives with a full 2D solution and Broyden's update. When using Broyden's update the Jacobian matrix needs to be reset every once in a while. We do this whenever the difference in data residual from the previous iteration is less than 5%. This combined inversion method reduces the computation time approximately a factor of 3 without losing model resolution.  相似文献   

CAPPS1和CAPPS2数值预报模式对比分析   总被引：6，自引：0，他引：6  

李雄 《气象科技》2004,32(6):414-416424

对CAPPS1模式和CAPPS2模式预报广西污染物进行了对比分析。主要对比了其污染扩散原理、气象预报模式、操作系统、运行时间、预报时效及预报效果。分析表明：CAPPS1和CAPPS2模式对各种污染物预报的相关系数均在0．44以上，两种模式预报广西空气污染都是可行的；两种模式均对SO2的预报效果最好，NO2次之，PM10相对最差；CAPPS2模式对各种污染物的预报值与实测值相关系数均比CAPPS1高，其中NO2的相关系数提高幅度最大；CAPPS2模式对NO2和PM10的预报准确率比CAPPS1偏高，而对SO2的预报准确率比CAPPS1偏低；CAPPS2预报首要污染物正确率比CAPPS1高出11个百分点，而预报首要污染物等级正确率则高出9个百分点。  相似文献   

赣南中侏罗世流纹岩地球化学及成因研究：上地壳成因的的微量元素和Pb-Nd-Sr同位素地球化学制约   总被引：4，自引：1，他引：4  

章邦桐凌洪飞  孔兴功 《岩石学报》2004,20(3):511-520

江西南部龙南-寻邬地区余田群菖蒲组的流纹岩形成于中侏罗世(164．8．1Ma～174．9Ma)，为双峰式火山岩组合的酸性端元，属弱过铝质(ACNK平均为1．04)钙碱性系列(σ平均为1．90)。稀土元素含量高(平均为278．9μg/g)，轻稀土富集(IMH平均为7．94)，Eu亏损中等(δEu平均为0．41)。微量元素以富集Rb，Ba，Th，U，Zr，Y，Pb，而亏损Sr，Nb为特征。Pb，Nd，Sr同位素组成为：(^206Pb/^204Pb)t=17．89～18．58，(^207Pb/^204Pb)t=15．58～15．70，(^208Pb/^204Pb)t=37．94～38．82；εNd(t)=-7．44～-11．9；ISr=0．71126～0．71228。对赣南流纹岩微量元素采用蛛网图和模糊聚类分析方法进行对比研究的结果及Pb-Nd-Sr同位素相关特征为其上地壳成因提供了地球化学佐证。赣南流纹岩是由进入上地壳内高位岩浆房的石英拉斑玄武岩浆的底侵作用直接使上地壳部分熔融形成的产物。  相似文献   

用于新生代定年的Ar-Ar法标准样品候选样品初测结果   总被引：5，自引：1，他引：5  

桑海清  丁林  戴檀谟 《矿物岩石地球化学通报》2004,23(1):24-27

Ar Ar法定年的特点是必须要有一套年龄从小到大的标准样。迄今用于新生代定年的国内标准样极少。为满足新生代矿物Ar Ar法定年的需要,初选了一个标准样候选者BT 1透长石。样品总重366g,粒级6080目,纯度100%,缩分为100瓶,每瓶3 66g。Ar Ar阶段加热法初测结果为:全部12个阶段给出的总气体年龄为30 8±0 9Ma,412阶段视年龄十分接近,年龄谱平坦,对应的39Ar析出量达96%,坪年龄为29 6±0 4Ma,等时年龄为29 6±0 6Ma,MSWD=1 01。40Ar/36Ar初始值为293 6±3 9,与尼尔值295 5相当。重复测定结果为:全部气体年龄为31 0±0 9Ma(全部9个阶段),坪年龄为29 5±0 4Ma,等时年龄为29 4±0 6Ma,MSWD=1 94,40Ar/36Ar初始值为282 2±6 3。这些结果表明,BT 1透长石不含过剩氩,作为新生代定年的Ar Ar法标准样品候选者是较为理想的。  相似文献   

Paleozoic serpentinite-enclosed chromitites from Tehuitzingo (Acatlán Complex, southern Mexico): a petrological and mineralogical study     

J. A. Proenza  F. Ortega-Gutirrez  A. Camprubí  J. Tritlla  M. Elías-Herrera  M. Reyes-Salas 《Journal of South American Earth Sciences》2004,16(8):129

The serpentinites and associated chromitite bodies in Tehuitzingo (Acatlán Complex, southern Mexico) are in close relationship with eclogitic rocks enclosed within a metasedimentary sequence, suggesting that the serpentinites, chromitites and eclogitic rocks underwent a common metamorphic history.Primary chromites from the chromitite bodies at Tehuitzingo are of refractory-grade (Al-rich) and have a chemical composition similar to that expected to be found in an ophiolitic environment. The chromite grains in chromitites and serpentinites are systematically altered to ‘ferritchromite’. The alteration trend is usually characterized by a decrease in the Al, Mg and Cr contents coupled by an increase in Fe³⁺ and Fe²⁺.The Tehutizingo chromitites have low Platinum Group Elements (PGE) contents, ranging from 102 to 303 ppb. The chondrite-normalized PGE patterns are characterized by an enrichment in the Ir-subgroup elements (IPGE=Os, Ir, Ru) relative to the Pd-subgroup elements (PPGE=Rh, Pt, Pd). In addition, all chromitite samples display a negative slope from Ru to Pd [(Os+Ir+Ru)/(Pt+Pd)=4.78−14.13]. These patterns, coupled with absolute PGE abundances, are typical of ophiolitic chromitites elsewhere. Moreover, all the analyzed samples exhibit chondrite-normalized PGE patterns similar to those found for non-metamorphosed ophiolitic chromitites. Thus, the PGE distribution patterns found in the Tehuitzingo chromitites have not been significantly affected by any subsequent Paleozoic high-pressure (eclogite facies) metamorphic event.The chondrite-normalized PGE patterns of the enclosing serpentinites also indicate that the PGE distribution in the residual mantle peridotites exposed in Tehuitzingo was unaffected by high-pressure metamorphism, or subsequent hydrothermal alteration since the serpentinites show a similar pattern to that of partially serpentinized peridotites present in mantle sequences of non-metamorphosed ophiolites.Our main conclusion is that the chromitites and serpentinites from Tehuizingo experienced no significant redistribution (or concentration) of PGE during the serpentinization process or the high-pressure metamorphic path, or during subsequent alteration processes. If any PGE mobilization occurred, it was restricted to individual chromitite bodies without changing the bulk-rock PGE composition.Our data suggest that the Tehuitzingo serpentinites and associated chromitites are a fragment of oceanic lithosphere formed in an arc/back-arc environment, and represent an ophiolitic mantle sequence from a supra-subduction zone, the chemical composition of which remained essentially unchanged during the alteration and metamorphic events that affected the Acatlán Complex.  相似文献   

Petrology and Geochemistry of Early Cretaceous Bimodal Continental Flood Volcanism of the NW Etendeka, Namibia. Part 1: Introduction, Mafic Lavas and Re-evaluation of Mantle Source Components   总被引：9，自引：4，他引：9  

EWART  A.; MARSH  J. S.; MILNER  S. C.; DUNCAN  A. R.; KAMBER  B. S.; ARMSTRONG  R. A. 《Journal of Petrology》2004,45(1):59-105

The bimodal NW Etendeka province is located at the continentalend of the Tristan plume trace in coastal Namibia. It comprisesa high-Ti (Khumib type) and three low-Ti basalt (Tafelberg,Kuidas and Esmeralda types) suites, with, at stratigraphicallyhigher level, interstratified high-Ti latites (three units)and quartz latites (five units), and one low-Ti quartz latite.Khumib basalts are enriched in high field strength elementsand light rare earth elements relative to low-Ti types and exhibittrace element affinities with Tristan da Cunha lavas. The unradiogenic²⁰⁶Pb/²⁰⁴Pb ratios of Khumib basalts are distinctive, most plottingto the left of the 132 Ma Geochron, together with elevated ²⁰⁷Pb/²⁰⁴Pbratios, and Sr–Nd isotopic compositions plotting in thelower ¹⁴³Nd/¹⁴⁴Nd part of mantle array (EM1-like). The low-Tibasalts have less coherent trace element patterns and variable,radiogenic initial Sr (  相似文献   

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media     

E. Chénier  R. Eymard  X. Nicolas 《Computational Geosciences》2004,8(2):163-172

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.  相似文献   

Regional moment tensor uncertainty due to mismodeling of the crust     

Jan ílený 《Tectonophysics》2004,383(3-4):133-147

The retrieval of earthquake moment tensor (MT) requires the response of the medium, in which seismic waves travel from the hypocenter to the stations, to be known. In inverting long-period (LP) seismic data (teleseismic and LP regional records), a gross earth model is sufficient; with decreasing periods, a more detailed model is needed. This is the case when waveforms of weak earthquakes at regional distances are to be inverted. Regional moment tensors (RMTs) of mostly Mediterranean earthquakes are determined on a routine basis by the Swiss Seismological Survey (SED) by using averaged models of the earth's crust. By inverting broad-band records of the M_w=4.8 earthquake near Udine, N Italy, on Feb. 14, 2002, we tested the sensitivity of the MT solution with respect to possible errors in the earth model used and in the location of the hypocenter depth. We perturbed the P and S velocities and the thickness in the 1-D earth model in the range from 3% to 30% of the parameter values and constructed estimates of confidence regions of the MT and error bars of the source time function (STF) and scalar moment in three frequency bands. Similarly, these error characteristics were determined assuming a mislocation in the hypocenter depth. We found that, in the band of periods from 25 to 50 s, the mechanism is resolved well (at the confidence level 95% at least) up to an earth model uncertainty of 30%, in the passband 10–25 s up to about 10%, but it is undetermined completely at periods of 5–10 s. An error in hypocenter depth of as much as double the value reported by the location procedure does not destroy the resolution of the mechanism at periods above 10 s. In the RMT catalog of the SED, earthquakes of M_w greater than about 3.5 are processed at periods above 30 s; thus, the solutions for these events are robust with respect to a possible uncertainty in the earth model used. Mechanisms of weaker earthquakes, retrieved from short periods, should be interpreted with caution.  相似文献   

HJ-1复合型钻井液的研究及应用          下载免费PDF全文

王政敏  牛振中  颜志华 《探矿工程》2004,31(5):48-50

分析和探讨了HJ-1复合型钻井液对复杂地层护壁的作用机理及护壁性能试验情况，介绍了其在寨上矿区复杂地层钻探中的成功应用。  相似文献   

MODIS 1B影像几何纠正方法研究及软件实现   总被引：24，自引：0，他引：24  

蒋耿明  刘容高  牛铮  庄大方 《遥感学报》2004,8(2):158-164

MODIS影像是一种新型和重要的数据。对MODIS 1B影像几何畸变原因进行了深入分析，选择了一种合适的纠正方法。对于1km分辨率MODIS 1B影像，直接采用1km分辨率的空间坐标进行几何纠正；对于250m和500m分辨率的MODIS 1B影像，先用三次样条曲线对坐标进行插值生成同分辨率的坐标，然后利用坐标插值结果对其进行几何纠正。由于MODIS影像在空问分布上的特殊性，采用前向和后向映射相结合的方式确定纠正后某一像素点在原始影像中的位置。根据该位置的条带重叠度，可以确定参与计算像素个数、搜索窗口的大小以及窗口的精确位置；采用归一化反距离加权插值法计算纠正后像素点的属性值。上述优化算法不仅保证了纠正后影像的质量，而且提高了数据处理速度。作者在Visual C 6．0环境下开发实现了上述算法。从坐标插值和几何纠正结果分析，无论是数据处理速度还是纠正后的影像质量和精度均达到要求。  相似文献