基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

中国地下工程氡污染健康影响与防护研究-以内蒙古自治区为例

由于对人体健康的危害，氡越来越受到世界各国的广泛关注。经济的发展，使我国地下公共场所大量涌现，但其空气中氡浓度的大规模调查及系统研究还属空白。本论文采取瞬时法（美国SUN公司的1027型连续测氰仪）与累积法（CR一39SSNTD）相结合的测量方法，对我国23个城市的地下工程进行了为期一年的监测，获得了实测数据700多个。在监测的同时，收集了每个测点的各种相关情况资料。论文以内蒙古自治区的呼和浩特、包头、集宁、鄂尔多斯四个城市为例，从环境地球化学角度出发，在研究了各个测点的各种情况之后，探讨了氧在不同类型工程和不同城市的浓度规律。对长期有人在内工作的地下工程内氛所致工作人员的有效剂量当量进行了估算。最后对氨异常的测点进行了原因分析，并提出了地下工程内氰污染的治理和防范措施。本次研究主要得出以下结论：1、地下工程中氰浓度的季节变化规律为：从总平均值看，夏高冬低。与前人对坑道的调查结果一致，与地上氨浓度的变化规律相反。从两季浓度差来看，两季浓度差从大到小依次为：隧道类工程）单建式铜室类工程＞结建式嗣室类工程。2、总体来看内蒙三类工程的氰浓度值：隧道类工程）结建式铜室类工程）单建式铜室类工程。不管是单季还是年平均值，各城市三类工程的氧浓度都遵循隧道类．工程远大于结建式铜室工程和单建式铜室工程的规律。通过对结建式桐室工程和单建式铜室工程的对比分析，发现在地下工程内工程上部有无建筑物阻隔，并不是影响氨浓度高低的重要因素，即氰的纵向扩散作用小于横向扩散作用。3、呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度年平均值分别为：190Bq／m3、700Bq／m3、190Bq／m3、270Bq／m3。去掉仅有一季结果的侧点值，包头、鄂尔多斯的氨浓度年平均值为：62oBq／m3、250Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度单季值的范围分别为：31-830Bq／m3、57-3500Bq／m3、26一65oBq／m，、64一69OBq／m'，；年平均值的范围分别为：39一580Bq／m：'、94一210oBq／m：，、35-460Bq／m3、110-390Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市分别有30.0%、44.4%、33.3%、50.0%的工程年平均氧浓度在200Bq／时以上。四个城市氦浓度的普遍水平：包头）鄂尔多斯）集宁）呼和浩特。4、对所有单建式铜室类地下工程和结建式铜室类工程内工作人员所受剂量当量进行了估算。单建式地下工程内工作人员所受剂量当量的范围是1.247mSv-6．SOgmsv，平均值是2.610Sv。结建式地下工程内工作人员所受剂量当量的范围是1．195Sv-8.531mSv，平均值是2.572mSv。两类地下工程内工作人员所受氨照射剂量当量都很高，有一半的工程剂量当量是地面上公众人员剂量当量0.924mSv的两倍以上。5、不管是从浓度角度还是从剂量估算角度分析，地下工程内的氛水平都比较偏高，需要采取一定防护措施，对现有状况进行改善。6、通风是地下工程防治氡污染最简便易行、见效最快的措施。使用防氡涂料是最经济的措施。

Cartilha dia de campo: suinocultura e meio ambiente: termo de compromisso de ajustamento de conduta da suinocultura.

Apresenta um breve resumo dos principais aspectos que fora abordados nas dez estações que compuseram a programação do dia de campo. Além disso, ela foi elaborada levando-se em considerações as dúvidas mais freqüentes que os suinocultores têm demonstrado em relação aos aspectos legais e técnicos do TAC da suinocultura.

Métodos de avaliação da estabilidade implantar

Projeto de Pós-Graduação/Dissertação apresentado à Universidade Fernando Pessoa como parte dos requisitos para obtenção do grau de Mestre em Medicina Dentária

Application of signal detection methods to fisheries management

The abundance of many commercially important fish stocks are declining and this has led to widespread concern on the performance of traditional approach in fisheries management. Quantitative models are used for obtaining estimates of population abundance and the management advice is based on annual harvest levels (TAC), where only a certain amount of catch is allowed from specific fish stocks. However, these models are data intensive and less useful when stocks have limited historical information. This study examined whether empirical stock indicators can be used to manage fisheries. The relationship between indicators and the underlying stock abundance is not direct and hence can be affected by disturbances that may account for both transient and persistent effects. Methods from Statistical Process Control (SPC) theory such as the Cumulative Sum (CUSUM) control charts are useful in classifying these effects and hence they can be used to trigger management response only when a significant impact occurs to the stock biomass. This thesis explores how empirical indicators along with CUSUM can be used for monitoring, assessment and management of fish stocks. I begin my thesis by exploring various age based catch indicators, to identify those which are potentially useful in tracking the state of fish stocks. The sensitivity and response of these indicators towards changes in Spawning Stock Biomass (SSB) showed that indicators based on age groups that are fully selected to the fishing gear or Large Fish Indicators (LFIs) are most useful and robust across the range of scenarios considered. The Decision-Interval (DI-CUSUM) and Self-Starting (SS-CUSUM) forms are the two types of control charts used in this study. In contrast to the DI-CUSUM, the SS-CUSUM can be initiated without specifying a target reference point (‘control mean’) to detect out-of-control (significant impact) situations. The sensitivity and specificity of SS-CUSUM showed that the performances are robust when LFIs are used. Once an out-of-control situation is detected, the next step is to determine how much shift has occurred in the underlying stock biomass. If an estimate of this shift is available, they can be used to update TAC by incorporation into Harvest Control Rules (HCRs). Various methods from Engineering Process Control (EPC) theory were tested to determine which method can measure the shift size in stock biomass with the highest accuracy. Results showed that methods based on Grubb’s harmonic rule gave reliable shift size estimates. The accuracy of these estimates can be improved by monitoring a combined indicator metric of stock-recruitment and LFI because this may account for impacts independent of fishing. The procedure of integrating both SPC and EPC is known as Statistical Process Adjustment (SPA). A HCR based on SPA was designed for DI-CUSUM and the scheme was successful in bringing out-of-control fish stocks back to its in-control state. The HCR was also tested using SS-CUSUM in the context of data poor fish stocks. Results showed that the scheme will be useful for sustaining the initial in-control state of the fish stock until more observations become available for quantitative assessments.

A new proposition on the martingale representation theorem and on the approximate hedging of contingent claim in mean-variance criterion

In this work we revisit the problem of the hedging of contingent claim using mean-square criterion. We prove that in incomplete market, some probability measure can be identified so that becomes -martingale under .This is in fact a new proposition on the martingale representation theorem. The new results also identify a weight function that serves to be an approximation to the Radon-Nikodým derivative of the unique neutral martingale measure.

Level of beta-adrenergic receptor kinase 1 inhibition determines degree of cardiac dysfunction after chronic pressure overload-induced heart failure.

BACKGROUND: Heart failure is characterized by abnormalities in beta-adrenergic receptor (betaAR) signaling, including increased level of myocardial betaAR kinase 1 (betaARK1). Our previous studies have shown that inhibition of betaARK1 with the use of the Gbetagamma sequestering peptide of betaARK1 (betaARKct) can prevent cardiac dysfunction in models of heart failure. Because inhibition of betaARK activity is pivotal for amelioration of cardiac dysfunction, we investigated whether the level of betaARK1 inhibition correlates with the degree of heart failure. METHODS AND RESULTS: Transgenic (TG) mice with varying degrees of cardiac-specific expression of betaARKct peptide underwent transverse aortic constriction (TAC) for 12 weeks. Cardiac function was assessed by serial echocardiography in conscious mice, and the level of myocardial betaARKct protein was quantified at termination of the study. TG mice showed a positive linear relationship between the level of betaARKct protein expression and fractional shortening at 12 weeks after TAC. TG mice with low betaARKct expression developed severe heart failure, whereas mice with high betaARKct expression showed significantly less cardiac deterioration than wild-type (WT) mice. Importantly, mice with a high level of betaARKct expression had preserved isoproterenol-stimulated adenylyl cyclase activity and normal betaAR densities in the cardiac membranes. In contrast, mice with low expression of the transgene had marked abnormalities in betaAR function, similar to the WT mice. CONCLUSIONS: These data show that the level of betaARK1 inhibition determines the degree to which cardiac function can be preserved in response to pressure overload and has important therapeutic implications when betaARK1 inhibition is considered as a molecular target.

On-board Robotic Multi-pinhole SPECT System for Region-of-interest (ROI) Imaging

On-board image guidance, such as cone-beam CT (CBCT) and kV/MV 2D imaging, is essential in many radiation therapy procedures, such as intensity modulated radiotherapy (IMRT) and stereotactic body radiation therapy (SBRT). These imaging techniques provide predominantly anatomical information for treatment planning and target localization. Recently, studies have shown that treatment planning based on functional and molecular information about the tumor and surrounding tissue could potentially improve the effectiveness of radiation therapy. However, current on-board imaging systems are limited in their functional and molecular imaging capability. Single Photon Emission Computed Tomography (SPECT) is a candidate to achieve on-board functional and molecular imaging. Traditional SPECT systems typically take 20 minutes or more for a scan, which is too long for on-board imaging. A robotic multi-pinhole SPECT system was proposed in this dissertation to provide shorter imaging time by using a robotic arm to maneuver the multi-pinhole SPECT system around the patient in position for radiation therapy.

A 49-pinhole collimated SPECT detector and its shielding were designed and simulated in this work using the computer-aided design (CAD) software. The trajectories of robotic arm about the patient, treatment table and gantry in the radiation therapy room and several detector assemblies such as parallel holes, single pinhole and 49 pinholes collimated detector were investigated. The rail mounted system was designed to enable a full range of detector positions and orientations to various crucial treatment sites including head and torso, while avoiding collision with linear accelerator (LINAC), patient table and patient.

An alignment method was developed in this work to calibrate the on-board robotic SPECT to the LINAC coordinate frame and to the coordinate frames of other on-board imaging systems such as CBCT. This alignment method utilizes line sources and one pinhole projection of these line sources. The model consists of multiple alignment parameters which maps line sources in 3-dimensional (3D) space to their 2-dimensional (2D) projections on the SPECT detector. Computer-simulation studies and experimental evaluations were performed as a function of number of line sources, Radon transform accuracy, finite line-source width, intrinsic camera resolution, Poisson noise and acquisition geometry. In computer-simulation studies, when there was no error in determining angles (α) and offsets (ρ) of the measured projections, the six alignment parameters (3 translational and 3 rotational) were estimated perfectly using three line sources. When angles (α) and offsets (ρ) were provided by Radon transform, the estimation accuracy was reduced. The estimation error was associated with rounding errors of Radon transform, finite line-source width, Poisson noise, number of line sources, intrinsic camera resolution and detector acquisition geometry. The estimation accuracy was significantly improved by using 4 line sources rather than 3 and also by using thinner line-source projections (obtained by better intrinsic detector resolution). With 5 line sources, median errors were 0.2 mm for the detector translations, 0.7 mm for the detector radius of rotation, and less than 0.5° for detector rotation, tilt and twist. In experimental evaluations, average errors relative to a different, independent registration technique were about 1.8 mm for detector translations, 1.1 mm for the detector radius of rotation (ROR), 0.5° and 0.4° for detector rotation and tilt, respectively, and 1.2° for detector twist.

Simulation studies were performed to investigate the improvement of imaging sensitivity and accuracy of hot sphere localization for breast imaging of patients in prone position. A 3D XCAT phantom was simulated in the prone position with nine hot spheres of 10 mm diameter added in the left breast. A no-treatment-table case and two commercial prone breast boards, 7 and 24 cm thick, were simulated. Different pinhole focal lengths were assessed for root-mean-square-error (RMSE). The pinhole focal lengths resulting in the lowest RMSE values were 12 cm, 18 cm and 21 cm for no table, thin board, and thick board, respectively. In both no table and thin board cases, all 9 hot spheres were easily visualized above background with 4-minute scans utilizing the 49-pinhole SPECT system while seven of nine hot spheres were visible with the thick board. In comparison with parallel-hole system, our 49-pinhole system shows reduction in noise and bias under these simulation cases. These results correspond to smaller radii of rotation for no-table case and thinner prone board. Similarly, localization accuracy with the 49-pinhole system was significantly better than with the parallel-hole system for both the thin and thick prone boards. Median localization errors for the 49-pinhole system with the thin board were less than 3 mm for 5 of 9 hot spheres, and less than 6 mm for the other 4 hot spheres. Median localization errors of 49-pinhole system with the thick board were less than 4 mm for 5 of 9 hot spheres, and less than 8 mm for the other 4 hot spheres.

Besides prone breast imaging, respiratory-gated region-of-interest (ROI) imaging of lung tumor was also investigated. A simulation study was conducted on the potential of multi-pinhole, region-of-interest (ROI) SPECT to alleviate noise effects associated with respiratory-gated SPECT imaging of the thorax. Two 4D XCAT digital phantoms were constructed, with either a 10 mm or 20 mm diameter tumor added in the right lung. The maximum diaphragm motion was 2 cm (for 10 mm tumor) or 4 cm (for 20 mm tumor) in superior-inferior direction and 1.2 cm in anterior-posterior direction. Projections were simulated with a 4-minute acquisition time (40 seconds per each of 6 gates) using either the ROI SPECT system (49-pinhole) or reference single and dual conventional broad cross-section, parallel-hole collimated SPECT. The SPECT images were reconstructed using OSEM with up to 6 iterations. Images were evaluated as a function of gate by profiles, noise versus bias curves, and a numerical observer performing a forced-choice localization task. Even for the 20 mm tumor, the 49-pinhole imaging ROI was found sufficient to encompass fully usual clinical ranges of diaphragm motion. Averaged over the 6 gates, noise at iteration 6 of 49-pinhole ROI imaging (10.9 µCi/ml) was approximately comparable to noise at iteration 2 of the two dual and single parallel-hole, broad cross-section systems (12.4 µCi/ml and 13.8 µCi/ml, respectively). Corresponding biases were much lower for the 49-pinhole ROI system (3.8 µCi/ml), versus 6.2 µCi/ml and 6.5 µCi/ml for the dual and single parallel-hole systems, respectively. Median localization errors averaged over 6 gates, for the 10 mm and 20 mm tumors respectively, were 1.6 mm and 0.5 mm using the ROI imaging system and 6.6 mm and 2.3 mm using the dual parallel-hole, broad cross-section system. The results demonstrate substantially improved imaging via ROI methods. One important application may be gated imaging of patients in position for radiation therapy.

A robotic SPECT imaging system was constructed utilizing a gamma camera detector (Digirad 2020tc) and a robot (KUKA KR150-L110 robot). An imaging study was performed with a phantom (PET CT PhantomTM), which includes 5 spheres of 10, 13, 17, 22 and 28 mm in diameter. The phantom was placed on a flat-top couch. SPECT projections were acquired with a parallel-hole collimator and a single-pinhole collimator both without background in the phantom, and with background at 1/10th the sphere activity concentration. The imaging trajectories of parallel-hole and pinhole collimated detectors spanned 180 degrees and 228 degrees respectively. The pinhole detector viewed a 14.7 cm-diameter common volume which encompassed the 28 mm and 22 mm spheres. The common volume for parallel-hole was a 20.8-cm-diameter cylinder which encompassed all five spheres in the phantom. The maneuverability of the robotic system was tested by navigating the detector to trace the flat-top table while avoiding collision with the table and maintaining the closest possible proximity to the common volume. For image reconstruction, detector trajectories were described by radius-of-rotation and detector rotation angle θ. These reconstruction parameters were obtained from the robot base and tool coordinates. The robotic SPECT system was able to maneuver the parallel-hole and pinhole collimated SPECT detectors in close proximity to the phantom, minimizing impact of the flat-top couch on detector to center-of-rotation (COR) distance. In no background case, all five spheres were visible in the reconstructed parallel-hole and pinhole images. In with background case, three spheres of 17, 22 and 28 mm diameter were readily observed with the parallel-hole imaging, and the targeted spheres (22 and 28 mm diameter) were readily observed in the pinhole ROI imaging.

In conclusion, the proposed on-board robotic SPECT can be aligned to LINAC/CBCT with a single pinhole projection of the line-source phantom. Alignment parameters can be estimated using one pinhole projection of line sources. This alignment method may be important for multi-pinhole SPECT, where relative pinhole alignment may vary during rotation. For single pinhole and multi-pinhole SPECT imaging onboard radiation therapy machines, the method could provide alignment of SPECT coordinates with those of CBCT and the LINAC. In simulation studies of prone breast imaging and respiratory-gated lung imaging, the 49-pinhole detector showed better tumor contrast recovery and localization in a 4-minute scan compared to parallel-hole detector. On-board SPECT could be achieved by a robot maneuvering a SPECT detector about patients in position for radiation therapy on a flat-top couch. The robot inherent coordinate frames could be an effective means to estimate detector pose for use in SPECT image reconstruction.

Opportunistic Control Over Shared Wireless Channels

© 2015 IEEE.We consider a wireless control architecture with multiple control loops over a shared wireless medium. A scheduler observes the random channel conditions that each control system experiences over the shared medium and opportunistically selects systems to transmit at a set of non-overlapping frequencies. The transmit power of each system also adapts to channel conditions and determines the probability of successfully receiving and closing the loop. We formulate the optimal design of channel-aware scheduling and power allocation that minimize the total power consumption while meeting control performance requirements for all systems. In particular, it is required that for each control system a given Lyapunov function decreases at a specified rate in expectation over the random channel conditions. We develop an offline algorithm to find the optimal communication design, as well as an online protocol which selects scheduling and power variables based on a random observed channel sequence and converges almost surely to the optimal operating point. Simulations illustrate the power savings of our approach compared to other non-channel-aware schemes.

CB1 cannabinoid receptor deficiency promotes cardiac remodeling induced by pressure overload in mice

Background: The endocannabinoid system is known to play a role in regulating myocardial contractility, but the influence of cannabinoid receptor 1 (CB1) deficiency on chronic heart failure (CHF) remains unclear. In this study we attempted to investigate the effect of CB1 deficiency on CHF induced by pressure overload and the possible mechanisms involved. Methods and results: A CHF model was created by transverse aortic constriction (TAC) in both CB1 knockout mice and wild-type mice. CB1 knockout mice showed a marked increase of mortality due to CHF from 4 to 8 weeks after TAC (p = 0.021). Five weeks after TAC, in contrast to wild-type mice, CB1 knockout mice had a higher left ventricular (LV) end-diastolic pressure, lower rate of LV pressure change (± dp/dt max), lower LV contractility index, and a larger heart weight to body weight ratio and lung weight to body weight ratio compared with wild-type mice (all p < 0.05-0.001). Phosphorylation of the epidermal growth factor receptor (EGFR) and mitogen-activated protein kinases (P38 and ERK) was higher in CB1 knockout mice than that in wild-type mice. In cultured neonatal rat cardiomyocytes, a CB1 agonist reduced cAMP production stimulated by isoproterenol or forskolin, and suppressed phosphorylation of the EGFR, P38, and ERK, while the inhibitory effect of a CB1 agonist on EGFR phosphorylation was abrogated by CB1 knockdown. Conclusion: These findings indicate that cannabinoid receptor 1 inactivation promotes cardiac remodeling by enhancing the activity of the epidermal growth factor receptor and mitogen-activated protein kinases. © 2012 Elsevier Ireland Ltd.

A Fast Nonlinear Model Identification Method

The identification of nonlinear dynamic systems using linear-in-the-parameters models is studied. A fast recursive algorithm (FRA) is proposed to select both the model structure and to estimate the model parameters. Unlike orthogonal least squares (OLS) method, FRA solves the least-squares problem recursively over the model order without requiring matrix decomposition. The computational complexity of both algorithms is analyzed, along with their numerical stability. The new method is shown to require much less computational effort and is also numerically more stable than OLS.

Blogs as virtual environments for teaching and learning in Higher Education

The transformation of the Information and Communication Technologies (ICT) into knowledge and learning technologies is increasingly becoming a matter of concern. Teaching settings associated with the use of blogs in Higher Education are presented in this paper, proceeding from an innovative learning experience of projects carried out by a group of professors between 2009 and 2013. Both, teachers and students who took part in the subjects that implemented the blog, considered it as helpful resource to create a virtual and learning-teaching environment due to the multiple potentialities it offers. Among some of these potentialities, some stand out: it makes easier the access to knowledge, promotes a more active and reflective learning, expands the social experience of learning, provides evidence about the students’ progress which helps to reorient the teaching-learning process, and encourages the critical judgment. Nevertheless, several problems related with the students’ participation and the teacher’s blog management have been identified.