基于正演方法的井间地震响应特征研究

地震勘探的基本任务是测量由人工震源激发的各种地震波，通过对所得到的地震波场的数据记录进行必要的处理和解释而提取所需要的地质信息。不同的地震波场特征反应了不同的地质构造等信息，因此，对地震波场特征进行研究就显得非常重要。井间地震是在一口井中激发、另一口井中接收的一种特殊的地震勘探方法，其实际记录要比常规的地震记录复杂的多，其波场中的各种波也更丰富复杂，如何识别不同类型的波并研究它们的响应特征是一个非常必要的任务。正演模拟可以为地震数据采集、处理、解释提供理论依据，是井间地震波场特征分析的一种有效方法。 本文主要研究了井间地震的数值正演模拟问题，系统地介绍了射线追踪正演模拟和波动方程正演模拟方法，利用Visual Basic 6.0语言编写了相应的射线追踪正演模拟程序，而且分两种方法详细分析了井间地震数值模型的地震响应特征，为直观的井间地震波场识别和勘探地震数据反演提供了必要的依据。本文取得的主要成果如下： 1)在前人的VSP正演模拟的基础上，改进并总结出了一套井间地震射线追踪正演模型的制作方法。 2)编制了井间地震射线追踪正演模拟程序V1.0。 3)利用编写的井间地震正演模拟程序，作了两种射线追踪数值模拟试算；利用Tesseral声波方程正演模拟试算了含有楔形体和丘状体的理论模型；取得了较好的模拟数据。

大陆地壳岩石结构重建：方法与实例

The composition of the continental crust has long been a subject of interest to earth scientists as it can provide key information about the crustal growth and evolution of the continents. In this paper we make a comparative study on the lithological discrimination schemes featuring with the use of different seismic attributes, such as P-wave velocity, P- to S-wave velocity ratio, acoustic or elastic impedances, Lame impedances and high-sensitive identification factors. The results demonstrate that Lame impedances have more powerful constrains than other seismic attributes. In order to fully take the advantage of make the best of the different seismic response of crustal rock, we firstly use seismic attribute that have weak distinguish power to construct loose constrained lithological model, then use seismic attributes that have stronger distinguish power to tighten the constrains of lithological discrimination. We propose a joint scheme (chain constrain technique) by combing all available constrains to reduce the non-uniqueness in mapping rock distribution. We adopt chain constrain technique to construct lithological model beneath Tunxi-Wenzhou transect, Southeastern China, Manzhouli-Suifenhe transect, Northeastern China, and geophysical profile in Bohai Bay Basin, North China. The results can be suumarized as the follows: (1) We compare the sensitivity of different seismic factor constraints on rock types, and conclude that Lame impedances have tighter constrains than seismic velocity, Vp/Vs, density. (2) We propose chain constrains to construct lithological model from integrated geophysical data, and reduce the non-uniqueness in mapping rock distribution. (3) We reconstruct crustal lithological model beneath Tunxi-Wenzhou transect, Southeastern China. The results suggested that Jiangshan-Shaoxing fault is a crust-scale, and it is the boundary between Cathaysia and Yanthze blocks. (4) We construct crustal lithological model beneath Manzhouli-Suifenhe transect, Northeastern China. (5) We map the petrologic distribution along a geophysical profile in Bohai Bay Basin, North China, and construct a three-layered petrology model from the depth 2 km to about 10 km, consisted of basalt (the first layer), pelitic siltstone (the second layer), and silty mudstone and fine sandstone (the third layer).

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

Transcriptome analysis of sulfate deficiency response in marine microalga Emiliania huxleyi

Summary The response to sulfate deficiency of plants and freshwater green algae has been extensively analysed by system biology approaches. By contrast, seawater sulfate concentration is high and very little is known about the sulfur metabolism of marine organisms. Here, we used a combination of metabolite analysis and transcriptomics to analyse the response of the marine microalga Emiliania huxleyi as it acclimated to sulfate limitation. Lowering sulfate availability in artificial seawater from 25 to 5 mM resulted in significant reduction in growth and intracellular concentrations of dimethylsulfoniopropionate and glutathione. Sulfate-limited E. huxleyi cells showed increased sulfate uptake but sulfate reduction to sulfite did not seem to be regulated. Sulfate limitation in E. huxleyi affected expression of 1718 genes. The vast majority of these genes were upregulated, including genes involved in carbohydrate and lipid metabolism, and genes involved in the general stress response. The acclimation response of E. huxleyi to sulfate deficiency shows several similarities to the well-described responses of Arabidopsis and Chlamydomonas, but also has many unique features. This dataset shows that even though E. huxleyi is adapted to constitutively high sulfate concentration, it retains the ability to re-program its gene expression in response to reduced sulfate availability.

Characterisation of the cyclic behaviour of reinforced concrete elements with plain bars

Sismos recentes comprovam a elevada vulnerabilidade dos edifícios existentes de betão armado. A resposta das estruturas aos sismos é fortemente condicionada pelas características da aderência aço-betão, que exibe degradação das propriedades iniciais quando sujeitas a carregamentos cíclicos e alternados. Este fenómeno é ainda mais gravoso para elementos com armadura lisa, predominantes na maioria das estruturas construídas até à década de 70 nos países do sul da Europa. A prática corrente de conceção, dimensionamento e pormenorização das estruturas antigas leva a que tenham características de comportamento e níveis de segurança associados não compatíveis com as exigências atuais. Os estudos realizados sobre o comportamento cíclico de elementos estruturais de betão armado com armadura lisa são ainda insuficientes para a completa caracterização deste tipo de elementos. Esta tese visou a caraterização da relação tensão de aderência versus escorregamento para elementos estruturais com armadura lisa e o estudo da resposta cíclica de pilares e nós viga-pilar de betão armado com armadura lisa. Foram realizados dez séries de ensaios de arrancamento (nove monotónicos e um cíclico) em provetes com varões lisos. Os resultados destes ensaios permitiram propor novas expressões empíricas para a estimativa dos parâmetros usados num modelo disponível na literatura para representação da relação tensão de aderência versus escorregamento. É ainda proposto um novo modelo monotónico para a relação tensão de aderência versus escorregamento que representa melhor a resposta após a resistência máxima de aderência. Uma campanha de ensaios unidirecionais em pilares e nós viga-pilar foi também realizada com o objetivo principal de caracterizar o comportamento cíclico deste tipo de elementos. No total foram realizados oito ensaios em pilares, sete ensaios em nós viga-pilar interiores e seis ensaios em nós viga-pilar exteriores representativos de estruturas antigas de betão armado com armadura lisa. Os resultados experimentais permitiram avaliar a influência do escorregamento e estudar o mecanismo de corte em nós e a evolução dos danos para elementos com armadura lisa. Com base nos resultados experimentais foi proposta uma adaptação na expressão do Eurocódigo 8-3 para o cálculo da capacidade última de rotação de elementos com armadura lisa. Foi também desenvolvido um estudo paramétrico, com diferentes estratégias de modelação não linear, para a simulação da resposta de pilares considerando o escorregamento da armadura lisa. Por último, foi proposto um novo modelo simplificado trilinear para o aço que contempla o efeito do escorregamento da armadura lisa.

Shaking-table tests of flat-bottom circular silos containing grain-like material

According to Eurocode 8, the seismic design of flat-bottom circular silos containing grain-like material is based on a rough estimate of the inertial force imposed on the structure by the ensiled content during an earthquake: 80% of the mass of the content multiplied by the peak ground acceleration. A recent analytical consideration of the horizontal shear force mobilised within the ensiled material during an earthquake proposed by some of the authors has resulted in a radically reduced estimate of this load suggesting that, in practice, the effective mass of the content is significantly less than that specified. This paper describes a series of laboratory tests that featured shaking table and a silo model, which were conducted in order to obtain some experimental data to verify the proposed theoretical formulations and to compare with the established code provisions. Several tests have been performed with different heights of ensiled material – about 0.5 mm diameter Ballotini glass – and different magnitudes of grain–wall friction. The results indicate that in all cases, the effective mass is indeed lower than the Eurocode specification, suggesting that the specification is overly conservative, and that the wall–grain friction coefficient strongly affects the overturning moment at the silo base. At peak ground accelerations up to around 0.35 g, the proposed analytical formulation provides an improved estimate of the inertial force imposed on such structures by their contents.

Novel hybrid columns made of ultra-high performance concrete and fiber reinforced polymers

The application of advanced materials in infrastructure has grown rapidly in recent years mainly because of their potential to ease the construction, extend the service life, and improve the performance of structures. Ultra-high performance concrete (UHPC) is one such material considered as a novel alternative to conventional concrete. The material microstructure in UHPC is optimized to significantly improve its material properties including compressive and tensile strength, modulus of elasticity, durability, and damage tolerance. Fiber-reinforced polymer (FRP) composite is another novel construction material with excellent properties such as high strength-to-weight and stiffness-to-weight ratios and good corrosion resistance. Considering the exceptional properties of UHPC and FRP, many advantages can result from the combined application of these two advanced materials, which is the subject of this research. The confinement behavior of UHPC was studied for the first time in this research. The stress-strain behavior of a series of UHPC-filled fiber-reinforced polymer (FRP) tubes with different fiber types and thicknesses were tested under uniaxial compression. The FRP confinement was shown to significantly enhance both the ultimate strength and strain of UHPC. It was also shown that existing confinement models are incapable of predicting the behavior of FRP-confined UHPC. Therefore, new stress-strain models for FRP-confined UHPC were developed through an analytical study. In the other part of this research, a novel steel-free UHPC-filled FRP tube (UHPCFFT) column system was developed and its cyclic behavior was studied. The proposed steel-free UHPCFFT column showed much higher strength and stiffness, with a reasonable ductility, as compared to its conventional reinforced concrete (RC) counterpart. Using the results of the first phase of column tests, a second series of UHPCFFT columns were made and studied under pseudo-static loading to study the effect of column parameters on the cyclic behavior of UHPCFFT columns. Strong correlations were noted between the initial stiffness and the stiffness index, and between the moment capacity and the reinforcement index. Finally, a thorough analytical study was carried out to investigate the seismic response of the proposed steel-free UHPCFFT columns, which showed their superior earthquake resistance, as compared to their RC counterparts.

Response of tall buildings with symmetric setbacks under blast loading

This study explores three-dimensional nonlineardynamic responses of typical tall buildings with and without setbacks under blast loading. These 20 storey reinforced concrete buildings have been designed for normal (dead, live and wind)loads. The influence of the setbacks on the lateral load response due to blasts in terms of peak deflections, accelerations, inter-storey drift and bending moments at critical locations (including hinge formation) were investigated. Structural response predictions were performed with a commercially available three-dimensional finite element analysis programme using non-linear direct integration time history analyses. Results obtained for buildings with different setbacks were compared and conclusions made. The comparisons revealed that buildings have setbacks that protect the tower part above the setback level from blast loading show considerably better response in terms of peak displacement and interstorey drift, when compared to buildings without setbacks. Rotational accelerations were found to depend on the periods of the rotational modes. Abrupt changes in moments and shears are experienced near the levels of the setbacks. Typical twenty storey tall buildings with shear walls and frames that are designed for only normaln loads perform reasonably well, without catastrophic collapse, when subjected to a blast that is equivalent to 500 kg TNT at a standoff distance of 10 m.

Structural damage detection using frequency response functions and neural networks

This paper illustrates the damage identification and condition assessment of a three story bookshelf structure using a new frequency response functions (FRFs) based damage index and Artificial Neural Networks (ANNs). A major obstacle of using measured frequency response function data is a large size input variables to ANNs. This problem is overcome by applying a data reduction technique called principal component analysis (PCA). In the proposed procedure, ANNs with their powerful pattern recognition and classification ability were used to extract damage information such as damage locations and severities from measured FRFs. Therefore, simple neural network models are developed, trained by Back Propagation (BP), to associate the FRFs with the damage or undamaged locations and severity of the damage of the structure. Finally, the effectiveness of the proposed method is illustrated and validated by using the real data provided by the Los Alamos National Laboratory, USA. The illustrated results show that the PCA based artificial Neural Network method is suitable and effective for damage identification and condition assessment of building structures. In addition, it is clearly demonstrated that the accuracy of proposed damage detection method can also be improved by increasing number of baseline datasets and number of principal components of the baseline dataset.

Blast response and vulnerability assessment of piled foundations

This research treated the response of pile foundations to blast loads. The influence of important parameters was investigated. The research techniques and the results will enable safer design of pile foundations that are vulnerable to blast loads.

Effect of Oil Price on Emissions Trading Market: VAR Analysis

In this study, we investigated the relationship of European Union carbon dioxide CO2 allowances EUAs prices and oil prices by employing a VAR analysis, Granger causality test and impulse response function. If oil price continues increasing, companies will decrease dependency on fossil fuels because of an increase in energy costs. Therefore, the price of EUAs may be affected by variations in oil prices if the greenhouse gases discharged by the consumption of alternative energy are less than that of fossil fuels. There are no previous studies that investigated these relationships. In this study, we analyzed eight types of EUAs EUA05 to EUA12 with a time series daily data set during 2005-2007 collected from a European Climate Exchange time series data set. Differentiations in these eight types were redemption period. We used the New York Mercantile Exchange light sweet crude price as an oil price. From our examination, we found that only the EUA06 and EUA07 types of EUAs Granger-cause oil prices and vice versa and other six types of EUAs do not Granger-cause oil price. These results imply that the earlier redemption period types of EUAs are more sensitive to oil price. In employing the impulse response function, the results showed that a shock to oil price has a slightly positive effect on all types of EUAs for a very short period. On the other hand, we found that a shock to price of EUA has a slightly negative effect on oil price following a positive effect in only EUA06 and EUA07 types. Therefore, these results imply that fluctuations in EUAs prices and oil prices have little effect on each other. Lastly, we did not consider the substitute energy prices in this study, so we plan to include the prices of coal and natural gas in future analyses.

Novel non-linear elastic structural analysis with generalised transverse element loads using a refined finite element

In the finite element modelling of structural frames, external loads such as wind loads, dead loads and imposed loads usually act along the elements rather than at the nodes only. Conventionally, when an element is subjected to these general transverse element loads, they are usually converted to nodal forces acting at the ends of the elements by either lumping or consistent load approaches. In addition, it is especially important for an element subjected to the first- and second-order elastic behaviour, to which the steel structure is critically prone to; in particular the thin-walled steel structures, when the stocky element section may be generally critical to the inelastic behaviour. In this sense, the accurate first- and second-order elastic displacement solutions of element load effect along an element is vitally crucial, but cannot be simulated using neither numerical nodal nor consistent load methods alone, as long as no equilibrium condition is enforced in the finite element formulation, which can inevitably impair the structural safety of the steel structure particularly. It can be therefore regarded as a unique element load method to account for the element load nonlinearly. If accurate displacement solution is targeted for simulating the first- and second-order elastic behaviour on an element on the basis of sophisticated non-linear element stiffness formulation, the numerous prescribed stiffness matrices must indispensably be used for the plethora of specific transverse element loading patterns encountered. In order to circumvent this shortcoming, the present paper proposes a numerical technique to include the transverse element loading in the non-linear stiffness formulation without numerous prescribed stiffness matrices, and which is able to predict structural responses involving the effect of first-order element loads as well as the second-order coupling effect between the transverse load and axial force in the element. This paper shows that the principle of superposition can be applied to derive the generalized stiffness formulation for element load effect, so that the form of the stiffness matrix remains unchanged with respect to the specific loading patterns, but with only the magnitude of the loading (element load coefficients) being needed to be adjusted in the stiffness formulation, and subsequently the non-linear effect on element loadings can be commensurate by updating the magnitude of element load coefficients through the non-linear solution procedures. In principle, the element loading distribution is converted into a single loading magnitude at mid-span in order to provide the initial perturbation for triggering the member bowing effect due to its transverse element loads. This approach in turn sacrifices the effect of element loading distribution except at mid-span. Therefore, it can be foreseen that the load-deflection behaviour may not be as accurate as those at mid-span, but its discrepancy is still trivial as proved. This novelty allows for a very useful generalised stiffness formulation for a single higher-order element with arbitrary transverse loading patterns to be formulated. Moreover, another significance of this paper is placed on shifting the nodal response (system analysis) to both nodal and element response (sophisticated element formulation). For the conventional finite element method, such as the cubic element, all accurate solutions can be only found at node. It means no accurate and reliable structural safety can be ensured within an element, and as a result, it hinders the engineering applications. The results of the paper are verified using analytical stability function studies, as well as with numerical results reported by independent researchers on several simple frames.

A nationwide web-based automated system for early outbreak detection and rapid response in China

Timely reporting, effective analyses and rapid distribution of surveillance data can assist in detecting the aberration of disease occurrence and further facilitate a timely response. In China, a new nationwide web-based automated system for outbreak detection and rapid response was developed in 2008. The China Infectious Disease Automated-alert and Response System (CIDARS) was developed by the Chinese Center for Disease Control and Prevention based on the surveillance data from the existing electronic National Notifiable Infectious Diseases Reporting Information System (NIDRIS) started in 2004. NIDRIS greatly improved the timeliness and completeness of data reporting with real time reporting information via the Internet. CIDARS further facilitates the data analysis, aberration detection, signal dissemination, signal response and information communication needed by public health departments across the country. In CIDARS, three aberration detection methods are used to detect the unusual occurrence of 28 notifiable infectious diseases at the county level and to transmit that information either in real-time or on a daily basis. The Internet, computers and mobile phones are used to accomplish rapid signal generation and dissemination, timely reporting and reviewing of the signal response results. CIDARS has been used nationwide since 2008; all Centers for Disease Control and Prevention (CDC) in China at the county, prefecture, provincial and national levels are involved in the system. It assists with early outbreak detection at the local level and prompts reporting of unusual disease occurrences or potential outbreaks to CDCs throughout the country.

Evaluation of the matabolic response of osteosarcoma cells to conventional and new anticancer drugs by NMR metabolomics

The main scope of this work was to evaluate the metabolic effects of anticancer agents (three conventional and one new) in osteosarcoma (OS) cells and osteoblasts, by measuring alterations in the metabolic profile of cells by nuclear magnetic resonance (NMR) spectroscopy metabolomics. Chapter 1 gives a theoretical framework of this work, beginning with the main metabolic characteristics that globally describe cancer as well as the families and mechanisms of action of drugs used in chemotherapy. The drugs used nowadays to treat OS are also presented, together with the Palladium(II) complex with spermine, Pd2Spm, potentially active against cancer. Then, the global strategy for cell metabolomics is explained and the state of the art of metabolomic studies that analyze the effect of anticancer agents in cells is presented. In Chapter 2, the fundamentals of the analytical techniques used in this work, namely for biological assays, NMR spectroscopy and multivariate and statistical analysis of the results are described. A detailed description of the experimental procedures adopted throughout this work is given in Chapter 3. The biological and analytical reproducibility of the metabolic profile of MG-63 cells by high resolution magic angle spinning (HRMAS) NMR is evaluated in Chapter 4. The metabolic impact of several factors (cellular integrity, spinning rate, temperature, time and acquisition parameters) on the 1H HRMAS NMR spectral profile and quality is analysed, enabling the definition of the best acquisition parameters for further experiments. The metabolic consequences of increasing number of passages in MG-63 cells as well as the duration of storage are also investigated. Chapter 5 describes the metabolic impact of drugs conventionally used in OS chemotherapy, through NMR metabolomics studies of lysed cells and aqueous extracts analysis. The results show that MG-63 cells treated with cisplatin (cDDP) undergo a strong up-regulation of lipid contents, alterations in phospholipid constituents (choline compounds) and biomarkers of DNA degradation, all associated with cell death by apoptosis. Cells exposed to doxorubicin (DOX) or methotrexate (MTX) showed much slighter metabolic changes, without any relevant alteration in lipid contents. However, metabolic changes associated with altered Krebs cycle, oxidative stress and nucleotides metabolism were detected and were tentatively interpreted at the light of the known mechanisms of action of these drugs. The metabolic impact of the exposure of MG-63 cells and osteoblasts to cDDP and the Pd2Spm complex is described in Chapter 6. Results show that, despite the ability of the two agents to bind DNA, the metabolic consequences that arise from exposure to them are distinct, namely in what concerns to variation in lipid contents (absent for Pd2Spm). Apoptosis detection assays showed that, differently from what was seen for MG-63 cells treated with cDDP, the decreased number of living cells upon exposure to Pd2Spm was not due to cell death by apoptosis or necrosis. Moreover, the latter agent induces more marked alterations in osteoblasts than in cancer cells, while the opposite seemed to occur upon cDDP exposure. Nevertheless, the results from MG-63 cells exposure to combination regimens with cDDP- or Pd2Spm-based cocktails, described in Chapter 7, revealed that, in combination, the two agents induce similar metabolic responses, arising from synergy mechanisms between the tested drugs. Finally, the main conclusions of this thesis are summarized in Chapter 8, and future perspectives in the light of this work are presented.