Tang X.-D.,CAS Institute of Geology and Geophysics |
Tang X.-D.,University of Chinese Academy of Sciences |
Liu H.,CAS Institute of Geology and Geophysics |
Zhang H.,MLR Key Laboratory of Marine Mineral Resources
Chinese Journal of Geophysics (Acta Geophysica Sinica) | Year: 2015
Full-waveform inversion (FWI) can supply accurate velocity model for seismic imaging and is becoming a very important technology in seismic exploration. The frequency domain modeling is the base of frequency domain FWI and can be used for each iteration. In this case, accurate and fast frequency domain modeling is very necessary. The traditional frequency domain modeling methods can only be used in square grid and the memory and computation requirements are huge. We present a high efficient and adaptable frequency domain modeling technique in this paper. We design a 17-point scheme to perform the frequency domain modeling and use the average-derivative method to overcome the disadvantage that the traditional frequency domain modeling method can not be used in rectangular grid. Our operator can not only be applied in the condition of square grid but also rectangular grid, which is very useful in the real world situation because the rectangular grid is very popular in field work. Therefore, our operator is of greater applicability than the traditional methods. In order to reduce the dispersion and the number of grid points per shortest wavelength, we combine the mass acceleration term into the point calculation and get the optimal coefficients through simulated annealing algorithm. As the program Matlab provides the large scale optimization method, we use fmincon in this program to further optimize the coefficients. At last, the number of grid points per shortest wavelength is reduced to 2.5 comparing to 4 in 9-point method which is the most common frequency domain modeling method with phase velocity errors no larger than 1%. Therefore, the requirements of calculation amount and memory are reduced greatly. There are two main ways in solving Helmholtz equation which is the algorithm of frequency domain acoustic modeling, including direct method and iterative method. The direct method is LU decomposition while the matrix of LU decomposition in our algorithm usually is huge in 2D data and is difficult to solve using current computer. Therefore, iterative method is an option. However, in iterative method the convergent rate usually is very low, which leads to huge computation. We use the incomplete LU decomposition as the precondition and accelerate the speed of convergence. For the same test, the traditional method needs more than 20000 times to converge to a stable value while our method needs only 6 times, which reduces the amount of calculation greatly. In order to accelerate the process of frequency domain modeling, we apply the Graphics Processing Unit (GPU) because GPU has powerful computing capability. In the GPU, fair use of video memory is very crucial and can lead to higher efficiency. Using the GPU, we obtain twice speed of the traditional CPU method in modeling. There are two main methods in frequency domain acoustic modeling, including 9-point method and 25-point method. We compare our method with those of two methods in accuracy and efficiency. We use the same velocity model and the dx=dz=25 m. The dominant frequency of the Ricker wavelet is 20 Hz. In order to remove all disturbances, we do not apply any boundary condition. Because our method only needs 2.4 grids, to reach the same accuracy our operator can use relatively coarse mesh and it leads to less computing work. Meanwhile, the storage is reduced, too. For example, if the 9-point method needs nx×ny grids, the nonzero elements of coefficient matrix is 9×nx×nx×nz×nz, the 25-point method is 9.76×nx×nx×nz×nz and our method is 6.12×nx×nx×nz×nz. Therefore, our method produces a great savings of memory. In order to test our method in the complex velocity model, we applied it in the Marmousi model the sampling intervals of which are 12.5 m and 4 m. The grid of Marmousi is rectangular and our method can be used in this kind of model while 9-point and 25-point can not. We also give the result of time domain method which is second-order in time and twelfth-order in space. Comparing the result from time-domain method with the result from our technique, we find that there is no significant difference between them. Besides, the boundary condition is very important in modeling, we derive the perfectly matched layer equation based on our 17-point scheme and verify it in the Marmousi model. We propose a 17-point scheme which has the advantage of wide application and high accuracy to perform frequency domain acoustic modeling. It overcomes the disadvantage of traditional methods which can only be used in the square grid. After optimization, the number of grid points per shortest wavelength is reduced to 2.5, which reduces huge memory requirement and computing time. Because of the slow convergence, we apply the incomplete LU decomposition as the precondition in solving Helmholtz equation and accelerate the speed of convergence. As the GPU has the powerful computing capability, we use it to accelerate the process of modeling and save much time on computing. ©, 2015, Science Press. All right reserved.
Sha Z.,MLR Key Laboratory of Marine Mineral Resources |
Sha Z.,Wuhan University |
Liang J.,MLR Key Laboratory of Marine Mineral Resources |
Zhang G.,MLR Key Laboratory of Marine Mineral Resources |
And 5 more authors.
Marine Geology | Year: 2015
The second gas hydrate expedition of the Guangzhou Marine Geological Survey drilled 13 sites and recovered a large amount of gas hydrate in the east of the Pearl River Mouth basin of the South China Sea. In this study, we examine three of these sites (sites 05, 08, and 16) using multichannel seismic data, LWD logging data, and multi-beam echo sounder data. The geophysical data suggest and core samples confirm that gas hydrates occur near the BSR at all three sites, whereas near-surface gas hydrates exist at sites 08 and 16 with high fluid flux.We present a conceptual model for the formation and accumulation of gas hydrates in our study area, controlled by the activity of dissolved and free gas-rich fluid. Interpreted abundant gas sources below the BSR contribute to the generation of overpressure, resulting in the movement of gas-rich fluid. Dissolved gas-rich fluid forms pore filling hydrates near the BSR, and the relatively deep hydrate layer thickens over time by making its way through the growing hydrate layer and adding to the top, if the fluid is supplied sufficiently. Free gas-rich fluid is more favorable to generate near-surface nodules or massive hydrates, which could be stable in natural conditions. As sediments deposit, the relatively shallow hydrates are buried deeper and keep their morphologies. Furthermore, gas hydrate formation and stability imply that gas hydrates may be increasingly considered as a potential energy resource. © 2015 Elsevier B.V.
Wang F.L.,MLR Key Laboratory of Marine Mineral Resources |
Wang F.L.,Sun Yat Sen University |
Zhao T.P.,CAS Guangzhou Institute of Geochemistry |
Wang Y.,CAS Guangzhou Institute of Geochemistry
Yanshi Xuebao/Acta Petrologica Sinica | Year: 2015
Syenitic dikes are widespread in the Panxi district and some of them host Nb and Ta deposits. The mineralized dikes intruded into gabbro intrusions and are spatially and temporally associated with basalts and felsic plutons (granite and syenite plutons) in the Panxi district, which are part of the ∼260Ma Emeishan large igneous province. These dikes are recently dated to be 258 ∼ 256Ma by LA-ICP-MS zircon U-Pb dating technique, contemporaneous with the gabbro and felsic plutons. The Sr-Nd isotopes of the mineralized syenitic dikes, barren syenitic dikes and syenitic plutons of Luku and Baicao Nb-Ta deposits in the Panxi district show similar and uniform initial εNd(t) values (-0.2 to +0.2 for mineralized syenitic dikes, -0.3 to +0.7 for barren syenitic dikes and -0.3 to +0.4 for syenitic plutons) but they have highly variable (87Sr/86Sr)i ratios (0.7049 to 0.7091 for the mineralized syenitic dikes and 0.7044 ∼ 0.7064 for barren syenitic dikes and 0.7032 to 0.7090 for barren syenitic plutons), comparable to the contemporaneous Hongge mafic/ultramafic layered intrusion εNd(t) = -2.7 to +1.0 and (86Sr/87Sr)i = 0.7058 to 0.7064) in the same region, suggesting that these lithologies may have evolved from a common basaltic magmas. According to the whole rock geochemistry and Sr-Nd isotopic characters, We consider that the syenitic plutons and dikes may be the products of low degree (5% ∼ 10%) partial remelting of underplated gabbro cumulates.
Wang F.,MLR Key Laboratory of Marine Mineral Resources |
Wang F.,CAS Guangzhou Institute of Geochemistry |
Wang F.,Sun Yat Sen University |
He G.,MLR Key Laboratory of Marine Mineral Resources |
And 3 more authors.
Yanshi Xuebao/Acta Petrologica Sinica | Year: 2016
The deep-sea sediments, especially the zeolite clays, are rich in rare earth elements and yttrium (REY, REE+Y) in the Pacific Ocean. The main components of the sediments include clay, zeolite, fish debris, micronodules and biological residues. However, the carriers of REY are under debate. We have conducted mineralogical investigations, and in-situ major and trace elemental analyses of constituent minerals of the REY-rich zeolite clay from the western and central areas in the Pacific Ocean in order to clarify the host-phase of REY in the mud. The EPMA and LA-ICP-MS data show that the fish debris are mainly composed of apatite with the highest REY contents of 2497 × 10-6 ∼ 18623 × 10-6, whereas the micronodules and zeolite have lower REY contents of 246 × 10-6 ∼ 333 × 10-6 and 29. 36 × 10-6 ∼ 256 × 10-6 We calculate that the apatites occupy above 90%of the total REY budget in the 63 ∼ 250μm fraction. However, the total REY value of apatites in the 63 ∼ 250μm fraction is still insignificant compared with the bulk sediment. The geochemical analyses of whole-rock and silt fraction (<63μm) and X-ray diffraction analysis show that REY mainly concentrates in the silt fraction and ∼ 70%of the whole REY budget concentartes in apatites in both samples. Furthermore, REY contents of the clay fraction (<2μm) in the two samples are 530 × 10-6 and 631 × 10-6 respectively and make up 2%-5%of the total REY. Therefore, we conclude that the main REY host phase in the zeolite clay is apatite in the western and central Pacific Ocean.
Li X.,Wuhan University |
Zhu P.,Wuhan University |
Kusky T.M.,Wuhan University |
Kusky T.M.,China Three Gorges University |
And 4 more authors.
Tectonophysics | Year: 2015
A review and comparison of the tectonic history and lithospheric structure of the Yangtze craton and North China craton (NCC) are made to test whether the Yangtze craton has lost part of its lithospheric root like the NCC has, and if so, is there a common driving tectonic force for root-loss. The results show that many similarities and some significant differences exist between these two cratons. The similarities include the following: (1) The eastern parts of the NCC and Yangtze craton have both been destroyed. The destroyed areas are bounded on the west by the North-South Gravity Lineament. (2) Fast velocity bodies located above depths between 200 and 400 km are detected beneath the NCC and Yangtze craton, implying that delamination may be an important mechanism of cratonic destruction for these two cratons. (3) The destruction of the NCC and Yangtze craton can be divided into several stages. The differences mainly lie in the tectonic evolution during Neoproterozoic-early Mesozoic. The tectonothermal events which are beneficial for the lithospheric weakening of the NCC dominantly affected the marginal regions. However, for the Yangtze craton, they are not only limited to the marginal regions, but also extend to the eastern and middle parts of the craton. This may be a key factor why the Yangtze craton shows more continuous slow velocity anomalies.We also make a study on the time and dynamic mechanism of the destruction of the Yangtze craton. The results show that the initial destruction started in the Neoproterozoic and may be caused by the aggregation and breakup of Rodinia. The large scale destruction occurred mostly in the Mesozoic (190-90. Ma) with a peak time of 130-120. Ma. The subduction and rollback of the Paleo-Pacific (Izanagi) plate may play a vital role in the destruction of the Yangtze craton. © 2015 Elsevier B.V.
Xi L.,MLR Key Laboratory of Marine Mineral Resources |
Li H.,South University of Science and Technology of China |
Xia Y.,Changjiang River Scientific Research Institute |
Qu W.,Water Resources University
Environmental Earth Sciences | Year: 2016
A mangrove transect (M–M′) and a bald beach transect (B–B′) were selected in Dongzhaigang National Nature Reserve, China for conducting a comparative study about groundwater chemical characteristics. It was found that the average concentration of Cl− in groundwater sampled from M–M′ was much lower than that of B–B′. However, concentrations of heavy metals Ba, Mn, Pb, Sb in groundwater of M–M′ were significantly higher than those of B–B′. The average concentrations of Ba, Mn, Pb, Sb were 0.13, 0.83, 0.25 and 0.32 mg/L in M–M′, respectively, and 0.09, 0.27, 0.19 and 0.21 mg/L in B–B′, respectively. Both transects were affected by the human activities such as aquacultures. The high concentration of Sb in the groundwater of transect M–M′ may be mainly due to high concentration of Sb in the Yanfeng River. The higher concentrations of Ba, Mn, Pb in groundwater in M–M′ might be attributed to the accumulating and enriching effects of the mangrove systems on heavy metals. The mechanism of accumulating and enriching need to be further investigated in the future. © 2016, Springer-Verlag Berlin Heidelberg.
Zhang G.-X.,MLR Key Laboratory of Marine Mineral Resources |
Xu H.-N.,MLR Key Laboratory of Marine Mineral Resources |
Liu X.-W.,Beijing University of Technology |
Zhang M.,MLR Key Laboratory of Marine Mineral Resources |
And 4 more authors.
Chinese Journal of Geophysics (Acta Geophysica Sinica) | Year: 2014
Based on integrated data acquired with 3-D seismic survey and OBS, both the images of subsurface formation and the information of P/S-wave velocities are obtained, which help to make certain the variation of P/S-wave velocity within the gas-hydrate bearing sediments in the study area. Comparison with the results of P/S-wave velocity calculated from travel time inversion indicates that some sediments have relatively high P/S-wave velocity within the depth of 500 m below seafloor. The abrupt rising of P-wave velocity means the presence of gas hydrate and the high velocity interval corresponds to the hydrate stable zone apparently, but the P-wave velocity decreasing as well slight S-wave velocity decreasing is very likely related to the occurrence of free gas. The presumed occurrence of free gas in the area does not agree with the viewpoints deduced from drilling and coring results within limited depth in this area in 2007.
Zhao Z.,MLR Key Laboratory of Marine Mineral Resources |
He Z.,Chengdu University of Technology |
Wan X.,MLR Key Laboratory of Marine Mineral Resources |
Shuai Q.,MLR Key Laboratory of Marine Mineral Resources
Xinan Shiyou Daxue Xuebao/Journal of Southwest Petroleum University | Year: 2016
The ability to identify gas-bearing and water-bearing sands of nine fluid identification factors has been compared, and the results show that the High-Sensitivity-Fluid-Identification-Factor has a strong ability to identify. The neural network and the Principal-Component-Analysis-neural-network are applied to high-quality 3-D data of HSFIF to perform the waveform analysis during the gas bearing interval in study area L of basin-S and good mapping effect has been achieved. The facies maps were analyzed and compared with the logging' interpretation. It proves the application of the PCA-neural network method can greatly reduce the difficulty of seismic facies interpretation of the map. In this paper, the application range of waveform classification is extended from seismic-sedimentary-facies analysis and reservoir prediction to the analysis and processing of the fluid factors in the target layers, thereby we can predict the fluid in layers. It is a new kind of valuable complement to fluid identification and prediction by using fluid factor only in profile and slice. It is the first time that the waveform classification techniques has been applied to fluid prediction. We believe that the method along with the results of other explanations has a guiding significance to reduce exploration risk and enhance drilling success rate. © 2016, Science Press. All right reserved.
Chen F.,MLR Key Laboratory of Marine Mineral Resources |
Hu Y.,Shanghai Ocean University |
Feng D.,CAS South China Sea Institute of Oceanology |
Zhang X.,MLR Key Laboratory of Marine Mineral Resources |
And 4 more authors.
Chemical Geology | Year: 2016
The application of molybdenum (Mo) enrichment to sediments can potentially provide unique constraints on the methane seepage dynamics at continental margins. We report herein elemental and isotopic measurements for authigenic carbonates and sediment samples in drill cores (~ 95 m long; Site GMGS2-08 at ~ 800 m water depth) from gas hydrate-bearing sediments in the northern South China Sea. Five intervals that were impacted by methane seepages in the sediment profile were identified by the presence of strongly 13C-depleted carbonates (δ13C values as low as − 56.8‰). The carbonates often cement chemosynthetic bivalve shell fragments and have anomalous δ18O positive values. These observations suggest that methane seepage was intense during those time intervals and that carbonate precipitation occurred in the uppermost sediments, probably due to gas hydrate dissociation. Interestingly, the five intervals were all characterized by strong Mo enrichments. We thus hypothesize that the presence of strong Mo enrichments is a good indicator of past episodes of methane seepage. Such an environment results in sulfidic conditions that occur in a narrow zone in proximity to seawater close to the seafloor, leading to resultant Mo enrichments in associated sediments. In contrast, environments with low methane flux would lead to the distribution of Mo throughout a wide sulfidic zone and, consequently, an absence of Mo enrichment in the sediments. The insights into the environment for Mo enrichments at seeps have significant implications for tracing past methane seepages, their intensities, and even possible occurrences of gas hydrate dissociation. © 2016 Elsevier B.V.
Gai S.H.,China University of Petroleum - Beijing |
Liu H.Q.,China University of Petroleum - Beijing |
He S.L.,China University of Petroleum - Beijing |
Lei G.,China University of Petroleum - Beijing |
And 4 more authors.
Geosystem Engineering | Year: 2016
Stress-dependent permeability and porosity of porous media is very important for the development of oil reservoirs, especially in low-permeability reservoirs. Porous structure can change in response to an increase in the effective stress to induce the decline of permeability and porosity, which ultimately affects the fluid flow property. In this study, a capillary module, considering the arrangement and the deformation of particles, has been established to express the stress sensitivity of porous media based on the particle packing model. The predictions of variation of permeability by the proposed model have been validated by comparing with the experimental data. The predictions by the current model show the same variation trend with the experimental data and are consistent with the compression experiments. The proposed variation of permeability with the increase in effective stress depends on the factors such as the arrangement of particles (contact angle), the deformation of particles (elastic module and Poisson’s ratio) and effective stress. The effects of these parameters on variation of permeability are discussed in detail. © 2015 The Korean Society of Mineral and Energy Resources Engineers (KSMER).