Determination of gold, platinum and palladium in geological samples by inductively coupled plasma-mass spectrometry after separation and preconcentration with C-410 anion exchange resin

A method for the determination of Au, Pt and Pd in geological samples is described. Au, Pt and Pd can be separated and concentrated quantitatively by C-410 anion-exchange resin in the condition of 1.5 mol/L HCl with the adsorption rates of 91.2%, 100.0% and 95.7% respectively. No interference exists from coexisting elements except for Ge(IV), Cr(VI),Ti(IV) in inductively coupled plasma-mass spectrometry. The detection limits are 0.27 mug/L, 0.40 mug/L and 0.19 mug/L for Au, Pt and Pd respectively. The results of these elements in standard geological materials are in agreement with certified values with precision of 19.2% RSD for Au (n = 8), 28.1% RSD for Pt (n=8), and 15.6% RSD for Pd (n=8).

Tectonic geomorphology of the Ryukyu Trench-Arc-Backarc System: geological-geophysical exploration and mapping

OKINAWA TROUGH; BASIN

Recent 2000-year geological records of mud in the inner shelf of the East China Sea and their climatic implications

AXIS(14)C dating and grain-size analysis for Core DD2, located at the north of the Yangtze River-derived mud off the Zhejiang-Fujian coasts in the inner shelf of the East China Sea, provide us a high-resolution grain-size distribution curve varying with depth and time. Data in the upper mud layer of Core DD2 indicate that there are at least 9 abrupt grain-size increasing in recent 2000 years, with each corresponding very well with the low-temperature events in Chinese history, which might result from the periodical strengthening of the East Asian Winter Monsoon (EAWM), including the first-revealed maximum temperature lowering event at around 990 a BP. At the same time, the finer grain size section in Core DD2 agrees well with the Sui-Tang Warming Period (600-1000 a AD) defined previously by Zhu Kezhen, during which the climate had a warm, cold and warm fluctuation, with a dominated cooling period of 750-850 a AD. The Little Ice Age (LIA) can also be identified in the core. It starts around 1450 a AD and was followed by a subsequent cooling events at 1510, 1670 and 1840 a AD. Timing of these cold events revealed here still needs to be further verified owing to some current uncertainty of dating we used in this study.

Simultaneous determination of four phases of copper in geological samples

According to the specific property of respective phases in copper ore, optimal extracting condition for, tour phases of copper with five infusion was examined. After a series of experiments, for example, different weight Of sample, artificial mixed single phases sample and repeated experiment, the extracting ratio constants foil all phases in different infusions were obtained. The mathematical model which is on the basis of the stable extracting constant is established. Based on the purpose-control genetic algorithms, we only need to determine the total copper in the sample in different infusions under the selected condition, then the value of respective phase could be calculated. The-relative standard deviation are (%) free oxidized copper phase:1.5, conjunction oxidized copper phase: 17.6,, secondary copper sulfide phase: 1.9, primary copper sulfide phase : 2.7,total copper: 0.8.

Geological records of marine environmental changes in the southern Okinawa Trough

Indexes of sediment grain size, sedimentation rates, geochemical composition, heavy minerals, benthic foraminiferal fauna, indicator species of the Kuroshio Current, paleo-SST and carbonate dissolution of core E017 conformably suggest a great marine environmental change occurring at about 10.1-9.2 cal. kaBP in the southern Okinawa Trough, which may correspond to the strengthening of the Kuroshio Warm Current and re-entering the Okinawa Trough through the sea area off northeast Taiwan. The invasion of Kuroshio current has experienced a process of gradual strengthening and then weakening, and its intensity became more fluctuation during the last 5000 years. Compared to the transition of sediment grain size, geochemical composition and heavy minerals, the foraminiferal faunas show a 900-year lag, which may indicate that the invasion of Kuroshio Current and the consequent sea surface and deep-water environmental changes is a gradual process, and fauna has an obvious lag compared to environment altering. The carbonate dissolution of the Okinawa Trough has had an apparent strengthening since 9.2 cal. kaBP, and reached a maximum in the late 3000 years, which may be caused by the deep-water environmental changes due to the invasion of Kuroshio Current.

Petroleum geological framework and hydrocarbon potential in the Yellow Sea

Sedimentary basins in the Yellow Sea can be grouped tectonically into the North Yellow Sea Basin (NYSB), the northern basin of the South Yellow Sea (SYSNB) and the southern basin of the South Yellow Sea (SYSSB). The NYSB is connected to Anju Basin to the east. The SYSSB extends to Subei Basin to the west. The acoustic basement of basins in the North Yellow Sea and South Yellow Sea is disparate, having different stratigraphic evolution and oil accumulation features, even though they have been under the same stress regime since the Late Triassic. The acoustic basement of the NYSB features China-Korea Platform crystalline rocks, whereas those in the SYSNB and SYSSB are of the Paleozoic Yangtze Platform sedimentary layers or metamorphic rocks. Since the Late Mesozoic terrestrial strata in the eastern of the NYSB (West Korea Bay Basin) were discovered having industrial hydrocarbon accumulation, the oil potential in the Mesozoic strata in the west depression of the basin could be promising, although the petroleum exploration in the South Yellow Sea has made no break-through yet. New deep reflection data and several drilling wells have indicated the source rock of the Mesozoic in the basins of South Yellow Sea, and the Paleozoic platform marine facies in the SYSSB and Central Rise could be the other hosts of oil or natural gas. The Mesozoic hydrocarbon could be found in the Mesozoic of the foredeep basin in the SYSNB that bears potential hydrocarbon in thick Cretaceous strata, and so does the SYSSB where the same petroleum system exists to that of oil-bearing Subei Basin.

东海陆架盆地复杂天然气藏地震勘探判别技术研究

Exploration study proves that East sea shelf basin embeds abundant hydrocarbon resources. However, the exploration knowledge of this area is very low. Many problems in exploration are encountered here. One of them is that the gas reservoir of this area, with rapid lateral variation, is deeply buried. Correlation of Impendence between sandstone, gas sand and shale is very poor. Another problem is that the S/N ratio of the seismic data is very low and multiples are relatively productive which seriously affect reservoir identification. Resolution of the seismic data reflected from 2500-3000 meter is rather low, which seriously affects the application of hydrocarbon direct identification (HDI) technology. This research established a fine geological & geophysical model based on drilling、well logging、geology&seismic data of East sea Lishui area. A Q value extraction method from seismic data is proposed. With this method, Q value inversion from VSP data and seismic data is performed to determine the subsurface absorption of this area. Then wave propagation and absorption rule are in control. Field acquisition design can be directed. And at the same time, with the optimization of source system, the performance of high resolution seismic acquisition layout system is enhanced. So the firm foundation is ensured for east sea gas reservoir exploration. For solving the multiple and amplitude preserving problems during the seismic data processing, wave equation pre-stack amplitude preservation migration and wave equation feedback iteratively multiple attenuation technologies are developed. Amplitude preservation migration technology can preserve the amplitude of imaging condition and wave-field extrapolation. Multiple removing technology is independent of seismic source wavelet and velocity model, which avoiding the weakness of Delft method. Aiming at the complicated formation condition of the gas reservoir in this area, with dissecting typical hydrocarbon reservoir, a series of pertinent advanced gas reservoir seismic identification technologies such as petrophysical properties analyzing and seismic modeling technology、pre-stack/post-stack joint elastic inversion, attribute extraction technology based on seismic non-stationary signal theory and formation absorption characteristic and so on are studied and developed. Integrated analysis of pre-stack/post-stack seismic data, reservoir information, rock physics and attribute information is performed. And finally, a suit of gas reservoir identification technology is built, according to the geological and geophysical characteristics of this area. With developed innovative technologies, practical application and intergrated interpretation appraisal researches are carried out in Lishui 36-1.The validity of these technologies is tested and verified. Also the hydrocarbon charging possibility and position of those three east sea gas exploration targets are clearly pointed out.

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

岩性油藏储层预测的地震理论与方法研究

At present the main object of the exploration and development (E&D) of oil and gas is not the structural oil-gas pools but the subtle lithological oil-gas reservoir. Since the last 90's, the ratio of this kind of pools in newly-added oil reserves is becoming larger and larger, so is the ratio in the eastern oilfields. The third oil-gas resource evaluation indicates the main exploration object of Jiyang depression is the lithological oil-gas pools in future. However, lack of effective methods that are applied to search for this kind of pool makes E&D difficult and the cost high. In view of the urgent demand of E&D, in this paper we deeply study and analyze the theory and application in which the seismic attributes are used to predict and describe lithological oil-gas reservoirs. The great results are obtained by making full use of abundant physics and reservoir information as well as the remarkable lateral continuity involved in seismic data in combination with well logging, drilling-well and geology. ①Based on a great deal of research and different geological features of Shengli oilfield, the great progresses are made some theories and methods of seismic reservoir prediction and description. Three kinds of extrapolation near well seismic wavelet methods-inverse distance interpolation, phase interpolation and pseudo well reflectivity-are improved; particularly, in sparse well area the method of getting pseudo well reflectivity is given by the application of the wavelet theory. The formulae for seismic attributes and coherent volumes are derived theoretically, and the optimal method of seismic attributes and improved algorithms of picking up coherent data volumes are put forward. The method of making sequence analysis on seismic data is put forward and derived in which the wavelet transform is used to analyze not only qualitatively but also quantitatively seismic characteristics of reservoirs.② According to geologic model and seismic forward simulation, from macro to micro, the method of pre- and post-stack data synthetic analysis and application is put forward using seismic in close combination with geology; particularly, based on making full use of post-stack seismic data, "green food"-pre-stack seismic data is as possible as utilized. ③ In this paper, the formative law and distributing characteristic of lithologic oil-gas pools of the Tertiary in Jiyang depression, the knowledge of geological geophysics and the feasibility of all sorts of seismic methods, and the applied knowledge of seismic data and the geophysical mechanism of oil-gas reservoirs are studied. Therefore a series of perfect seismic technique and software are completed that fit to E&D of different categories of lithologic oil-gas reservoirs. ④ This achievement is different from other new seismic methods that are put forward in the recent years, that is multi-wave multi-component seismic, cross hole seismic, vertical seismic, and time-lapse seismic etc. that need the reacquisition of seismic data to predict and describe the oil-gas reservoir. The method in this paper is based on the conventional 2D/3D seismic data, so the cost falls sharply. ⑤ In recent years this technique that predict and describe lithologic oil-gas reservoirs by seismic information has been applied in E&D of lithologic oil-gas reservoirs on glutenite fans in abrupt slop and turbidite fans in front of abrup slop, slump turbidite fans in front of delta, turbidite fans with channel in low slope and channel sanbody, and a encouraging geologic result has been gained. This achievement indicates that the application of seismic information is one of the most effective ways in solving the present problem of E&D. This technique is significant in the application and popularization, and positive on increasing reserves and raising production as well as stable development in Shengli oilfield. And it will be directive to E&D of some similar reservoirs

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.

地震反射系数成像方法研究

Reflectivity sequences extraction is a key part of impedance inversion in seismic exploration. Although many valid inversion methods exist, with crosswell seismic data, the frequency brand of seismic data can not be broadened to satisfy the practical need. It is an urgent problem to be solved. Pre-stack depth migration which developed in these years becomes more and more robust in the exploration. It is a powerful technology of imaging to the geological object with complex structure and its final result is reflectivity imaging. Based on the reflectivity imaging of crosswell seismic data and wave equation, this paper completed such works as follows: Completes the workflow of blind deconvolution, Cauchy criteria is used to regulate the inversion(sparse inversion). Also the precondition conjugate gradient(PCG) based on Krylov subspace is combined with to decrease the computation, improves the speed, and the transition matrix is not necessary anymore be positive and symmetric. This method is used to the high frequency recovery of crosswell seismic section and the result is satisfactory. Application of rotation transform and viterbi algorithm in the preprocess of equation prestack depth migration. In equation prestack depth migration, the grid of seismic dataset is required to be regular. Due to the influence of complex terrain and fold, the acquisition geometry sometimes becomes irregular. At the same time, to avoid the aliasing produced by the sparse sample along the on-line, interpolation should be done between tracks. In this paper, I use the rotation transform to make on-line run parallel with the coordinate, and also use the viterbi algorithm to complete the automatic picking of events, the result is satisfactory. 1. Imaging is a key part of pre-stack depth migration besides extrapolation. Imaging condition can influence the final result of reflectivity sequences imaging greatly however accurate the extrapolation operator is. The author does migration of Marmousi under different imaging conditions. And analyzes these methods according to the results. The results of computation show that imaging condition which stabilize source wave field and the least-squares estimation imaging condition in this paper are better than the conventional correlation imaging condition. The traditional pattern of "distributed computing and mass decision" is wisely adopted in the field of seismic data processing and becoming an obstacle of the promoting of the enterprise management level. Thus at the end of this paper, a systemic solution scheme, which employs the mode of "distributed computing - centralized storage - instant release", is brought forward, based on the combination of C/S and B/S release models. The architecture of the solution, the corresponding web technology and the client software are introduced. The application shows that the validity of this scheme.