共查询到20条相似文献，搜索用时 390 毫秒
1.
In this paper, we study the problem of constructing a smooth approximant of a surface defined by the equation z = f(x 1, x 2), the data being a finite set of patches on this surface. This problem occurs, for example, after geophysical processing such as migration of time-maps or depth-maps. The usual algorithms to solve this problem are picking points on the patches to get Lagrange's data or trying to get local junctions on patches. But the first method does not use the continuous aspect of the data and the second one does not perform well to get a global regular approximant (C 1 or more). As an approximant of f, a discrete smoothing spline belonging to a suitable piecewise polynomial space is proposed. The originality of the method consists in the fidelity criterion used to fit the data, which takes into account their particular aspect (surface's patches): the idea is to define a function that minimizes the volume located between the data patches and the function, and which is globally C k. We first demonstrate the new method on a theoretical aspect and numerical results on real data are given. 相似文献
2.
We introduce a surface approximation technique to address the problem of fitting a surface to a given set of curves. The originality of the method lies in its ability to take into account the continuous aspect of the data, and also in the possibility to arbitrarily select the regularity (C0, C1, or higher) of the approximant obtained. We demonstrate the efficiency of the approach by constructing a bathymetry map of the Marianas trench based upon a set of SONAR (SOnic Navigation And Ranging) bathymetry ship track data. 相似文献
3.
We introduce a surface approximation technique to address the problem of fitting a surface to a given set of curves. The originality of the method lies in its ability to take into account the continuous aspect of the data, and also in the possibility to arbitrarily select the regularity (C0, C1, or higher) of the approximant obtained. We demonstrate the efficiency of the approach by constructing a bathymetry map of the Marianas trench based upon a set of SONAR (SOnic Navigation And Ranging) bathymetry ship track data. 相似文献
4.
The problem of determining capillary pressure functions from centrifuge data leads to an integral equation of the form
a
x
K(x,t)f(t)dt=g(x),x[a,b],(1)where the kernel K is known exactly and given by the underlying mathematical model. g is only known with a limited degree of accuracy in a finite and discrete set of points x
1,...,x
M
. However, the sought function f(t) is continuous. By the nature of the right-hand side, g(x), equation (1) is a discrete inverse problem which is ill-posed in the sense of Hadamard [9]. By a parameterization of the sought function, equation (1) reduces to a system of linear equations of the form
Ac=b+ ,where b is the observation vector and A arises from discretization of the forward problem. is the error vector associated with b, and c contains the model parameters. The matrix A is usually ill-conditioned. The ill-conditioning is closely connected to the parameterization of the problem [23].In this paper a semi-iterative regularization method for solving the Volterra integral equation in the 2-norm, namely, Brakhage's -method [2], is investigated. The iterative method is tested on synthetically generated, and on experimental data. 相似文献
5.
In this study, we present ac susceptibility measurements for a synthetic and a natural hemo-ilmenite (HI) solid solution (x)FeTiO3–(1 − x)Fe2O3 with compositions x = 0.87(1) and 0.88(8), respectively. The focus of the investigation is the magnetic ordering at the Curie temperature T
C and the spin-glass-like freezing at the freezing temperature T
f. The sharpness of T
C for the synthetic solid solution with well-defined structure indicates the chemical homogeneity of the solution, whereas
the disperse magnetic ordering of the natural solid solution reveals inhomogeneities described as spin glass system ofvariations
in composition x. The frequency dispersion of T
f was determined between 10 Hz and 10 kHz and was found to obey a dynamic scaling power law. The relaxation rates deviate by
five orders of magnitude where the synthetic solid solution exhibits ω0 = 3(1) × 104 Hz and the natural one 5.5 × 109 Hz. The strong deviation is attributed to the difference in the ordered state above T
f. These findings provide an insight into the cooling-rate effects of natural solid solutions and how magnetic thermodynamics
can be used to probe the chemical homogeneity of such systems. 相似文献
6.
The Dharwar craton, Cuddapah basin and the Godavari graben characterise three diverse geological and tectonic settings in
the peninsular shield of India. Owing to their contiguous proximity, they offer a unique opportunity to document differences,
if any, in their seismic wave attenuation characteristics that might have a bearing on the seismogenic nature of the crust
in a craton, basin and a rift-like graben structure. An attempt is made here to bring out these differences using constraints
from coda-Q. We considered local earthquakes with epicentral distances ranging from␣14 to 150 km recorded at the digital broadband stations
at Dharwar (DHD), Cuddapah (CUD) and Kothagudem (KGD) regions to derive the frequency-dependent coda-Q relations. Using the single scattering method, we obtained the frequency-dependent Q
C relationship (Q
C = Q
0
f
n
)for each of the three geological units separately: DHD: Q
C = (730.62 ± 0.09)f
(0.54 ± 0.01); CUD: Q
C = (535.06 ± 0.13)f
(0.59 ± 0.01) and KGD: Q
C = (150.56 ± 0.08)f (0.91 ± 0.01). The Q
C values obtained for all the three sub-regions show moderate to strong frequency dependence and essentially reflect the level
of crustal heterogeneities to varying degrees. 相似文献
7.
Rhythmic zoning is ubiquitous in igneous rocks. Based on polymerization and regular solution models, a non-ideal, disequilibrium
and non-linear interface equation has been established to describe the process of crystal growth:f = 1/[1 + (β/x
s
)(1 -x
s
)exp(-W / RT)(1 - 2f)] wheref and xs represent respectively the mole fraction of a component in the crystal and melt at the interface;W the total exchange energy;R the gas constant; andT the temperature; β=k
B
/k
A
, withk
A
andk
B
representing respectively the rate constants of components A and B. Results of numerical simulation of the equation demonstrate
that a domain of triple valuedness exists ifW / RT < -2. This model, together with the mass balance equation, explains reasonably the rhythmic phenomena in silicate solid solutions,
indicating that self-organization is responsible for this process during mineral growth. 相似文献
8.
Cation tracer diffusion coefficients, DMe *, for Me=Fe, Mn, Co and Ti, were measured using radioactive isotopes in the spinel solid solution (Ti x Fe 1−x )3−δO4 as a function of the oxygen activity. Experiments were performed at different cationic compositions (x=0, 0.1, 0.2 and 0.3) at 1100, 1200, 1300 and 1400 °C. The oxygen activity dependence of all data for DMe * at constant temperature and cationic composition can be described by equations of the type DMe *=D∘ Me[V]. CV·a O2 2/3+DMe[I] ∘·a O2 −2/3·DMe[V] ∘ and DMe[I] ∘ are constants and CV is a factor of the order of unity which decreases with increasing δ. All log DMe * vs. loga O2 curves obtained for different values of x and for different temperatures go through a minimum due to a change in the type of point defects dominating the cation diffusion with oxygen activity. Cation vacancies prevail for the cation diffusion at high oxygen activities while cation interstitials become dominant at low oxygen activities. At constant values of x, DMe[V] ∘ decreases with increasing temperature while DMe[I] ∘ increases. 相似文献
9.
The practice of stochastic simulation for different environmental and earth sciences applications creates new theoretical
problems that motivate the improvement of existing algorithms. In this context, we present the implementation of a new version
of the direct sequential co-simulation (Co-DSS) algorithm. This new approach, titled Co-DSS with joint probability distributions,
intends to solve the problem of mismatch between co-simulation results and experimental data, i.e. when the final biplot of
simulated values does not respect the experimental relation known for the original data values. This situation occurs mostly
in the beginning of the simulation process. To solve this issue, the new co-simulation algorithm, applied to a pair of covariates
Z
1(x) and Z
2(x), proposes to resample Z
2(x) from the joint distribution F(z
1,z
2) or, more precisely, from the conditional distribution of Z
2(x
0), at a location x
0, given the previously simulated value
z1(l)(x0)z_{1}^{(l)}(x_{0})
(F(Z2|Z1=z1(l)(x0)F(Z_{2}|Z_{1}=z_{1}^{(l)}(x_{0})
). The work developed demonstrates that Co-DSS with joint probability distributions reproduces the experimental bivariate
cdf and, consequently, the conditional distributions, even when the correlation coefficient between the covariates is low. 相似文献
10.
Crystals of hydronium jarosite were synthesized by hydrothermal treatment of Fe(III)–SO4 solutions. Single-crystal XRD refinement with R1=0.0232 for the unique observed reflections (|Fo| > 4F) and wR2=0.0451 for all data gave a=7.3559(8) Å, c=17.019(3) Å, Vo=160.11(4) cm3, and fractional positions for all atoms except the H in the H3O groups. The chemical composition of this sample is described by the formula (H3O)0.91Fe2.91(SO4)2[(OH)5.64(H2O)0.18]. The enthalpy of formation (Hof) is –3694.5 ± 4.6 kJ mol–1, calculated from acid (5.0 N HCl) solution calorimetry data for hydronium jarosite, -FeOOH, MgO, H2O, and -MgSO4. The entropy at standard temperature and pressure (So) is 438.9±0.7 J mol–1 K–1, calculated from adiabatic and semi-adiabatic calorimetry data. The heat capacity (Cp) data between 273 and 400 K were fitted to a Maier-Kelley polynomial Cp(T in K)=280.6 + 0.6149T–3199700T–2. The Gibbs free energy of formation is –3162.2 ± 4.6 kJ mol–1. Speciation and activity calculations for Fe(III)–SO4 solutions show that these new thermodynamic data reproduce the results of solubility experiments with hydronium jarosite. A spin-glass freezing transition was manifested as a broad anomaly in the Cp data, and as a broad maximum in the zero-field-cooled magnetic susceptibility data at 16.5 K. Another anomaly in Cp, below 0.7 K, has been tentatively attributed to spin cluster tunneling. A set of thermodynamic values for an ideal composition end member (H3O)Fe3(SO4)2(OH)6 was estimated: Gof= –3226.4 ± 4.6 kJ mol–1, Hof=–3770.2 ± 4.6 kJ mol–1, So=448.2 ± 0.7 J mol–1 K–1, Cp (T in K)=287.2 + 0.6281T–3286000T–2 (between 273 and 400 K). 相似文献
11.
Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations 总被引：3，自引：0，他引：3
Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature‐induced grain coarsening. Despite this large‐scale trend, a considerable grain size scatter occurs on the outcrop‐scale indicating local influence of second‐order effects such as thermal perturbations, fluid flow and second‐phase particles. Second‐phase particles, whose sizes range from nano‐ to the micron‐scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron‐scale second‐phase particles, or to samples containing a large number of nano‐scale particles. The obtained data set suggests that the second phases induce a temperature‐controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D = C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second‐phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian‐type distributions of the pure samples. In contrast, fluid‐enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable δ18O and δ13C isotope ratios in fluid‐affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid‐induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1. To summarize, considerable grain size variations of up to one order of magnitude can locally result from second‐order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies. 相似文献
12.
Qiang He Xi Liu Xiaomin Hu Sicheng Li Hejing Wang 《Physics and Chemistry of Minerals》2011,38(10):741-752
The solid solution between lead fluorapatite and lead fluorvanadate apatite, Pb10[(PO4)6-x
(VO4)
x
]F2 with x equal to 0, 1, 2, 3, 4, 5 and 6, was synthesized by solid-state reaction at 1 atm and 700°C for 72 h and characterized by
scanning electronic microprobe, electronic microprobe analysis, micro-Raman spectroscopy, and powder X-ray diffraction. The
volume-composition relationship at ambient temperature does not show significant deviation from the Vegard’s Law. The Raman
spectrum data suggest that both P and V are identical on a C
s
site and both end-members show no apparent factor-group effect. The Raman frequency shift of the symmetric stretching vibration
is linearly dependent on the composition. High temperature X-ray diffraction data, up to 600°C, suggest that the thermal expansion
coefficients α
a
, α
c
, and α
V
also vary linearly with the compositions of the apatites. 相似文献
13.
Frédéric Hatert André-Mathieu Fransolet Walter V. Maresch 《Contributions to Mineralogy and Petrology》2006,152(4):399-419
In order to assess the geothermometric potential of the Na2(Mn2−2x
Fe1+2x
)(PO4)3 system (x = 0–1), which represents the compositions of natural weakly oxidized alluaudites, we performed hydrothermal experiments between 400 and 800°C, at 1 kbar, under an oxygen fugacity (f(O2)) controlled by the Ni–NiO (NNO), Fe2O3–Fe3O4 (HM), Cu2O–CuO (CT), and Fe–Fe3O4 (MI) buffers. When f(O2) is controlled by NNO, single-phase alluaudites crystallize at 400 and 500°C, whereas the association alluaudite + marićite appears between 500 and 700°C. The limit between these two fields corresponds to the maximum temperature that can be reached by alluaudites in granitic pegmatites, because marićite has never been observed in these geological environments. Because alluaudites are very sensitive to variations of oxygen fugacity, the field of hagendorfite, Na2MnFe2+Fe3+(PO4)3, has been positioned in the f(O2)–T diagram, and provides a tool that can be used to estimate the oxygen fugacity conditions that prevailed in granitic pegmatites during the crystallization of this phosphate. 相似文献
14.
Wenjun Yong E. Dachs A. C. Withers E. J. Essene 《Physics and Chemistry of Minerals》2007,34(2):121-127
A multi-anvil device was used to synthesize 24 mg of pure γ-Fe2SiO4 crystals at 8.5 GPa and 1,273 K. The low-temperature heat capacity (C
p) of γ-Fe2SiO4 was measured between 5 and 303 K using the heat capacity option of a physical properties measurement system. The measured
heat capacity data show a broad λ-transition at 11.8 K. The difference in the C
p between fayalite and γ-Fe2SiO4 is reduced as the temperature increases in the range of 50–300 K. The gap in C
p data between 300 and 350 K of γ-Fe2SiO4 is an impediment to calculation of a precise C
p equation above 298 K that can be used for phase equilibrium calculations at high temperatures and high pressures. The C
p and entropy of γ-Fe2SiO4 at standard temperature and pressure (S°298) are 131.1 ± 0.6 and 140.2 ± 0.4 J mol−1 K−1, respectively. The Gibbs free energy at standard pressure and temperature (ΔG°
f,298) is calculated to be −1,369.3 ± 2.7 J mol−1 based on the new entropy data. The phase boundary for the fayalite–γ-Fe2SiO4 transition at 298 K based on current thermodynamic data is located at 2.4 ± 0.6 GPa with a slope of 25.4 bars/K, consistent
with extrapolated results of previous experimental studies. 相似文献
15.
Fe-Ti oxide geothermometry: thermodynamic formulation and the estimation of intensive variables in silicic magmas 总被引：3，自引：0，他引：3
A new thermodynamic formulation of the Fe–Ti oxide geothermometer/oxygen barometer is developed. The method is based upon recently calibrated models for spinel solid solutions in the quinary system (Fe2+, Mg)(Al,Fe3+,Cr)2O4–(Fe2+, Mg)2TiO4 by Sack and Ghiorso, and rhombohedral oxides in the quaternary system (Fe2+,Mg,Mn)TiO3–Fe2O3 (this paper). The formulation is internally consistent with thermodynamic models for (Fe2+,Mg)-olivine and -orthopyroxene solid solutions and end-member thermodynamic properties tabulated by Berman. The constituent expressions account for compositional and temperature dependent cation ordering and reproduce miscibility gap features in all of the component binaries. The calibration does not account for the excess Gibbs energy resulting from compositional and temperature dependent magnetic ordering in either phase. This limits application of the method to assemblages that equilibrated at temperatures above 600° C. Practical implementation of the proposed geothermometer/oxygen barometer requires minimal use of projection algorthms in accommodating compositions of naturally occurring phases. The new formulation is applied to the estimation of temperature and oxygen fugacity in a wide variety of intermediate to silicic volcanic rocks. In combination with previous work on olivine and orthopyroxene thermodynamics, equilibration pressures are computed for a subset of these volcanics that contain the assemblage quartz, oxides and either ferromagnesian silicate. The calculated log10
f
O
2-T relations are reflected in coexisting ferromagnesian mineral assemblages. Volcanics with the lowest relative oxygen fugacity (log10
f
O
2) are characterized by the assemblage olivine-quartz, those with slightly higher log10
f
O
2 s, by the assemblage orthopyroxene-quartz. The sequence proceeds with the necessary phases biotite-feldspar, then hornblende-quartz-clinopyroxene, and finally at the highest log10
f
O
2 s, sphene-quartz-clinopyroxene. Quantitative analysis of these trends, utilizing thermodynamic data for the constituent phases, establishes that, in most cases, the T-log10
f
O
2value computed from the oxides is consistent with the compositions of coexisting silicate phases, indicating that phenocryst equilibrium was achieved prior to eruption. There is, however, considerable evidence of oxide-silicate disequilibrium in samples collected from more slowly cooled domes and obsidians. In addition, T-log10
f
O
2trends from volcanic rocks that contain biotite and orthopyroxene are interpreted to imply a condition of Fe2+–Mg exchange disequilibrium between orthopyroxene and coexisting ferromagnesian silicates and melt. It is suspected that many biotite-feldspar-quartz-orthopyroxene bearing low temperature volcanic rocks inherit orthopyroxene xenocrysts which crystallized earlier in the cooling history of the magma body.The problem is probably at least as complex as that of the feldspars... A.F. Buddington (1956) 相似文献
16.
Laser Ablation ICPMS study of trace element partitioning between plagioclase and basaltic melts: an experimental approach 总被引：3，自引：1，他引：3
Mario Aigner-Torres Jon Blundy Peter Ulmer Thomas Pettke 《Contributions to Mineralogy and Petrology》2007,153(6):647-667
Plagioclase-melt partition coefficients (D) for 34 trace elements at natural concentration levels were determined experimentally in a natural MORB composition at atmospheric
pressure using thin Pt-wire loops. Experiments were carried out at three temperatures (1,220, 1,200, and 1,180°C), and at
three different oxygen fugacities (fO2 = IW, QFM, air) in order to assess the effect of fO2 on the partitioning of elements with multiple valence (Fe, Eu, Cr). Run products were analyzed by laser-ablation ICP-MS.
Most trace element Ds increase slightly as temperature decreases, except for D
Zr, D
Fe, D
Eu and D
Cr that vary systematically with fO2. Applying the Lattice Strain Model to our data suggests the presence of Fe2+ entirely in the octahedral site at highly to moderate reducing conditions, while Fe3+ was assigned wholly to the tetrahedral site of the plagioclase structure. Furthermore, we provide a new quantitative framework
for understanding the partitioning behaviour of Eu, which occurs as both 2+ and 3+ cations, depending on fO2and confirm the greater compatibility of Eu2+, which has an ionic radius similar to Sr, relative to Eu3+ in plagioclase and the higher Eu2+/ Eu3+ under reducing conditions. For petrogenetic basaltic processes, a combined fractionation of Eu2+–Sr and Fe–Mg by plagioclase has considerable potential as an oxybarometer for natural magmatic rocks. 相似文献
17.
The full set of elastic constants for plagioclase end-member phase albite (NaAlSi3O8) is reported for the first time. Velocities of surface acoustic waves (both Rayleigh and pseudo-surface waves) were measured using impulsively stimulated light scattering on polished surfaces having six different orientations (three normal to the Cartesian axes and three lying on diagonals). Data were inverted and results tested using several non-linear optimization techniques. Compliance moduli determined under hydrostatic compression provided additional constraints and reduced covariance in the reported constants. The Cartesian coordinate system associated with the constants (using the
unit cell) has the y-axis parallel to the crystal b axis, the x-axis parallel to a* (perpendicular to b and c) and the z-axis consistent with a right-handed coordinate system. The values of the moduli C11, C12, C13, C14, C15, C16, C22, C23, C24, C25, C26, C33, C34, C35, C36, C44, C45, C46, C55, C56, C66 and their 2σ uncertainties (in parentheses) are, respectively, 69.1(0.6), 34.0(0.7), 30.8(0.5), 5.1(0.1), −2.4(0.1), −0.9(0.1), 183.5(2.7), 5.5(2.2), −3.9(0.5), −7.7(0.7), −5.8(0.7), 179.5(2.3), −8.7(0.4), 7.1(0.6), −9.8(0.6), 24.9(0.1), −2.4(0.1), −7.2(0.1), 26.8 (0.2), 0.5(0.1), 33.5(0.2). These constants differ significantly from the previously reported pseudo-monoclinic constants that were based on velocity measurements on polysynthetic twinned crystal aggregates. Differences are consistent with systematic errors in the earlier study associated with sparse data and the presence of cracks and other imperfections. 相似文献
18.
Lin-gun Liu Chien-chih Chen Chung-Cherng Lin Yi-jong Yang 《Physics and Chemistry of Minerals》2005,32(2):97-102
The elastic constants of natural single-crystal aragonite (CaCO3) have been measured by Brillouin spectroscopy at ambient conditions. The elastic constants C11, C22, C33, C44, C55, C66, C12, C13 and C23 are 171.1±1.0, 110.1±0.9, 98.4±1.2, 39.3±0.6, 24.2±0.4, 40.2±0.6, 60.3±1.0, 27.8±1.6 and 41.9±2.0 GPa, respectively, for aragonite. The linear compressibilities of the a-, b- and c-axis for aragonite at ambient conditions were derived from our measured data to be 3.0±0.2, 4.2±0.2 and 7.3±0.6×10–3 GPa–1, respectively. The aggregate bulk and shear moduli for aragonite using the Voigt-Reuss-Hill (VRH) scheme are thus calculated to be 68.9±1.4 and 35.8±0.2 GPa, respectively. The value of bulk modulus is in remarkable contrast to the literature value of 46.9 GPa measured almost a century ago. Our new datum, however, is closer to that derived from recent atomistic simulation and static compression studies. 相似文献
19.
The high-pressure stability of zoisite and phase relationships of zoisite-bearing assemblages 总被引：5，自引：0，他引：5
The fluid-absent reaction 12 zoisite = 3 lawsonite + 7 grossular + 8 kyanite + 1 coesite was experimentally reversed in the
model system CaO-Al2O3-SiO2-H2O (CASH) using a multi-anvil apparatus. The upper pressure stability limit for zoisite was found to extend to 5.0 GPa at 700 °C
and to 6.6 GPa at 950 °C. Additional experiments both in the H2O-SiO2-saturated and in the H2O-Al2O3-saturated portions of CASH provide further constraints on high pressure phase relationships of lawsonite, zoisite, grossular,
kyanite, coesite, and an aqueous fluid. Consistency of the present experiments with the H2O-saturated breakdown of lawsonite is demonstrated by thermodynamic analysis using linear programming techniques. Two sets
of data consistent with databases of Berman (1988) and Holland and Powell (1990) were retrieved combining experimental phase
relationships, calorimetric constraints, and recently measured elastic properties of solid phases. The best fits result in
G
f
,1,298
∘,zoisite=−6,499,400 J and S
1,298
∘,zoisite=302 J/K, and G
f
,1,298
∘,lawsonite=−4,514,600 J and S
1,298
∘,lawsonite=220 J/K for the dataset of Holland and Powell, and G
f
,1,298
∘,zoisite=−6,492,120 J and S
1,298
∘,zoisite=304 J/K, and G
f
,1,298
∘,lawsonite=−4,513,000 J and S
1,298
∘,lawsonite= 218 J/K for the dataset of Berman. Examples of the usage of zoisite as a geohygrometer and as a geobarometer in rocks metamorphosed
at eclogite facies conditions are worked, profiting from the thermodynamic properties retrieved here.
Received: 23 December 1996 / Accepted: 29 August 1997 相似文献
20.
Stable carbon isotope compositions and source rock geochemistry of the giant gas accumulations in the Ordos Basin, China 总被引：16，自引：0，他引：16
Jinxing Dai Jian Li Xia Luo Wenzheng Zhang Guoyi Hu Chenghua Ma Jianmin Guo Shouguo Ge 《Organic Geochemistry》2005,36(12):1617
Ordos Basin, the second largest sedimentary basin in China, contains enormous natural gas resources. Each of the four giant gas fields discovered so far in this basin (i.e., Sulige, Yulin, Wushenqi and Jingbian) has over 100 billion cubic meters (bcm) or 3.53 trillion cubic feet (tcf) of proven gas reserves. This study examines the stable carbon isotope data of 125 gas samples collected from the four giant gas fields in the Ordos Basin. Source rocks in the Upper Paleozoic coal measures are suggested by the generally high δ13C values of C1–C4 gaseous hydrocarbons in the gases from the Sulige, Yulin and Wushenqi gas fields. While the δ13CiC4 value is higher than that of the δ13CnC4, the dominant ranges for the δ13C1, δ13C2, and δ13C3 values in these Upper Paleozoic reservoired gases are −34 to −32‰, −27 to −23‰, and −25 to −24‰, respectively. The δ13C values of methane, benzene and toluene in gases from the Lower Paleozoic reservoirs of the Jingbian field indicate a significant contribution from humic source rocks, as they are similar to those in the Upper Paleozoic reservoirs of the Sulige, Yulin and Wushenqi gas fields. However, the wide variation and reversal in the δ13C1, δ13C2 and δ13C3 values in the Jinbian gases cannot be explained using a single source scenario, thus the gases were likely derived dominantly from the Carboniferous-Permian coal measures with some contribution from the carbonates in the Lower Permian Taiyuan Formation. The gas isotope data and extremely low total organic carbon contents (<0.2% TOC) suggest that the Ordovician Majiagou Formation carbonates are unlikely to be a significant gas source rock, thus almost all of the economic gas accumulations in the Ordos Basin were derived from Upper Paleozoic source rocks. 相似文献