Multi-component induction resistivity and sonic measurements acquired in high-angle wells can be strongly influenced by shoulder-bed effects, anisotropy resulting from sand-shale laminations, and presence of mud-filtrate invasion. Understanding the corresponding biasing effects aids in the interpretation of resistivity and sonic measurements and subsequently leads to more accurate and reliable formation evaluation. This dissertation describes numerical simulation studies examining the effects on multi-component induction and sonic measurements in a variety of complex formation models. Subsequently, a joint inversion scheme is presented that combines resistivity and sonic measurements to estimate in situ petrophysical and elastic properties in the presence of mud-filtrate invasion. To facilitate the simulation study of multi-component induction logs, I develop a new finite-difference algorithm for the numerical simulation of frequency-domain electromagnetic borehole measurements. The algorithm~uses a coupled scalar-vector potential formulation for arbitrary three-dimensional inhomogeneous and electrically anisotropic media. Simulations show that shoulder-bed anisotropy: enhances shoulder-bed effects across sand layers; and impacts invasion sensitivities to significantly alter the assessment of invasion in terms of invaded- and virgin-zone resistivities, radial length, and front shape. For the simulation study of sonic logs, I develop a three-dimensional, finite-difference time-domain algorithm that models elastic wave propagation in a fluid-filled borehole. Simulations show that presence of anisotropy not only alters the degree of dispersion observed in flexural and Stoneley waves, but also alters their responses to invasion. In addition, presence of a dipping shoulder bed can significantly distort flexural dispersion, making it difficult to identify the low frequency asymptote corresponding to formation shear wave velocity. Lastly, I consider a radial one-dimensional model in the development of a joint resistivity and sonic inversion algorithm. This scheme simultaneously inverts array-induction apparent conductivities and sonic flexural and Stoneley dispersions for the rock's elastic moduli and water saturation in the presence of mud-filtrate invasion. Inversions are performed on numerically simulated data for a variety of models reflecting soft and hard rock formations with presence of water- and oil-based mud-filtrate invasion. Results show the estimated invasion profiles display excellent agreement with the true models, and the elastic moduli are estimated to within a few percent of the true values. / text

Pulsed-neutron borehole measurements involve a physical process in which a source emits energetic neutrons that lose energy upon collisions with formation nuclei, and are eventually captured by a nucleus to form a heavier, excited state. The excited nucleus decays to its ground state by the emission of gamma rays. Both thermal-neutron and gamma-ray populations decay with time at a rate defined by Sigma, which is a nuclear property that quantifies a material’s ability to capture thermal neutrons. The large contrast in Sigma between hydrocarbon and salty connate water enables calculations of water saturation directly from pulsed-neutron measurements. Sigma logs have proven useful in the assessment of thinly bedded formations because they exhibit a small volume of investigation, and have been deemed superior to resistivity logs in the petrophysical evaluation of carbonate formations. The recognized potential of Sigma logs in formation evaluation initiated the development of multi-detector Logging-While-Drilling (LWD) Sigma measurements. These measurements are acquired using one thermal-neutron and two gamma-ray detectors at different spacings from the source. Such a design is aimed at providing distinct radial depths of investigation to detect filtrate invasion in the near-wellbore zone. Despite their formation-evaluation potential, multi-detector time-decay measurements commonly remain affected by invasion, shoulder-bed, and well-deviation effects. The purpose of this dissertation is to develop a fast-forward simulation method to reproduce multi-detector time decays and combine the method with inversion techniques to improve the petrophysical interpretation of LWD Sigma measurements. First-order perturbation theory and a library of pre-calculated Monte Carlo detector-specific sensitivity functions and time decays are used to numerically simulate borehole Sigma measurements in realistic logging environments. The new simulation method is one hundred thousand times faster than rigorous Monte Carlo calculations and remains within two capture units of disparity. Next, the fast-forward simulation method is embedded within inversion algorithms to estimate layer-by-layer radial length of invasion and formation Sigma corrected for shallow invasion, shoulder-bed, and well-deviation effects. Both fast-forward and inverse modeling algorithms are benchmarked against laboratory and synthetic time decays. The improvement of formation Sigma obtained with inversion-based interpretation leads to an improvement in the estimation of Sigma-derived water saturation. Likewise, the estimated radial length of invasion is combined with neutron and density measurements to correct the latter for invasion effects. Results indicate that the inversion-based interpretation method is well suited for the evaluation of high-porosity formations invaded by salty mud filtrate. Inversion-based interpretation of field LWD time decays enables the estimation of lower values of water saturation when compared to conventional Sigma interpretation or resistivity methods. Estimated values of water saturation are as much as fifty percent lower than predicted by conventional interpretation of Sigma logs in the case of measurements affected by shoulder-bed effects, and as much as one hundred percent lower than predicted by the conventional interpretation method for measurements additionally affected by salty filtrate invasion. The key attributes of the combined petrophysical interpretation of multi-detector Sigma, neutron, and density measurements developed in this dissertation are that it explicitly enforces the physics of all nuclear measurements, honors the pressure and temperature dependency of reservoir fluid nuclear properties, and takes into account a-priori information such as mud-filtrate salinity, connate-water salinity, and bed-boundary locations. / text