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

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.

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

滩海过渡带地震资料拼接处理方法研究

The technology of seismic acquirement has been more and more intact in both land area and marine area, however, since the acquirement equipment can not handle the tasks in this both two kinds of prospects at the same time, it is very hard to acquire in the prospect of seashore area, and also the quality of the raw data is poor, the related work in seashore area has been a tough task. With the development of the level of exploration, acquired equipment and technology, a certain series of processing methods has been built up, which enhanced the resolving power, the ratio of signal to noise and the quality of imaging, these methods could meet most of the requirements of exploration and the development. What is more, more and more nowadays prospects are composed by two or more 3D sub-prospect, the differences between frequencies, energy levels and the ratios of signal to noise are definitely obvious, we can not perform finely construct analysis and interpretation of reservoir without solving these problems in matching different raw data. As for the prospect that composed by both seashore and marine ones, through the analysis on the characteristics on frequency, amplitude, polarity and multiples of land prospect and marine prospect, the methods on how to eliminate those differences mentioned above were researched, and this paper gives a series of methods includes the computing of matching factor, recognizing of time difference, recognizing of phase difference, matching frequency, collapsing of abnormal energy, frequency compensating and so on, to fulfill the attenuations of energy difference, frequency difference, time difference and phase difference, and what is more, improve the innovative methods: includes the quantities method to synthetically and automatically recognize time difference and phase difference on the basis of signal match technology and according to statistic rule. And distinct the differences on phase and frequency characteristics of received raw data between velocity seismometer and piezoelectricity seismometer, and deduces that there are 90 degrees gap between data from these two seismometers, thirdly, synthetically apply zone abnormal frequency processing, surface consistent abnormal frequency compensate, residual abnormal frequency compensate, stack time equilibrium and optimize the processing seis-flow, deal with the would-be problems, such as low signal and noise rate, low resolution, poor imaging and poor state of consistent. Gain excellent effect in Da-gang oil field.

基于波动方程的层间多次波预测方法研究

In exploration geophysics，velocity analysis and migration methods except reverse time migration are based on ray theory or one-way wave-equation. So multiples are regarded as noise and required to be attenuated. It is very important to attenuate multiples for structure imaging, amplitude preserving migration. So it is an interesting research in theory and application about how to predict and attenuate internal multiples effectively. There are two methods based on wave-equation to predict internal multiples for pre-stack data. One is common focus point method. Another is inverse scattering series method. After comparison of the two methods, we found that there are four problems in common focus point method: 1. dependence of velocity model; 2. only internal multiples related to a layer can be predicted every time; 3. computing procedure is complex; 4. it is difficult to apply it in complex media. In order to overcome these problems, we adopt inverse scattering series method. However, inverse scattering series method also has some problems: 1. computing cost is high; 2. it is difficult to predict internal multiples in the far offset; 3. it is not able to predict internal multiples in complex media. Among those problems, high computing cost is the biggest barrier in field seismic processing. So I present 1D and 1.5D improved algorithms for reducing computing time. In addition, I proposed a new algorithm to solve the problem which exists in subtraction, especially for surface related to multiples. The creative results of my research are following: 1. derived an improved inverse scattering series prediction algorithm for 1D. The algorithm has very high computing efficiency. It is faster than old algorithm about twelve times in theory and faster about eighty times for lower spatial complexity in practice; 2. derived an improved inverse scattering series prediction algorithm for 1.5D. The new algorithm changes the computing domain from pseudo-depth wavenumber domain to TX domain for predicting multiples. The improved algorithm demonstrated that the approach has some merits such as higher computing efficiency, feasibility to many kinds of geometries, lower predictive noise and independence to wavelet； 3. proposed a new subtraction algorithm. The new subtraction algorithm is not used to overcome nonorthogonality, but utilize the nonorthogonality's distribution in TX domain to estimate the true wavelet with filtering method. The method has excellent effectiveness in model testing. Improved 1D and 1.5D inverse scattering series algorithms can predict internal multiples. After filtering and subtracting among seismic traces in a window time, internal multiples can be attenuated in some degree. The proposed 1D and 1.5D algorithms have demonstrated that they are effective to the numerical and field data. In addition, the new subtraction algorithm is effective to the complex theoretic models.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

基于结构分析地震资料去噪方法研究

With the development of oil and gas exploration, the exploration of the continental oil and gas turns into the exploration of the subtle oil and gas reservoirs from the structural oil and gas reservoirs in China. The reserves of the found subtle oil and gas reservoirs account for more than 60 percent of the in the discovered oil and gas reserves. Exploration of the subtle oil and gas reservoirs is becoming more and more important and can be taken as the main orientation for the increase of the oil and gas reserves. The characteristics of the continental sedimentary facies determine the complexities of the lithological exploration. Most of the continental rift basins in East China have entered exploration stages of medium and high maturity. Although the quality of the seismic data is relatively good, this areas have the characteristics of the thin sand thickness, small faults, small range of the stratum. It requests that the seismic data have high resolution. It is a important task how to improve the signal/noise ratio of the high frequency of seismic data. In West China, there are the complex landforms, the deep embedding the targets of the prospecting, the complex geological constructs, many ruptures, small range of the traps, the low rock properties, many high pressure stratums and difficulties of boring well. Those represent low signal/noise ratio and complex kinds of noise in the seismic records. This needs to develop the method and technique of the noise attenuation in the data acquisition and processing. So that, oil and gas explorations need the high resolution technique of the geophysics in order to solve the implementation of the oil resources strategy for keep oil production and reserves stable in Ease China and developing the crude production and reserves in West China. High signal/noise ratio of seismic data is the basis. It is impossible to realize for the high resolution and high fidelity without the high signal/noise ratio. We play emphasis on many researches based on the structure analysis for improving signal/noise ratio of the complex areas. Several methods are put forward for noise attenuation to truly reflect the geological features. Those can reflect the geological structures, keep the edges of geological construction and improve the identifications of the oil and gas traps. The ideas of emphasize the foundation, give prominence to innovate, and pay attention to application runs through the paper. The dip-scanning method as the center of the scanned point inevitably blurs the edges of geological features, such as fault and fractures. We develop the new dip scanning method in the shap of end with two sides scanning to solve this problem. We bring forward the methods of signal estimation with the coherence, seismic wave characteristc with coherence, the most homogeneous dip-sanning for the noise attenuation using the new dip-scanning method. They can keep the geological characters, suppress the random noise and improve the s/n ratio and resolution. The rutine dip-scanning is in the time-space domain. Anew method of dip-scanning in the frequency-wavenumber domain for the noise attenuation is put forward. It use the quality of distinguishing between different dip events of the reflection in f-k domain. It can reduce the noise and gain the dip information. We describe a methodology for studying and developing filtering methods based on differential equations. It transforms the filtering equations in the frequency domain or the f-k domain into time or time-space domains, and uses a finite-difference algorithm to solve these equations. This method does not require that seismic data be stationary, so their parameters can vary at every temporal and spatial point. That enhances the adaptability of the filter. It is computationally efficient. We put forward a method of matching pursuits for the noise suppression. This method decomposes any signal into a linear expansion of waveforms that are selected from a redundant dictionary of functions. These waveforms are chosen in order to best match the signal structures. It can extract the effective signal from the noisy signal and reduce the noise. We introduce the beamforming filtering method for the noise elimination. Real seismic data processing shows that it is effective in attenuating multiples and internal multiples. The s/n ratio and resolution are improved. The effective signals have the high fidelity. Through calculating in the theoretic model and applying it to the real seismic data processing, it is proved that the methods in this paper can effectively suppress the random noise, eliminate the cohence noise, and improve the resolution of the seismic data. Their practicability is very better. And the effect is very obvious.

Consequências do trabalho por turnos e nocturno em profissionais de enfermagem (estudo realizado em dois hospitais do norte de Portugal)

Dissertação de Mestrado apresentada à Universidade Fernando Pessoa como parte dos requisitos para obtenção do grau de Mestre em Psicologia.

A construção ideográfica nas relações internacionais portuguesas 1945-1960

Tese apresentada à Universidade Fernando Pessoa como parte dos requisitos para obtenção do grau de Doutor em Ciências Sociais, especialidade em Relações Internacionais

Determinación de la producción de forraje y la calidad nutricional de diferentes secuencias de cultivos anuales invernales

El contexto actual que enfrentan los productores lecheros, caracterizado por un incremento de los costos de producción y una fuerte presión de la agricultura sobre la actividad tambo, hace que deban incrementar la productividad para poder mantener la rentabilidad de sus explotaciones. Aumentos en la productividad pueden lograrse mediante el incremento de la cantidad de forraje producido y utilizado por hectárea, por ejemplo implementando rotaciones forrajeras complementarias (CFR). Los objetivos del siguiente trabajo fueron comparar la oferta forrajera y la calidad nutricional de los cultivos anuales invernales que componen las CFR, frente a otros de mayor tradicion en la pampa húmeda. Se evaluó la cantidad de materia seca producida por hectárea (tn MS ha-1) de tres tratamientos: colza forrajera (Brassica rapa) - trebol persa (Trifolium repesinatum), avena (Avena sativa) y raigras anual (Lolium multiflorum). Se tomaron muestras de cada cultivo para la determinacion de la calidad nutricional: proteína cruda (PC), fibra detergente neutro (FDN), fibra detergente acido (FDA) y se estimó la digestibilidad de la materia seca (DMS). El diseño experimental utilizado fue de bloques completos aleatorizados con tres repeticiones. Los cultivos se sembraron a fines de marzo, se realizaron tres cortes con destino a pastoreo y un cuarto con destino a silaje para la avena y el raigras anual. En el caso de la colza forrajera, solo fue factible realizar dos cosechas debido a un pobre segundo rebrote posiblemente atribuído a la falta de precipitaciones durante ese período. Asimismo, el trébol persa intersembrado luego de la segunda cosecha de la colza, no logró germinar debido a las mismas condiciones expuestas anteriormente. Los tratamientos se analizaron estadísticamente mediante un análisis de varianza (ANOVA) y cuando correspondió, se realizó una prueba de comparaciones multiples (Tukey p . 0,05) para determinar las diferencias entre las medias. Los resultados obtenidos en el presente trabajo demostraron que la oferta forrajera total del tratamiento colza forrajera fue significativamente menor a los tratamientos avena y raigras anual (p . 0,05) debido principalmente al deficit hídrico que generó la pérdida de la tercer cosecha de colza forrajera y la falla en la siembra del trebol persa. En cuanto a la calidad forrajera, se observaron menores valores de FDN, FDA y consecuentemente mayor DMS para el tratamiento colza forrajera, sin embargo para PC no se hallaron resultados consistentes como para los casos anteriores. Los resultados han sido promisorios considerando que la oferta forrajera y calidad nutricional de los dos primeros cortes del tratamiento colza forrajera igualaron a la avena y raigras. Serán necesarias nuevas experiencias en secano y también bajo riego con el fin de generar mas informacion sobre estas especies.

Phenotypical characterization in "Plasmopara halstedii" (sunflower downy mildew) isolates of several races

Phenotypic variation (morphological and pathogenic characters), and genetic variability were studied in 50 isolates of seven Plasmopara halstedii (sunflower downy mildew) races 100, 300, 304, 314, 710, 704 and 714. There were significant morphological, aggressiveness, and genetic differences for pathogen isolates. However, there was no relationship between morphology of zoosporangia and sporangiophores and pathogenic and genetic characteristics for the races used in our study. Also, our results provided evidence that no relation between pathogenic traits and multilocus haplotypes may be established in P. halstedii. The hypothesis explaining the absence of relationships among phenotypic and genetic characteristics is discussed.

Searching for the IRA 'disappeared': Ground-penetrating radar investigation of a churchyard burial site, Northern Ireland

A search for the body of a victim of terrorist abduction and murder was made in a graveyard on the periphery of a major conurbation in Northern Ireland. The area is politically sensitive and the case of high profile. This required non-invasive, completely non-destructive and rapid assessment of the scene. A MALA RAMAC ground-penetrating radar system was used to achieve these objectives. Unprocessed and processed 400MHz data shows the presence of a collapse feature above and around a known 1970s burial with no similar collapse above the suspect location. In the saturated, clay-rich sediments of the site, 200MHz data offered no advantage over 400MHz data. Unprocessed 100MHz data shows a series of multiples in the known burial with no similar features in the suspect location. Processed 100MHz lines defined the shape of the collapse around the known burial to 2m depth, together with the geometry of the platform (1m depth) the gravedigger used in the 1970s to construct the site. In addition, processed 100MHz data showed both the dielectric contrast in and internal reflection geometry of the soil imported above the known grave. Thus the sequence, geometry, difference in infill and infill direction of the grave was reconstructed 30 years after burial. The suspect site showed no evidence of shallow or deep inhumation. Subsequently, the missing person������¢���¯���¿���½���¯���¿���½s body was found some distance from this site, vindicating the results and interpretation from ground-penetrating radar. The acquisition, processing, collapse feature and sequence stratigraphic interpretation of the known burial and empty (suspect) burial site may be useful proxies for other, similar investigations. GPR was used to evaluate this site within 3 hours of the survey commencing, using unprocessed data. An additional day of processing established that the suspect body did not reside here, which was counter to police and community intelligence.

Los nuevos territorios de la gráfica: imagen, proceso y distribución

La gráfica se encuentra en un proceso de simbiosis con la cultura visual en el que la tecnología de la imagen, la descentralización de la matriz y la adaptación del formato expositivo se funden con el resto de las actividades artísticas. El análisis de las comunidades formadas alrededor de esta práctica configura un puzzle, aparentemente bien encajado, que se divide entre el formalismo y el conceptualismo y entre la idea y el proceso. Las bienales y trienales, a la vanguardia del discurso, cuestionan cualquier concepto que se haya podido asentar con solidez: ante la decadencia del proceso artesanal el mensaje se profundiza nutriéndose de la invectiva provocada por las contradicciones sociales, territoriales y medioambientales para convertirse en un reflejo de la sociedad contemporánea. Comprender la gráfica actual pasa por un ejercicio de aperturismo y un proceso de adaptación al cambio técnico e iconográfico de los sectores implicados. El grabado es ahora un territorio en el que la integración de todas las artes permite cruzar sus fronteras con fluidez.

Remarks on common hypercyclic vectors

We treat the question of existence of common hypercyclic vectors for families of continuous linear operators. It is shown that for any continuous linear operator T on a complex Fréchet space X and a set ? ? R+ × C which is not of zero three-dimensional Lebesgue measure, the family {a T + b I : (a, b) ? ?} has no common hypercyclic vectors. This allows to answer negatively questions raised by Godefroy and Shapiro and by Aron. We also prove a sufficient condition for a family of scalar multiples of a given operator on a complex Fréchet space to have a common hypercyclic vector. It allows to show that if D = {z ? C : | z | < 1} and f ? H8 (D) is non-constant, then the family {z Mf{star operator} : b- 1 < | z | < a- 1} has a common hypercyclic vector, where Mf : H2 (D) ? H2 (D), Mf f = f f, a = inf {| f (z) | : z ? D} and b = sup {| f (z) | : | z | ? D}, providing an affirmative answer to a question by Bayart and Grivaux. Finally, extending a result of Costakis and Sambarino, we prove that the family {a Tb : a, b ? C {set minus} {0}} has a common hypercyclic vector, where Tb f (z) = f (z - b) acts on the Fréchet space H (C) of entire functions on one complex variable.

Les deux André et le concept de la durée: André Chamson et André Malraux dans le même combat pour les valeurs spirituelles.

Les vies et les carrières d’André Chamson et d’André Malraux se ressemblent par de multiples aspects. Contemporains, ils se connaissaient bien, menaient souvent les mêmes combats, souffraient des mêmes préoccupations, Mais c’est surtout par rapport à l’esprit qui les animait, en tant que grandes personnalités marquantes du XXe siècle, que d’intéressants parallèles peuvent être établis entre eux, que nous pouvons avec profit examiner les valeurs universelles et spirituelles qu’ils promouvaient, valeurs emblématiques pour tant de penseurs, d’acteurs, et d’écrivains de ce siècle. Ainsi pourrions-nous par ailleurs interroger et même évaluer la solidité, la résistance et le potentiel de telles valeurs pour le siècle à venir.