In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can be...In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling (LWD), in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
According to the positive correlation of coal ash content and natural gamma, using a new coal core reposition method, which is ordered by global and local extreme, coal samples from medium-thickness seam are reasonabl...According to the positive correlation of coal ash content and natural gamma, using a new coal core reposition method, which is ordered by global and local extreme, coal samples from medium-thickness seam are reasonably located. Inte- grated the data of coal macrostructure characteristics, coal petrography analysis and coal gas production test, it studies the rela- tionship between coalbody structure and amplitude variation of different well logging data, and the tectonic coal recognition method with well logging data in fresh-water mud invasion. The results show that: the anomalous response of natural gamma ray, neutron, density and apparent resistivity does not reflect the coalbody structure type. In fresh-water drilling mud invasion, using the crossplot technique of dual-lateral, RXO resistivity response and the coalbody structure can classify granulated coal accurately; the proposed method is of good practicability and high reliability.展开更多
In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The...In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The result shows that the resistivity logs become smoother and lower as the borehole diameter increases, the increase of the contrast between mud resistivity and formation resistivity induce the logs to be more pointed. When the formation thickness is less than lm, the two-peak on the logs for resistive invasion vanished, and for thickness between 1 m and 4 m, the form of logs does not vary significantly. If the formation thickness is greater than 4 m, a platform appears on the logs at the middle of the formation. The thinner the invaded zone is, the more obvious the invasion feature on the laterologs is. For thick invaded zone the form of logs tend to be that of an uninvaded resistive formation. The form and amplitude of logs depend on the resistivity contrast between invaded zone, uninvaded formation and adjacentlayers.展开更多
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML (H-PML) combining the normal absorption factor of convolutional PML (C-PML) with tangential absorption factor of Mutiaxial PML (M-PML). The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling (LWD), in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
文摘According to the positive correlation of coal ash content and natural gamma, using a new coal core reposition method, which is ordered by global and local extreme, coal samples from medium-thickness seam are reasonably located. Inte- grated the data of coal macrostructure characteristics, coal petrography analysis and coal gas production test, it studies the rela- tionship between coalbody structure and amplitude variation of different well logging data, and the tectonic coal recognition method with well logging data in fresh-water mud invasion. The results show that: the anomalous response of natural gamma ray, neutron, density and apparent resistivity does not reflect the coalbody structure type. In fresh-water drilling mud invasion, using the crossplot technique of dual-lateral, RXO resistivity response and the coalbody structure can classify granulated coal accurately; the proposed method is of good practicability and high reliability.
文摘In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The result shows that the resistivity logs become smoother and lower as the borehole diameter increases, the increase of the contrast between mud resistivity and formation resistivity induce the logs to be more pointed. When the formation thickness is less than lm, the two-peak on the logs for resistive invasion vanished, and for thickness between 1 m and 4 m, the form of logs does not vary significantly. If the formation thickness is greater than 4 m, a platform appears on the logs at the middle of the formation. The thinner the invaded zone is, the more obvious the invasion feature on the laterologs is. For thick invaded zone the form of logs tend to be that of an uninvaded resistive formation. The form and amplitude of logs depend on the resistivity contrast between invaded zone, uninvaded formation and adjacentlayers.