滇西一些新生代飞来峰的构造演化

在印度和欧亚大陆相互作用下,滇西三江造山带在新生代中期发生了强烈的东西向缩短,地壳物质的大规模推覆作用形成众多的飞来峰。野外调查发现:这些飞来峰中有的在就位后没有再发生过变形,如分布在兰坪—思茅盆地中段的无量山东北山麓一带的一系列上三叠统灰岩飞来峰;有的飞来峰在就位后再次遭受到强烈的变形,如分布在兰坪盆地金顶一带的上三叠统灰岩飞来峰,在渐新世时与下伏岩层一同卷入了强烈的褶皱变形;分布在大理地区鹤庆盆地西侧的中三叠统灰岩飞来峰因受到后期东西向伸展构造的破坏,已被改造成分选极差的重力滑塌沉积。这些古近纪飞来峰由于在就位后再次发生了变形,所以不应称之为飞来峰,应称为异地岩或者是外来岩块。

Seismic AVAZ inversion for orthorhombic shale reservoirs in the Longmaxi area, Sichuan

Seismic AVAZ inversion method based on an orthorhombic model can be used to invert anisotropy parameters of the Longmaxi shale gas reservoir in the Sichuan Basin..As traditional seismic inversion workfl ow does not suffi ciently consider the infl uence of fracture orientation,we predict fracture orientation using the method based on the Fourier series to correct pre-stacked azimuth gathers to guarantee the accuracy of input data,and then conduct seismic AVAZ inversion based on the VTI constraints and Bayesian framework to predict anisotropy parameters of the shale gas reservoir in the study area.We further analyze the rock physical relation between anisotropy parameters and fracture compliance and mineral content for quantitative interpretation of seismic inversion results.Research results reveal that the inverted anisotropy parameters are related to P-and S-wave respectively,and thus can be used to distinguish the effect of fracture and fl uids by the joint interpretation.Meanwhile high values of anisotropy parameters correspond to high values of fracture compliance,so the anisotropy parameters can refl ect the development of fractures in reservoir.There is two sets of data from different sources,including the content of brittle mineral quartz obtained from well data and the anisotropy parameters inverted from seismic data,also show the positive correlation.This further indicates high content of brittle mineral makes fractures developing in shale reservoir and enhances seismic anisotropy of the shale reservoir.The inversion results demonstrate the characterization of fractures and brittleness for the Longmaxi shale gas reservoir in the Sichuan Basin.

Heterogenous zircons from Mesozoic igneous rocks in the North China Craton and their tectonic implications

Based on a synthesis of zircon inheritance from Mesozoic igneous rocks in the eastern portion of the North China Craton （NCC） , it is proposed that inherited zircons with Neoproterozoie ages identified from these rocks are of a heterogenous derivation from the Yangtze/South China block, rather than from the NCC itself. The mechanism that introduces these zircons incorporated into the NCC is likely by tectonic underplating during the Triassic continental subduction of the Yangtze block beneath the NCC. Tectonic addition of abundant crustal materials represented by the heterogenous zircons into the NCC, probably along Moho or weak interfaces within the NCC＇s crust, led to the crustal thickening in the NCC. These heterogenous materials, either as （partial） source rocks or as contaminants of the magmas generated during an extension environment following the crustal thickening, were reworked and therefore have significant contribution to petrogenesis of the Mesozoic igneous rocks. The crustal thickening resulted from the tectonic underplating, as indicated by the distribution extent of the heterogenous zircons, is spatially similar to that of the lithospheric thinning, with both mainly occurring in the eastern segment of the NCC. This is probably suggestive of an intrinsic relation between the thickening and thinning events during the Mesozoic evolution of the NCC.

Continental crust formation at arcs, the arclogite ‘‘delamination''cycle, and one origin for fertile melting anomalies in the mantle

The total magmatic output in modern arcs,where continental crust is now being formed, is believed to derive from melting of the mantle wedge and is largely basaltic. Globally averaged continental crust, however, has an andesitic bulk composition and is hence too silicic to have been derived directly from the mantle. It is well known that one way this imbalance can be reconciled is if the parental basalt differentiates into a mafic garnet pyroxenitic residue/cumulate(‘‘arclogite'') and a complementary silicic melt, the former foundering or delaminating into the mantle due to its high densities and the latter remaining as the crust.Using the Sierra Nevada batholith in California as a case study, the composition of mature continental arc crust is shown in part to be the product of a cyclic process beginning with the growth of an arclogite layer followed by delamination of this layer and post-delamination basaltic underplating/recharge into what remains of the continental crust.A model is presented, wherein continuous arc magmatism and production of arclogites in continental arcs are periodically punctuated by a delamination event and an associated magmatic pulse every *10–30 My. The recycling flux of arclogites is estimated to be *5 %–20 % that of oceanic crust recycling by subduction. Delaminated arclogites have the necessary trace-element compositions to yield time-integrated isotopic compositions similar to those inferred toexist as reservoirs in the mantle. Because of their low melting temperatures, such pyroxenites may be preferentially melted, possibly forming a component of some hotspot magmas.

Numerical modeling of the anomalous post-rift subsidence in the Baiyun Sag,Pearl River Mouth Basin

The Baiyun Sag is the deepest sag in the Pearl River Mouth Basin in northern continental margin of South China Sea,with the maximum sediment thickness over 12.5 km above the basement including >6.5 km sediments above the 30 Ma breakup unconformity.According to the theoretical models for the rifted basins,the post-rift subsidence is driven solely by the thermal contraction and can be calculated as the function of the lithospheric stretching factor.A method combining the forward modeling and reverse backstripping was designed to estimate lithospheric stretching factor.Using the 2D forward modeling based on the flexural cantilever model,we simulated the multi-rifting process of the Baiyun Sag with constrain of the backstripped profiles.By doing this the lithospheric stretching factor was obtained,and then the theoretical post-rift subsidence was calculated.The calculated theoretical subsidence was much smaller than the observed subsidence given by backstripping.Along the 1530 line in the Baiyun Sag,the anomalous post-rift subsidence is over 2 km in the sag center,and varies slightly to the north and south edges of the sag.This suggests that the anomalous post-rift subsidence continues beyond the sag both in the continental shelf to the north and in the continental slope to the south.The sensitivity tests in the forward modeling process indicate that only the use of low-angle faults (≤13°) can we simulate the shape of the backstripped profile.

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/Tm, where T is temperature and Tm is the melting （solidus） temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system, （Mg,Fe）2SiO4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite （Fe2SiO4）, the triple point of fayalite-Fe2SiO4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/Tm of olivine with depth allows us to compare the strength of the up- per mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/Tm of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/Tm of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/Tm 〉 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/Tm is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/Tm of olivine need further investigation. Below 6.4 GPa （〈200 kin）, T/TIn of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10-12-10-15 s-1 and low stress in the upper mantle, the [100]（010） slip system （A-type fabric） becomes dominant when T/TIn〉 0.55-0.60. When T/Tm〈 0.55-0.60, [001] slip is easier and low T/Tm favors the operation of [001]（100） slip system （C-type fabric）. This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/Tm. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experi- ments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.

Platinum-group element geochemistry of Cenozoic basalts from the North China Craton: Implications for mantle heterogeneity

Forty-two Cenozoic （mostly Miocene） basalt samples from Jining, Chifeng, Fansi, Xiyang, and Zuoquan areas of the North China Craton （the NCC basalts hereafter） were analyzed for platinum-group elements （PGE, including Os, It, Ru, Rh, Pt, and Pd）. Most of them are alkaline basalts and tholeiites and all of them display little crustal contamination. The total PGE contents of the NCC basalts vary from 0.1 to 0.9 ppb, much lower than those of the primitive mantle values of 23.5 ppb. Primitive man- tie-normalized PGE patterns of these basalts define positive slopes and Pd/Ir ratios vary from 1.2 to 25. In terms of both PGE contents and Pd/Ir ratios, they are quite similar to the mid-ocean ridge basalts. There are no obvious negative correlations be- tween PGE vs. MgO, Ni, and Cu in the NCC basalts, indicating that fractional crystallization of olivine, pyroxene, and/or sul- fides during magmatic process cannot be the controlling factor for the observed PGE variation. The observed Pd/Ir variations of the NCC basalts require involvement of non-chondritic heterogeneous mantle sources. Based on Sr-Nd-Pb-Hf isotopic sys- tematics and incompatible-element signatures, a mixing of partial melts from both asthenospheric peridotites and enclosed mantle eclogites at the top of asthenosphere was proposed for the origin of these NCC basalts. The lenses of eclogites are de- rived from upwelling of recycled continental crust during the westward subduction of the Pacific plate from the -600 km dis- continuity zone. The PGE geochemistry of these basalts provides independent evidence to support this conclusion and the ob- served Pd/Ir variations may reflect variations in proportions of tapped peridotitic and eclogitic melts.

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.