天然气藏地震识别技术研究

The Research of Seismic Recognition Techniques for Gas Reservoir Shang Yong_sheng(Geophysics) Directed by Yang Chang-chun Abstract Gas reservior is one of the most important nature resources. Someone forecast that the output will exceed crude oil in 2015 and become the largest energy source. Recently,more and more gas reservior are discovered as the oil field and gas filed exploration go deep into development. Although the gas proved reserves rise greatly the explorative degree of natural gas resource in our country is still very low. The potential of gas exploration is very great and our task is so hard. How to recognise and discover new gas reservoir is the first task based on the great gas reservior resources foreground. the gas reservior in different oil and gas field have its special gas generation, reservoiring, physical property conditions. However,it may have the same geophysical characters. So,it is very important to analyse, research, summarizing the geophysical characters of the gas reservior and make use of the characters to identify the gas layer effectively. This paper start with modeling,and it discuss the geophysical characters of the gas reservior response. It analyse the seismical wave characters of the gas reservoir. Furthermore, it summarize the method of using the seismica profile to identify the gas reservior directly. The paper discuss the research of extracting diffraction wave for mass diffraction wave grow at the edge of the gas reservoir at the seismic section. Making use of the technique of extracting diffraction wave to identify the gas reservior is the first experiment of the gas reservoir prediction technique. The avo technology is a new geophysical method. From the pre-stack analysis, this paper discuss the technique of using the rich information to identify the gas reservoir. Based on the case study of the Qidam basin and the Hailaer basin it discuss the method of predicating gas reservoir using pre-stack information. It include pre-stack amplitude preserve process, AVO modeling, fluid replacement technique, AVO analysis and interpretation technique. The paper summarize a gas reservoir prediction procedure focusing on the pre-stack information. The seismic wave will cause great attenuation when it pass through the gas layer and the high frequency component loss. This paper discuss the technique of extracting seismic attributes to represent the attenuation degree of seismica wave. Based on the attenuation attributes,it does the research of the gas reservor identification and prediction. At last,the paper discuss the method of calculating the azimuthal anisotropy to predict the fracture reservoir. Keyword: gas reservoir, diffraction wave, AVO, attenuation attribute,fracture prediction

残留盆地分布综合地球物理研究——以大港地区为例

Macro-distribution of residual basins is a basic question in residual basin research，the main object of macro-distribution study is to build strata framework, compute thickness of residual strata and analyze characteristics of residual basins. With the guidance of the theory of integrated geology and geophysical research, the paper assembled series of methods and established the technical chart based on gravity and magnetic data, with restriction of geology, seismic and drilling data. Based on potential field data processing and analysis, forward and inverse computation, region potential field analysis and potential field separation, etc. it computed depth of gravity/magnetic basement and got strata framework. It had got effective results in the research of macro-distribution of residual basin research in the Dagang area. It did the wavelet transform of gravity/magnetic data with multi-kind of wavelet basis using a trou algorithm. From comparison of processing result and their spectral of wavelet analysis, up continuation and filter method, the wavelet approximation is better to fit the regional potential field, and it is an effective method to separate gravity/magnetic effect caused by deep geology bodies. The experiment of matching pursuit shows that te transform domain methods have great advantage in potential data analysis. From the integrated geophysical study of rock property study, gravity/magnetic basement inversion and fault system analysis of the Dagang area, it gets the strata framework and the thickness of pre-Cenozoic residual strata. Comprehensive study with gravity and magnetotelluric profile inversion and interpretation, three prospect plays of macro-distribution of residual basins are fingered out. It has great residual strata thickness in the northern part of Chengning Uplift and there is thrust fault in the deep zone and good up-Paleozoic hydrocarbon source rocks in this area. With integrated analysis, this area will be the most prospective hydrocarbon location of pre-Cenozoic residual basins.

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

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.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

非自喷层测试及评价解释技术研究与应用

Now low porosity and low permeability reservoir is one of the main targets of exploration for the onshore oilfields of China. Most of the reservoirs are none flowing because of bad formation percolation condition, poor gas oil ratio , low formation pressure coefficient and other factors. In the recent years, a number of domestic oilfields have carried out some research work and achieved some success on oil testing and production technology in such formation. But by now, there is still no systematic and mature technology, particularly testing technology in none flowing formation is still needed further study. Based on study the key problem of well testing and interpretation technology in none flowing formation, solve the important problems in well testing technology, continuously improve and innovate geological information acquisition technology for none flowing reservoir, accurately acquire boundary information and evaluate reservoir flow characteristics. Its wide application remarkable result has shown. The main results and cognitions obtained from research are as follows: 1. This new technology research results help solve the occurrent problems in well testing process for none flowing formations, such as small investigation radius, poor representative of interpretation results from the poor data, low level application of interpretation results. This new technology helps create favorable conditions for early precise reservoir evaluation and reduction of the risk of exploration. 2. The technological difficulties for none flowing well testing are successfully solved by using none flowing formation combined mechanical tool string .This method has been proved by its applications to be able to improve the efficiency of the testing and the quantity of the acquired test data ,and so as to enhance the application of the interpretation results of the test in development of oil fields. 3. The application of the rotary formation tester, selective test valve, well testing string and their allier tools help to resolve problems such as the operation of opening and shutting-in the well under different well conditions, to broaden the scope of well test technology for none flowing formations. 4. Refined Testing Technique for production Wells has greatly shortened the testing dwration and improved the efficiency and accuracy of operation, enriched test results, and at the same time created conditions for conducting multi-well interference well testing.

青藏高原东部长周期大地电磁测深研究

As a complement to conventional MT, Long-period Magnetotellurics (LMT) has been developed at low frequency for soundings of deep electric structures. Eastern Himalayan Syntaxis (EHS) and surrounding area is a key place for the study of dynamics of the uplift of Tibetan plateau. Experiments in the pioneer studies for EHS3D project showed that the study area shares an unusual low resistive crust and upper mantle. Conventional MT could not provide sufficient information about the deep structures of the study area that requested long period MT measurement to be complemented. This thesis presents the LMT studies in eastern Tibet along the EHS3D-3 Profile from Xiachayu to Yushu including data acquisition, processing, inversion and interpretation. The effective period of the measured LMT signals extend from 10s up to 30000s for the duration more than one week measurement. The resulting model shows that the LMT sounding coincides with the MT data in overlapped periods. Especially the induction arrows and tippers derived from LMT data provide more information about the base of the conductors beneath the plateau with higher resolution. Anomalous induction coefficients and 2-D model suggest extensive conductive bodies beneath Lhasa block and Qiangtang terrain which would be a possible evidence for partial melt and fluids at depth.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

基于可复制性的独立成分分析——应用于功能磁共振成像数据

This dissertation systematically depicted and improved the application of Independent Component Analysis (ICA) to Functional Magnetic Resonance Imaging (fMRI), following the logic of verification, improvement, extension, and application. The concept of “reproducibility” was the philosophy throughout its four concluded studies. In the “verification” study, ICA was applied to the resting-state fMRI data, verified the resultant components with reproducibility, and examined the consistency of the results from ICA and traditional “seed voxel” method. At the meantime, the limitation of ICA application on fMRI data analysis was presented. In the “improvement” study, an improved ICA algorithm based on reproducibility, RAICAR, was developed to aid some of the limitations of ICA application. RAICAR was able to rank ICA components by reproducibility, determine the number of reliable components, and obtain more stable results. RAICAR provided useful tools for validation and interpretation of ICA results. In the “extension” study, RAICAR as well as the concept of “reproducibility” was extended to multi-subject ICA analysis, and gRAICAR algorithm was developed. gRAICAR allows some variation across subjects, examining common components among subjects. gRAICAR is also capable to detect potential subject grouping on some components. It is a new way for exploratory group analysis on fMRI. In the “application” study, two newly developed methods, RAICAR and gRAICAR, were used to investigate the effect of early music training on the brain mechanism of memory and learning. The results showed brain mechanism difference in memory retrieval and learning process between two groups of subjects. This study also verified the usefulness and importance of the new methods.

An artificial immune system for data analysis

Timmis J and Neal M J. An artificial immune system for data analysis. In Proceedings of 3rd international workshop on information processing in cells and tissues (IPCAT), Indianapolis, U.S.A., 1999.

Accurate Range Image Registration: Eliminating or Modelling Outliers

Liu, Yonghuai, Liu, Honghai, Li, Longzhuang, Wei, Baogang. Accurate Range Image Registration: Eliminating or Modelling Outliers. Proceedings of 12th IEEE Conference on Emerging Technologies and Factory Automation, 2007, pp. 1316-1323. Sponsorship: IEEE

Społeczna interpretacja wybranych wczesnośredniowiecznych cmentarzysk rzędowych w perspektywie gender studies i lifecycle analysis

Wydział Historyczny: Instytut Prahistorii