Seismic damage estimation through measurable dynamic characteristics

Ductility based design of reinforced concrete structures implicitly assumes certain damage under the action of a design basis earthquake. The damage undergone by a structure needs to be quantified, so as to assess the post-seismic reparability and functionality of the structure. The paper presents an analytical method of quantification and location of seismic damage, through system identification methods. It may be noted that soft ground storied buildings are the major casualties in any earthquake and hence the example structure is a soft or weak first storied one, whose seismic response and temporal variation of damage are computed using a non-linear dynamic analysis program (IDARC) and compared with a normal structure. Time period based damage identification model is used and suitably calibrated with classic damage models. Regenerated stiffness of the three degrees of freedom model (for the three storied frame) is used to locate the damage, both on-line as well as after the seismic event. Multi resolution analysis using wavelets is also used for localized damage identification for soft storey columns.

Seismic hazard analysis of Lucknow considering local and active seismic gaps

The Himalayas are one of very active seismic regions in the world where devastating earthquakes of 1803 Bihar-Nepal, 1897 Shillong, 1905 Kangra, 1934 Bihar-Nepal, 1950 Assam and 2011 Sikkim were reported. Several researchers highlighted central seismic gap based on the stress accumulation in central part of Himalaya and the non-occurrence of earthquake between 1905 Kangra and 1934 Bihar-Nepal. The region has potential of producing great seismic event in the near future. As a result of this seismic gap, all regions which fall adjacent to the active Himalayan region are under high possible seismic hazard due to future earthquakes in the Himalayan region. In this study, the study area of the Lucknow urban centre which lies within 350 km from the central seismic gap has been considered for detailed assessment of seismic hazard. The city of Lucknow also lies close to Lucknow-Faizabad fault having a seismic gap of 350 years. Considering the possible seismic gap in the Himalayan region and also the seismic gap in Lucknow-Faizabad fault, the seismic hazard of Lucknow has been studied based on deterministic and the probabilistic seismic hazard analysis. Results obtained show that the northern and western parts of Lucknow are found to have a peak ground acceleration of 0.11-0.13 g, which is 1.6- to 2.0-fold higher than the seismic hazard compared to the other parts of Lucknow.

Physical modelling of flexible retaining walls under seismic actions

A series of dynamic centrifuge tests on reduced scale models of flexible retaining structures were conducted on the Turner beam centrifuge at the Schofield Centre of the University of Cambridge. The paper illustrates the main results of the experimental work in terms of observed amplifications of ground motion and mobilised shear stiffness and damping ratio for all tests. The experimental results for one test on a pair of cantilevered walls in dense sand are also presented in terms of measured bending moments and horizontal displacements of the walls during (maximum values) and at the end of (residual values) each seismic event. Finally, the experimental data are discussed in the light of the results obtained from dynamic numerical analyses of the behaviour of cantilevered walls under real seismic actions. © 2010 Taylor & Francis Group, London.

Factors affecting seismic response of submarine slopes

The response of submerged slopes on the continental shelf to seismic or storm loading has become an important element in the risk assessment for offshore structures and "local" tsunami hazards worldwide. The geological profile of these slopes typically includes normally consolidated to lightly overconsolidated soft cohesive soils with layer thickness ranging from a few meters to hundreds of meters. The factor of safety obtained from pseudo-static analyses is not always a useful measure for evaluating the slope response, since values less than one do not necessarily imply slope failure with large movements of the soil mass. This paper addresses the relative importance of different factors affecting the response of submerged slopes during seismic loading. The analyses use a dynamic finite element code which includes a constitutive law describing the anisotropic stress-strain-strength behavior of normally consolidated to lightly overconsolidated clays. The model also incorporates anisotropic hardening to describe the effect of different shear strain and stress histories as well as bounding surface principles to provide realistic descriptions of the accumulation of the plastic strains and excess pore pressure during successive loading cycles. The paper presents results from parametric site response analyses on slope geometry and layering, soil material parameters, and input ground motion characteristics. The predicted maximum shear strains, permanent deformations, displacement time histories and maximum excess pore pressure development provide insight of slope performance during a seismic event. © 2006 Author(s). This work is licensed under a Creative Commons License.

Seismic response of submarine slopes

The geological profile of submerged slopes on the continental shelf typically includes soft cohesive soils with thicknesses ranging from a few meters to tens or hundreds of meters. The response of these soils in simple shear tests is largely influenced by the presence of an initial consolidation shear stress, inducing anisotropic stress-strain-strength properties which depend also on the direction of shear. In this paper, a new simplified effective-stress-based model describing the behavior of normally to lightly overconsolidated cohesive soils is used in conjunction with a one-dimensional seismic site response analysis computer code to illustrate the importance of accounting for anisotropy and small strain nonlinearity. In particular, a simple example is carried out to compare results for different slope inclinations. Depth profiling of the maximum shear strains and permanent deformations provide insight into the mechanisms of deformation during a seismic event, and the effects of sloping ground conditions.

Seismic vulnerability and risk assessment of historic masonry buildings

Seismic risk evaluation of built-up areas involves analysis of the level of earthquake hazard of the region, building vulnerability and exposure. Within this approach that defines seismic risk, building vulnerability assessment assumes great importance, not only because of the obvious physical consequences in the eventual occurrence of a seismic event, but also because it is the one of the few potential aspects in which engineering research can intervene. In fact, rigorous vulnerability assessment of existing buildings and the implementation of appropriate retrofitting solutions can help to reduce the levels of physical damage, loss of life and the economic impact of future seismic events. Vulnerability studies of urban centresshould be developed with the aim of identifying building fragilities and reducing seismic risk. As part of the rehabilitation of the historic city centre of Coimbra, a complete identification and inspection survey of old masonry buildings has been carried out. The main purpose of this research is to discuss vulnerability assessment methodologies, particularly those of the first level, through the proposal and development of a method previously used to determine the level of vulnerability, in the assessment of physical damage and its relationship with seismic intensity.

A comprehensive study of the 26th December 2004 Sumatra earthquake: possible implications on Earth rotation and investigations on the coseismic and postseismic stress diffusion associated with the seismic rupture

In this work a multidisciplinary study of the December 26th, 2004 Sumatra earthquake has been carried out. We have investigated both the effect of the earthquake on the Earth rotation and the stress field variations associated with the seismic event. In the first part of the work we have quantified the effects of a water mass redistribution associated with the propagation of a tsunami wave on the Earth’s pole path and on the length-of-day (LOD) and applied our modeling results to the tsunami following the 2004 giant Sumatra earthquake. We compared the result of our simulations on the instantaneous rotational axis variations with some preliminary instrumental evidences on the pole path perturbation (which has not been confirmed yet) registered just after the occurrence of the earthquake, which showed a step-like discontinuity that cannot be attributed to the effect of a seismic dislocation. Our results show that the perturbation induced by the tsunami on the instantaneous rotational pole is characterized by a step-like discontinuity, which is compatible with the observations but its magnitude turns out to be almost one hundred times smaller than the detected one. The LOD variation induced by the water mass redistribution turns out to be not significant because the total effect is smaller than current measurements uncertainties. In the second part of this work of thesis we modeled the coseismic and postseismic stress evolution following the Sumatra earthquake. By means of a semi-analytical, viscoelastic, spherical model of global postseismic deformation and a numerical finite-element approach, we performed an analysis of the stress diffusion following the earthquake in the near and far field of the mainshock source. We evaluated the stress changes due to the Sumatra earthquake by projecting the Coulomb stress over the sequence of aftershocks taken from various catalogues in a time window spanning about two years and finally analyzed the spatio-temporal pattern. The analysis performed with the semi-analytical and the finite-element modeling gives a complex picture of the stress diffusion, in the area under study, after the Sumatra earthquake. We believe that the results obtained with the analytical method suffer heavily for the restrictions imposed, on the hypocentral depths of the aftershocks, in order to obtain the convergence of the harmonic series of the stress components. On the contrary we imposed no constraints on the numerical method so we expect that the results obtained give a more realistic description of the stress variations pattern.

Multiresolution spherical wavelet analysis in global seismic tomography

Every seismic event produces seismic waves which travel throughout the Earth. Seismology is the science of interpreting measurements to derive information about the structure of the Earth. Seismic tomography is the most powerful tool for determination of 3D structure of deep Earth's interiors. Tomographic models obtained at the global and regional scales are an underlying tool for determination of geodynamical state of the Earth, showing evident correlation with other geophysical and geological characteristics. The global tomographic images of the Earth can be written as a linear combinations of basis functions from a specifically chosen set, defining the model parameterization. A number of different parameterizations are commonly seen in literature: seismic velocities in the Earth have been expressed, for example, as combinations of spherical harmonics or by means of the simpler characteristic functions of discrete cells. With this work we are interested to focus our attention on this aspect, evaluating a new type of parameterization, performed by means of wavelet functions. It is known from the classical Fourier theory that a signal can be expressed as the sum of a, possibly infinite, series of sines and cosines. This sum is often referred as a Fourier expansion. The big disadvantage of a Fourier expansion is that it has only frequency resolution and no time resolution. The Wavelet Analysis (or Wavelet Transform) is probably the most recent solution to overcome the shortcomings of Fourier analysis. The fundamental idea behind this innovative analysis is to study signal according to scale. Wavelets, in fact, are mathematical functions that cut up data into different frequency components, and then study each component with resolution matched to its scale, so they are especially useful in the analysis of non stationary process that contains multi-scale features, discontinuities and sharp strike. Wavelets are essentially used in two ways when they are applied in geophysical process or signals studies: 1) as a basis for representation or characterization of process; 2) as an integration kernel for analysis to extract information about the process. These two types of applications of wavelets in geophysical field, are object of study of this work. At the beginning we use the wavelets as basis to represent and resolve the Tomographic Inverse Problem. After a briefly introduction to seismic tomography theory, we assess the power of wavelet analysis in the representation of two different type of synthetic models; then we apply it to real data, obtaining surface wave phase velocity maps and evaluating its abilities by means of comparison with an other type of parametrization (i.e., block parametrization). For the second type of wavelet application we analyze the ability of Continuous Wavelet Transform in the spectral analysis, starting again with some synthetic tests to evaluate its sensibility and capability and then apply the same analysis to real data to obtain Local Correlation Maps between different model at same depth or between different profiles of the same model.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Medium change monitoring using ambient seismic noise and coda wave interferometry: examples from intraplate NE Brazil and the Mid-Atlantic Ridge

This thesis presents and discusses the results of ambient seismic noise correlation for two different environments: intraplate and Mid-Atlantic Ridge. The coda wave interferometry method has also been tested for the intraplate data. Ambient noise correlation is a method that allows to retrieve the structural response between two receivers from ambient noise records, as if one of the station was a virtual source. It has been largely used in seismology to image the subsurface and to monitor structural changes associated mostly with volcanic eruptions and large earthquakes. In the intraplate study, we were able to detect localized structural changes related to a small earthquake swarm, which main event is mR 3.7, North-East of Brazil. We also showed that the 1-bit normalization and spectral whitening result on the loss of waveform details and that the phase auto-correlation, which is amplitude unbiased, seems to be more sensitive and robust for our analysis of a small earthquake swarm. The analysis of 6 months of data using cross-correlations detect clear medium changes soon after the main event while the auto-correlations detect changes essentially after 1 month. It could be explained by fluid pressure redistribution which can be initiated by hydromechanical changes and opened path ways to shallower depth levels due to later occurring earthquakes. In the Mid-Atlantic Ridge study, we investigate structural changes associated with a mb 4.9 earthquake in the region of the Saint Paul transform fault. The data have been recorded by a single broadband seismic station located at less than 200 km from the Mid-Atlantic ridge. The results of the phase auto-correlation for a 5-month period, show a strong co-seismic medium change followed by a relatively fast post-seismic recovery. This medium change is likely related to the damages caused by the earthquake’s ground shaking. The healing process (filling of the new cracks) that lasted 60 days can be decomposed in two phases, a fast recovery (70% in ~30 days) in the early post-seismic stage and a relatively slow recovery later (30% in ~30 days). In the coda wave interferometry study, we monitor temporal changes of the subsurface caused by the small intraplate earthquake swarm mentioned previously. The method was first validated with synthetics data. We were able to detect a change of 2.5% in the source position and a 15% decrease of the scatterers’ amount. Then, from the real data, we observed a rapid decorrelation of the seismic coda after the mR 3.7 seismic event. This indicates a rapid change of the subsurface in the fault’s region induced by the earthquake.

Medium change monitoring using ambient seismic noise and coda wave interferometry: examples from intraplate NE Brazil and the Mid-Atlantic Ridge

This thesis presents and discusses the results of ambient seismic noise correlation for two different environments: intraplate and Mid-Atlantic Ridge. The coda wave interferometry method has also been tested for the intraplate data. Ambient noise correlation is a method that allows to retrieve the structural response between two receivers from ambient noise records, as if one of the station was a virtual source. It has been largely used in seismology to image the subsurface and to monitor structural changes associated mostly with volcanic eruptions and large earthquakes. In the intraplate study, we were able to detect localized structural changes related to a small earthquake swarm, which main event is mR 3.7, North-East of Brazil. We also showed that the 1-bit normalization and spectral whitening result on the loss of waveform details and that the phase auto-correlation, which is amplitude unbiased, seems to be more sensitive and robust for our analysis of a small earthquake swarm. The analysis of 6 months of data using cross-correlations detect clear medium changes soon after the main event while the auto-correlations detect changes essentially after 1 month. It could be explained by fluid pressure redistribution which can be initiated by hydromechanical changes and opened path ways to shallower depth levels due to later occurring earthquakes. In the Mid-Atlantic Ridge study, we investigate structural changes associated with a mb 4.9 earthquake in the region of the Saint Paul transform fault. The data have been recorded by a single broadband seismic station located at less than 200 km from the Mid-Atlantic ridge. The results of the phase auto-correlation for a 5-month period, show a strong co-seismic medium change followed by a relatively fast post-seismic recovery. This medium change is likely related to the damages caused by the earthquake’s ground shaking. The healing process (filling of the new cracks) that lasted 60 days can be decomposed in two phases, a fast recovery (70% in ~30 days) in the early post-seismic stage and a relatively slow recovery later (30% in ~30 days). In the coda wave interferometry study, we monitor temporal changes of the subsurface caused by the small intraplate earthquake swarm mentioned previously. The method was first validated with synthetics data. We were able to detect a change of 2.5% in the source position and a 15% decrease of the scatterers’ amount. Then, from the real data, we observed a rapid decorrelation of the seismic coda after the mR 3.7 seismic event. This indicates a rapid change of the subsurface in the fault’s region induced by the earthquake.

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3-8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal-India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs. Crown Copyright (C) 2013 Published by Elsevier Ltd. All rights reserved.

Damage detection in civil structures using high-frequency seismograms

The dynamic properties of a structure are a function of its physical properties, and changes in the physical properties of the structure, including the introduction of structural damage, can cause changes in its dynamic behavior. Structural health monitoring (SHM) and damage detection methods provide a means to assess the structural integrity and safety of a civil structure using measurements of its dynamic properties. In particular, these techniques enable a quick damage assessment following a seismic event. In this thesis, the application of high-frequency seismograms to damage detection in civil structures is investigated.

Two novel methods for SHM are developed and validated using small-scale experimental testing, existing structures in situ, and numerical testing. The first method is developed for pre-Northridge steel-moment-resisting frame buildings that are susceptible to weld fracture at beam-column connections. The method is based on using the response of a structure to a nondestructive force (i.e., a hammer blow) to approximate the response of the structure to a damage event (i.e., weld fracture). The method is applied to a small-scale experimental frame, where the impulse response functions of the frame are generated during an impact hammer test. The method is also applied to a numerical model of a steel frame, in which weld fracture is modeled as the tensile opening of a Mode I crack. Impulse response functions are experimentally obtained for a steel moment-resisting frame building in situ. Results indicate that while acceleration and velocity records generated by a damage event are best approximated by the acceleration and velocity records generated by a colocated hammer blow, the method may not be robust to noise. The method seems to be better suited for damage localization, where information such as arrival times and peak accelerations can also provide indication of the damage location. This is of significance for sparsely-instrumented civil structures.

The second SHM method is designed to extract features from high-frequency acceleration records that may indicate the presence of damage. As short-duration high-frequency signals (i.e., pulses) can be indicative of damage, this method relies on the identification and classification of pulses in the acceleration records. It is recommended that, in practice, the method be combined with a vibration-based method that can be used to estimate the loss of stiffness. Briefly, pulses observed in the acceleration time series when the structure is known to be in an undamaged state are compared with pulses observed when the structure is in a potentially damaged state. By comparing the pulse signatures from these two situations, changes in the high-frequency dynamic behavior of the structure can be identified, and damage signals can be extracted and subjected to further analysis. The method is successfully applied to a small-scale experimental shear beam that is dynamically excited at its base using a shake table and damaged by loosening a screw to create a moving part. Although the damage is aperiodic and nonlinear in nature, the damage signals are accurately identified, and the location of damage is determined using the amplitudes and arrival times of the damage signal. The method is also successfully applied to detect the occurrence of damage in a test bed data set provided by the Los Alamos National Laboratory, in which nonlinear damage is introduced into a small-scale steel frame by installing a bumper mechanism that inhibits the amount of motion between two floors. The method is successfully applied and is robust despite a low sampling rate, though false negatives (undetected damage signals) begin to occur at high levels of damage when the frequency of damage events increases. The method is also applied to acceleration data recorded on a damaged cable-stayed bridge in China, provided by the Center of Structural Monitoring and Control at the Harbin Institute of Technology. Acceleration records recorded after the date of damage show a clear increase in high-frequency short-duration pulses compared to those previously recorded. One undamage pulse and two damage pulses are identified from the data. The occurrence of the detected damage pulses is consistent with a progression of damage and matches the known chronology of damage.

高邮凹陷陈堡及陈堡东地区油气成藏规律研究

With the development of petroleum exploration in Gaoyou Depression, both old and new areas have been the active exploration targets, so the study of petroleum accumulation is significant to the petroleum exploration in the study area and the integrated oil and gas accumulation theory. Based on hydrocarbon accumulation theory and systematical research methods and combined with the structural characteristics of Gaoyou Depression, Chenbao and East of Chenbao were selected as the study areas in this dissertation, oil and gas migration pathways, accumulation periods, as well as accumulation models were studied, and favorable exploration targets were proposed. There develop three sets source rocks, which are Tai-2 Member, Fu-2 Member and Fu-4 Member respectively. Tai-2 Member is the predominant source rock in the eastern part. Fu-2 Member mainly occurs in the northern slope, while Fu-4 Member develops in the deep depression. In the study area, oil mainly comes from Fu-2 Member of Liuwushe subsag. The lower limit of TOC is 0.4%, and active source rock mostly distributed in the south fault-step zone. The source rock in Liuwushe subsag began to generate hydrocarbon in the late of Dainan depositional stage and the threshold was 2300m. The macro and micro characteristics of reservoirs and the reservoir heterogeneity characteristics of the Fu-1 Member were studied systematicly. The results show that Fu-1 Member, which has better reservoir properties, are medium porosity-medium permeability reservoir. The reservoir permeability has good correlation with porosity connectivity. The reservoirs have strong dissolution, pores are mainly thin to medium throat, and throat radii are distributed concentratedly, the sorting is good and pore structures are homogeneous. Sandstone reservoirs whether in the plan view, interlayer or in layers have a certain degree of heterogeneity, in particular, the heterogeneity in layers directly affect and control the oil and gas migration and accumulation. By analyzing the lithology correlation of the fault walls, shale smear, cross section stress, the configuration of fracture active periods and hydrocarbon generation and expulsion periods and fuzzy comprehensive evaluation, the main faults sealing were evaluated. The results show that the faults in Chenbao and East of Chenbao had poor sealing properties in Sanduo period and could be used as the migration pathways at that time. After Sanduo period, the tectonic stress fields in the area changed largely, and, consequently, the fault properties converted from tensional shear to compressive shear, the faults changed progressively from close to open, so the faults sealing became better and were conducive to the preservation of oil and gas reservoirs. According to the seismic event suspension modes and profile configurations above and under the unconformities, combined with tectonic evolutions of the study areas, the unconformity types can be classified into truncation unconformity, overlapped unconformity and parallel unconformity and the distribution characteristics of unconformities in the plan view was also studied. The unconformity structure was divided into basal conglomerate, weathered clay and semi-weathered layer vertically in the study area and this kind of structure make unconformities to be effective oil and gas migration pathways and is significant to hydrocarbon accumulation in a parts of areas. With the analyses of typical oil and gas reservoirs in the study area, combined with the research results of pathway systems, hydrocarbon accumulation models were established and the oil and gas accumulation laws in Chenbao and East of Chenbao analyzed. The oil and gas came from Liuwushe subsag and Liuliushe subsag. The oil and gas from Liuwushe subsag mainly migrated from the structural high parts into the fault-step zone along strata in northeast direction, a part of them migrated upward into the fault-step zone and the Wubao low uplift along Wu-1 Fault in northeast direction. The oil and gas from Liuliushe subsag mainly migrated into the upper reservoirs through Wu-2 fault, and lesser oil and gas migrated into the fault-step zone because of the controls of cross-section orientation, depression center and the hydrocarbon formation tendency. The favorable exploration targets in Chenbao and East of Chenbao have been concluded: the southern fault-step zone is a favorable oil and gas accumulation zone of Liuwushe subsag, and they are fault block reservoirs where fault acted as the barriers, the main target intervals are Fu-1 Member and Fu-3 Member in palaeocene; Oil and gas in the middle and northern fault-step zone mainly laterally migrated from the south areas, and the main target interval is Fu-3 Member in palaeocene; Fu-1 Member and the reserviors above the Wubao subsag are still the focuses in future explorations. The results of this study have important guiding significance for the future oil and gas exploration.

The Cyrenaican Prehistory Project 2012: the sixth season of excavations in the Haua Fteah cave

The paper reports on the sixth season of fieldwork of the Cyrenaican Prehistory Project (CPP) undertaken in September 2012. As in the spring 2012 season, work focussed on the Haua Fteah cave and on studies of materials excavated in previous seasons, with no fieldwork undertaken elsewhere in the Gebel Akhdar. An important discovery, in a sounding excavated below the base of McBurney's 1955 Deep Sounding (Trench S), is of a rockfall or roof collapse conceivably dating to the cold climatic regime of Marine Isotope Stage (MIS) 6 (globally dated to c. 190-130 ka) but more likely the result of a seismic event within MIS 5 (globally dated to c. 130-80 ka). The sediments and associated molluscan fauna in Trench S and in Trench D, a trench being cut down the side of the Deep Sounding, indicate that this part of the cave was at least seasonally waterlogged during the accumulation, probably during MIS 5, of the -6.5 rn of sediment cut through by the Deep Sounding. Evidence for human fréquentation of the cave in this period is more or less visible depending on how close the trench area was to standing water as it fluctuated through time. Trench M, the trench being cut down the side of McBurney's Middle Trench, has now reached the depth of the latest Middle Stone Age or Middle Palaeolithic (Levalloiso-Mousterian) industries. The preliminary indications from its excavation are that the transition from the Levalloiso-Mousterian to the blade-based Upper Palaeolithic or Late Stone Age Dabban industry was complex and perhaps protracted, at a time when the climate was oscillating between warmstage stable environmental conditions and colder and more arid environments. The estimated age of the sediments, c. 50-40 ka, places these oscillations within the earlier part of MIS 3 (globally dated to 60-24 ka), when global climates experienced rapid fluctuations as part of an overall trend to increasing aridity and cold.