首页 > 最新文献

SPWLA 62nd Annual Online Symposium Transactions最新文献

英文 中文
A NOVEL METHOD FOR FORMATION DENSITY MEASUREMENT IN CASED WELLS 套管井地层密度测量新方法
Pub Date : 2021-05-17 DOI: 10.30632/spwla-2021-0065
Xinguang Wang, Dong Li, Lei Zhang, Feng Zhang, Wen-hao Wang, Huiping Gao, China Oilfield Services Limited, Xi 'an HuiNeng Electronic Equipment Co. Ltd
Density logging is a fundamental logging method of open-hole logging series. In cased wells, the formation is usually treated as a cylindrical multiple-layer model including borehole, casing pipe, cement, and original formation, producing more complicated logging response compared to open-hole wells. Several studies have indicated that it is feasible to measure formation density using a density logging tool in cased well under certain conditions, of which the key problem is to determine casing thickness, cement thickness and cement density. This study proposes a novel method to evaluate formation density behind casing using a three-detector density tool. Considering that the backscattered gamma is easily affected by the density behind the casing, especially when the casing thickness is small, a designed ratio of high and low energy window counts of the scattered gamma spectrum is used to calculate casing thickness, which improves the lower limit of the detection. After obtaining casing thickness, the cylindrical multiple-layer model of cased well can be simplified to a cylindrical three-layer model, consisting of cased borehole, cement, and the original formation. Subsequently, the analytical expressions of the thickness of the middle layer (cement) as well as the density of the outer layer (original formation) can be derived based on the logging geometric factor concept, G, which borrowing from electric logging. Consequently, the cement thickness and original formation density could be calculated independently. To demonstrate the feasibility of the proposed method, experimental data are measured in a group of calibration wells. The data processing results show that the proposed analytical model fit well with the experimental data. In addition, the application conditions of the novel method are discussed based on experimental data. Finally, several log examples illustrate that under suitable casing and cement conditions, the formation density calculated with the proposed method in cased borehole is in good agreement with the corresponding open-hole density.
密度测井是裸眼测井系列的一种基本测井方法。在套管井中,地层通常被视为一个圆柱形多层模型,包括井眼、套管、水泥和原始地层,与裸眼井相比,产生更复杂的测井响应。多项研究表明,在一定条件下,利用密度测井仪测量套管井地层密度是可行的,其中关键问题是确定套管厚度、水泥厚度和水泥密度。本研究提出了一种利用三探测器密度工具评估套管后地层密度的新方法。考虑到后向散射伽马容易受到套管后密度的影响,特别是当套管厚度较小时,采用设计的散射伽马能谱高低能窗计数比来计算套管厚度,提高了检测下限。在得到套管厚度后,将套管井的圆柱多层模型简化为由套管井眼、水泥和原始地层组成的圆柱三层模型。随后,借鉴电测井的测井几何因子G概念,推导出了中间层(水泥)厚度和外层(原始地层)密度的解析表达式。因此,可以独立计算水泥厚度和原始地层密度。为了证明该方法的可行性,在一组标定井中测量了实验数据。数据处理结果表明,所提出的分析模型与实验数据吻合较好。此外,根据实验数据对该方法的应用条件进行了讨论。最后,几个测井实例表明,在合适的套管和水泥条件下,用该方法计算的套管井地层密度与相应的裸眼密度吻合较好。
{"title":"A NOVEL METHOD FOR FORMATION DENSITY MEASUREMENT IN CASED WELLS","authors":"Xinguang Wang, Dong Li, Lei Zhang, Feng Zhang, Wen-hao Wang, Huiping Gao, China Oilfield Services Limited, Xi 'an HuiNeng Electronic Equipment Co. Ltd","doi":"10.30632/spwla-2021-0065","DOIUrl":"https://doi.org/10.30632/spwla-2021-0065","url":null,"abstract":"Density logging is a fundamental logging method of open-hole logging series. In cased wells, the formation is usually treated as a cylindrical multiple-layer model including borehole, casing pipe, cement, and original formation, producing more complicated logging response compared to open-hole wells. Several studies have indicated that it is feasible to measure formation density using a density logging tool in cased well under certain conditions, of which the key problem is to determine casing thickness, cement thickness and cement density. This study proposes a novel method to evaluate formation density behind casing using a three-detector density tool. Considering that the backscattered gamma is easily affected by the density behind the casing, especially when the casing thickness is small, a designed ratio of high and low energy window counts of the scattered gamma spectrum is used to calculate casing thickness, which improves the lower limit of the detection. After obtaining casing thickness, the cylindrical multiple-layer model of cased well can be simplified to a cylindrical three-layer model, consisting of cased borehole, cement, and the original formation. Subsequently, the analytical expressions of the thickness of the middle layer (cement) as well as the density of the outer layer (original formation) can be derived based on the logging geometric factor concept, G, which borrowing from electric logging. Consequently, the cement thickness and original formation density could be calculated independently. To demonstrate the feasibility of the proposed method, experimental data are measured in a group of calibration wells. The data processing results show that the proposed analytical model fit well with the experimental data. In addition, the application conditions of the novel method are discussed based on experimental data. Finally, several log examples illustrate that under suitable casing and cement conditions, the formation density calculated with the proposed method in cased borehole is in good agreement with the corresponding open-hole density.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"98 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124785833","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
AN ACCURATELY DETERMINING POROSITY METHOD FROM PULSED-NEUTRON ELEMENT LOGGING IN UNCONVENTIONAL RESERVOIRS 非常规储层脉冲中子元素测井孔隙度精确测定方法
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0115
Feng Zhang, Qiu Fei, Q. Fang, Xiaoyang Zhang, Hui Zhang, F. Tang, Fan Jilin
Unconventional reservoirs have low porosity and complex mineral composition containing quartz, feldspar, calcite, dolomite, pyrite and kerogen, which may seriously reduce the accuracy of the porosity measurement. The multi-detector pulsed neutron logging technique was already used for determining porosity through the combination of inelastic and capture gamma ray information in different spacing. In this paper, the new parameter, which is characterized by thermal neutron count ratio and lithology factor based on element content, is proposed to determine porosity from the three-detector pulsed neutron element logging in unconventional reservoir. To evaluate mineral composition, lithology, and gas/oil/water saturation in unconventional reservoir, a new multi-detector pulsed neutron logging tool was put out. The instrument consists of two He-3 thermal neutron detectors and a LaBr3 gamma detector. Therefore, the combination of thermal neutron count ratio between near detector and long detector with lithology factor of element content can measure neutron porosity and eliminate the influence of complex lithology. Based on some calibration pit data measured in laboratory, as well as the numerical simulation method, the influences of different lithological characters and mineral types on the neutron count ratio were studied. Meanwhile, large numbers of stratigraphic models with different lithological characters and different mineral compositions were established using Monte Carlo simulation method, and the content of silicon, calcium, hydrogen, oxygen, magnesium, aluminum and iron under different stratigraphic conditions was determined by the spectral element solution. A regression analysis was conducted to establish the relationship between the content of elements and the lithologic factor. The count ratio difference stemming from different lithological and mineral compositions was eliminated through a combination of lithological correction factor and thermal neutron count ratio. Different mineral compositions of stratigraphic simulation models were set up for verification. The absolute error of porosity measurement was less than 1.0p.u. in the formations with porosity less than 15p.u., which verified the accuracy of this method for neutron porosity evaluation in complex lithological characters formations. Two field examples were processed by this new parameter which in combination of thermal neutron count ratio and formation elements content information from the three-detector pulsed neutron instrument, which indicated a good accuracy for unconventional oil and gas reservoir evaluation.
非常规储层孔隙度低,矿物组成复杂,含石英、长石、方解石、白云石、黄铁矿、干酪根等,可能严重降低孔隙度测量的准确性。多探测器脉冲中子测井技术已被用于通过结合非弹性和捕获不同间距的伽马射线信息来确定孔隙度。本文提出了以热中子计数比和基于元素含量的岩性因子为特征的新参数来确定非常规储层的孔隙度。为评价非常规储层的矿物组成、岩性及气/油/水饱和度,研制了一种新型多探测器脉冲中子测井工具。该仪器由两台He-3热中子探测器和一台LaBr3伽马探测器组成。因此,将近探测器与长探测器之间的热中子计数比与元素含量的岩性因子相结合,可以测量中子孔隙度,消除复杂岩性的影响。基于实验室标定坑实测数据,结合数值模拟方法,研究了不同岩性特征和矿物类型对中子计数比的影响。同时,利用蒙特卡罗模拟方法建立了大量具有不同岩性特征和不同矿物组成的地层模型,利用谱元溶液测定了不同地层条件下硅、钙、氢、氧、镁、铝、铁的含量。通过回归分析,建立了元素含量与岩性因素的关系。通过结合岩性校正系数和热中子计数比,消除了不同岩性和矿物组成造成的计数比差异。建立了不同矿物组成的地层模拟模型进行验证。孔隙度测量的绝对误差小于1.0p.u。在孔隙度小于15p.u的地层中。验证了该方法在复杂岩性地层中子孔隙度评价中的准确性。该参数结合热中子计数比和三探测器脉冲中子仪的地层元素含量信息对2个油田实例进行了处理,结果表明该参数在非常规油气藏评价中具有较好的准确性。
{"title":"AN ACCURATELY DETERMINING POROSITY METHOD FROM PULSED-NEUTRON ELEMENT LOGGING IN UNCONVENTIONAL RESERVOIRS","authors":"Feng Zhang, Qiu Fei, Q. Fang, Xiaoyang Zhang, Hui Zhang, F. Tang, Fan Jilin","doi":"10.30632/SPWLA-2021-0115","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0115","url":null,"abstract":"Unconventional reservoirs have low porosity and complex mineral composition containing quartz, feldspar, calcite, dolomite, pyrite and kerogen, which may seriously reduce the accuracy of the porosity measurement. The multi-detector pulsed neutron logging technique was already used for determining porosity through the combination of inelastic and capture gamma ray information in different spacing. In this paper, the new parameter, which is characterized by thermal neutron count ratio and lithology factor based on element content, is proposed to determine porosity from the three-detector pulsed neutron element logging in unconventional reservoir. To evaluate mineral composition, lithology, and gas/oil/water saturation in unconventional reservoir, a new multi-detector pulsed neutron logging tool was put out. The instrument consists of two He-3 thermal neutron detectors and a LaBr3 gamma detector. Therefore, the combination of thermal neutron count ratio between near detector and long detector with lithology factor of element content can measure neutron porosity and eliminate the influence of complex lithology. Based on some calibration pit data measured in laboratory, as well as the numerical simulation method, the influences of different lithological characters and mineral types on the neutron count ratio were studied. Meanwhile, large numbers of stratigraphic models with different lithological characters and different mineral compositions were established using Monte Carlo simulation method, and the content of silicon, calcium, hydrogen, oxygen, magnesium, aluminum and iron under different stratigraphic conditions was determined by the spectral element solution. A regression analysis was conducted to establish the relationship between the content of elements and the lithologic factor. The count ratio difference stemming from different lithological and mineral compositions was eliminated through a combination of lithological correction factor and thermal neutron count ratio. Different mineral compositions of stratigraphic simulation models were set up for verification. The absolute error of porosity measurement was less than 1.0p.u. in the formations with porosity less than 15p.u., which verified the accuracy of this method for neutron porosity evaluation in complex lithological characters formations. Two field examples were processed by this new parameter which in combination of thermal neutron count ratio and formation elements content information from the three-detector pulsed neutron instrument, which indicated a good accuracy for unconventional oil and gas reservoir evaluation.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"4 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130668366","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
INVERSION-BASED MEASUREMENT INTERPRETATION OF A NEW THROUGH-THE-BIT OBM PHOTOREALISTIC BOREHOLE IMAGER 一种新型过钻头obm逼真井眼成像仪基于反演的测量解释
Pub Date : 2021-05-17 DOI: 10.30632/spwla-2021-0080
Yong-hua Chen, Schlumberger, R. Bloemenkamp, Peter Schlicht, Lin Liang, L. Comparon
A new 2 1/8-in. outer-diameter photorealistic imager for oil-based muds (OBM) has recently started field testing in unconventional formations in North America. To obtain the best interpretation of its measurements, a two-step quantitative inversion workflow has been developed with a performance similar to the existing inversion workflows for the regular high-definition OBM imagers. The new inversion workflow provides borehole resistivity images, borehole rugosity images, and borehole dielectric permittivity images as well as multiple quality curves. The modeling of the new borehole imager is performed with a 2D axisymmetric finite element code. An efficient forward model is developed by fitting the tool response tables into fourth-order polynomials in terms of the sensor standoff, formation, and mud impedivities for broad ranges of model parameters. The fast forward model based on the polynomial fitting is calibrated against the actual tool measurements in a laboratory setup and applied in the inversion algorithms. The inversion workflow is tested with synthetic data and the inverted model parameters are compared with their true values to study and analyze their corresponding measurement sensitivity and optimize the inversion input parameters. It is used to invert several field test datasets in unconventional wells. The results show that the inversion results provide critical added value for formation evaluation, showing geological features that would otherwise be missed, such as fracture properties. Projection-based formation impedivity images, as available for the regular high-definition OBM imagers, are ideal for conductive formations but suffer from a rollover effect in resistive formations. In comparison, the image formed from the inverted formation resistivity does not roll over and is more consistent for resistive formations. The image formed by the inverted standoff reflects surface conditions of the borehole and can be used to interpret whether the fractures and the faults are open, closed, or damaged in the drilling process. Multiple image examples are given from unconventional wells to demonstrate that the inverted standoff image can reveal fractures when there is insufficient or even no contrast in medium properties. The inverted standoff image also serves as a diagnostic tool for interpreting borehole and tool conditions during the measurements. The inverted permittivity may have a larger dynamic range than the resistivity especially for unconventional formations, thus providing an alternative and potentially clearer borehole image.
新的2 1/8英寸。油基泥浆(OBM)外径逼真成像仪最近开始在北美非常规地层进行现场测试。为了获得其测量结果的最佳解释,开发了一种两步定量反演工作流程,其性能与常规高清OBM成像仪的现有反演工作流程相似。新的反演工作流程可提供井眼电阻率图像、井眼粗糙度图像、井眼介电常数图像以及多条质量曲线。采用二维轴对称有限元程序对新型井眼成像仪进行建模。在广泛的模型参数范围内,通过将工具响应表拟合成基于传感器距离、地层和泥浆阻力的四阶多项式,建立了一个有效的正演模型。基于多项式拟合的快进模型在实验室设置中针对实际工具测量进行校准,并应用于反演算法。利用合成数据对反演工作流程进行测试,并将反演模型参数与其真实值进行比较,研究分析其对应的测量灵敏度,优化反演输入参数。它被用于反演非常规井的多个现场测试数据集。结果表明,反演结果为储层评价提供了重要的附加价值,显示了裂缝性质等可能被忽略的地质特征。常规高清晰度OBM成像仪可获得基于投影的地层阻抗图像,这对于导电地层来说是理想的,但在电阻地层中会受到翻转效应的影响。相比之下,由地层电阻率反演得到的图像不会翻转,并且对于电阻性地层来说更加一致。倒立成像反映了井眼的表面情况,可以用来解释裂缝和断层在钻井过程中是张开的、闭合的还是损坏的。以非常规井为例,说明在介质物性对比不足甚至没有对比的情况下,倒立成像也能显示裂缝。在测量过程中,倒立图像还可以作为解释井眼和工具状况的诊断工具。反向介电常数可能比电阻率具有更大的动态范围,特别是对于非常规地层,因此提供了一种替代的、可能更清晰的井眼图像。
{"title":"INVERSION-BASED MEASUREMENT INTERPRETATION OF A NEW THROUGH-THE-BIT OBM PHOTOREALISTIC BOREHOLE IMAGER","authors":"Yong-hua Chen, Schlumberger, R. Bloemenkamp, Peter Schlicht, Lin Liang, L. Comparon","doi":"10.30632/spwla-2021-0080","DOIUrl":"https://doi.org/10.30632/spwla-2021-0080","url":null,"abstract":"A new 2 1/8-in. outer-diameter photorealistic imager for oil-based muds (OBM) has recently started field testing in unconventional formations in North America. To obtain the best interpretation of its measurements, a two-step quantitative inversion workflow has been developed with a performance similar to the existing inversion workflows for the regular high-definition OBM imagers. The new inversion workflow provides borehole resistivity images, borehole rugosity images, and borehole dielectric permittivity images as well as multiple quality curves. The modeling of the new borehole imager is performed with a 2D axisymmetric finite element code. An efficient forward model is developed by fitting the tool response tables into fourth-order polynomials in terms of the sensor standoff, formation, and mud impedivities for broad ranges of model parameters. The fast forward model based on the polynomial fitting is calibrated against the actual tool measurements in a laboratory setup and applied in the inversion algorithms. The inversion workflow is tested with synthetic data and the inverted model parameters are compared with their true values to study and analyze their corresponding measurement sensitivity and optimize the inversion input parameters. It is used to invert several field test datasets in unconventional wells. The results show that the inversion results provide critical added value for formation evaluation, showing geological features that would otherwise be missed, such as fracture properties. Projection-based formation impedivity images, as available for the regular high-definition OBM imagers, are ideal for conductive formations but suffer from a rollover effect in resistive formations. In comparison, the image formed from the inverted formation resistivity does not roll over and is more consistent for resistive formations. The image formed by the inverted standoff reflects surface conditions of the borehole and can be used to interpret whether the fractures and the faults are open, closed, or damaged in the drilling process. Multiple image examples are given from unconventional wells to demonstrate that the inverted standoff image can reveal fractures when there is insufficient or even no contrast in medium properties. The inverted standoff image also serves as a diagnostic tool for interpreting borehole and tool conditions during the measurements. The inverted permittivity may have a larger dynamic range than the resistivity especially for unconventional formations, thus providing an alternative and potentially clearer borehole image.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"27 1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130927918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A NEW LOOK AT THE DUAL DEPTH OF INVESTIGATION OF LWD PROPAGATION RESISTIVITY LOGGING 随钻传播电阻率测井双深度研究新进展
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0119
G. Wang, D. Homan, D. Maggs, David G. Allen
It is well established that phase shift and attenuation measurements acquired by an electromagnetic propagation tool come with different depths of investigation (DOI). The attenuation measurement sees deeper into the formation than the phase shift measurement. This difference has been reported not only for the 2 MHz propagation resistivity tool, but also for the deep propagation tool that operates at 25 MHz. Although the difference has been demonstrated with modeling, test tank experiments and logs, a complete physical explanation has been notably absent since the introduction of the MHz-frequency propagation logging in 1980s. The question is so intriguing that it has been raised repeatedly over the past decades: what drives the difference of DOI for the two measurements that are acquired with the same electromagnetic field? In this paper, we revisit this problem with an aim of providing a physical insight to bridge the gap between theory and application. This is an extension of our recent work on the theory of apparent conductivity for propagation measurements. We address the problem by applying high-order geometric theory for low-frequency electromagnetic problems in lossy media in conjunction with the Taylor series expansion for the voltage ratio measured by a propagation tool. In so doing, we find that in a resistive formation where the dielectric effect is small: 1) the phase shift measurement is primarily due to the first-order eddy current induced in the formation; 2) in contrast, the leading source of the attenuation measurement is the second-order eddy current. Since the second-order eddy current is more spread out than the first-order eddy current, this explains why the DOI of attenuation resistivity is larger than that of phase shift resistivity. The difference in spatial distribution of two eddy currents is also the reason for the difference of vertical resolution between the two. The same root cause for the difference of DOI and vertical resolution also holds when comparing R-signal and X-signal from induction resistivity logging. Other properties shared by propagation and induction resistivity logging will be discussed, such as skin effect and dielectric effect, as well as their asymptotic properties in high-resistivity formations. We conclude that propagation and induction resistivity logging are essentially similar, even though the two measurement principles may seem rather different.
众所周知,通过电磁传播工具获得的相移和衰减测量值具有不同的调查深度(DOI)。衰减测量比相移测量看到更深的地层。据报道,这种差异不仅适用于2 MHz传播电阻率工具,也适用于25 MHz的深度传播工具。尽管这种差异已经通过建模、测试罐实验和测井证明,但自20世纪80年代引入mhz频率传播测井以来,一直没有一个完整的物理解释。这个问题非常有趣,在过去的几十年里,它被反复提出:是什么导致了在相同电磁场下获得的两个测量值的DOI差异?在本文中,我们重新审视这个问题,目的是提供一个物理洞察力,以弥合理论和应用之间的差距。这是我们最近关于传播测量的视电导率理论工作的延伸。我们将高阶几何理论应用于有耗介质中低频电磁问题,并结合传播工具测量电压比的泰勒级数展开来解决这个问题。在这样做的过程中,我们发现在介电效应较小的电阻地层中:1)相移测量主要是由于地层中产生的一阶涡流;2)相比之下,衰减测量的主要来源是二阶涡流。由于二阶涡流比一阶涡流更分散,这就解释了衰减电阻率的DOI比相移电阻率的DOI大。两种涡流在空间分布上的差异也是导致两者垂直分辨率差异的原因。在比较感应电阻率测井的r -信号和x -信号时,造成DOI和垂直分辨率差异的根本原因同样成立。将讨论传播电阻率和感应电阻率测井共有的其他性质,如集肤效应和介电效应,以及它们在高电阻率地层中的渐近性质。我们得出结论,传播电阻率测井和感应电阻率测井在本质上是相似的,尽管这两种测量原理可能看起来相当不同。
{"title":"A NEW LOOK AT THE DUAL DEPTH OF INVESTIGATION OF LWD PROPAGATION RESISTIVITY LOGGING","authors":"G. Wang, D. Homan, D. Maggs, David G. Allen","doi":"10.30632/SPWLA-2021-0119","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0119","url":null,"abstract":"It is well established that phase shift and attenuation measurements acquired by an electromagnetic propagation tool come with different depths of investigation (DOI). The attenuation measurement sees deeper into the formation than the phase shift measurement. This difference has been reported not only for the 2 MHz propagation resistivity tool, but also for the deep propagation tool that operates at 25 MHz. Although the difference has been demonstrated with modeling, test tank experiments and logs, a complete physical explanation has been notably absent since the introduction of the MHz-frequency propagation logging in 1980s. The question is so intriguing that it has been raised repeatedly over the past decades: what drives the difference of DOI for the two measurements that are acquired with the same electromagnetic field? In this paper, we revisit this problem with an aim of providing a physical insight to bridge the gap between theory and application. This is an extension of our recent work on the theory of apparent conductivity for propagation measurements. We address the problem by applying high-order geometric theory for low-frequency electromagnetic problems in lossy media in conjunction with the Taylor series expansion for the voltage ratio measured by a propagation tool. In so doing, we find that in a resistive formation where the dielectric effect is small: 1) the phase shift measurement is primarily due to the first-order eddy current induced in the formation; 2) in contrast, the leading source of the attenuation measurement is the second-order eddy current. Since the second-order eddy current is more spread out than the first-order eddy current, this explains why the DOI of attenuation resistivity is larger than that of phase shift resistivity. The difference in spatial distribution of two eddy currents is also the reason for the difference of vertical resolution between the two. The same root cause for the difference of DOI and vertical resolution also holds when comparing R-signal and X-signal from induction resistivity logging. Other properties shared by propagation and induction resistivity logging will be discussed, such as skin effect and dielectric effect, as well as their asymptotic properties in high-resistivity formations. We conclude that propagation and induction resistivity logging are essentially similar, even though the two measurement principles may seem rather different.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"115361996","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
ENHANCED MINERAL QUANTIFICATION AND UNCERTAINTY ANALYSIS FROM DOWNHOLE SPECTROSCOPY LOGS USING VARIATIONAL AUTOENCODERS 利用变分自动编码器增强了井下光谱测井的矿物定量和不确定度分析
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0069
P. Craddock, Prakhar Srivastava, H. Datir, D. Rose, T. Zhou, L. Mosse, Lalitha Venkataramanan
This paper describes an innovative machine learning application, based on variational autoencoder frameworks, to quantify the concentrations and associated uncertainties of common minerals in sedimentary formations using the measurement of atomic element concentrations from geochemical spectroscopy logs as inputs. The algorithm comprises an input(s), encoder, decoder, output(s), and a novel cost function to optimize the model coefficients during training. The input to the algorithm is a set of dry-weight concentrations of atomic elements with their associated uncertainty. The first output is a set of dry-weight fractions of fourteen minerals, and the second output is a set of reconstructed dry-weight concentrations of the original elements. Both sets of outputs include estimates of uncertainty on their predictions. The encoder and decoder are multilayer feed-forward artificial neural networks (ANN), with their coefficients (weights) optimized during calibration (training). The cost function simultaneously minimizes error (the accuracy metric) and variance (the precision or robustness metric) on the mineral and reconstructed elemental outputs. Training of the weights is done using a set of several-thousand core samples with independent, high-fidelity elemental and mineral (quartz, potassium-feldspar, plagioclase-feldspar, illite, smectite, kaolinite, chlorite, mica, calcite, dolomite, ankerite, siderite, pyrite, and anhydrite) data. The algorithm provides notable advantages over existing methods to estimate formation lithology or mineralogy relying on simple linear, empirical, or nearest-neighbor functions. The ANN numerically capture the multi-dimensional and nonlinear geochemical relationship (mapping) between elements and minerals that is insufficiently described by prior methods. Training is iterative via backpropagation and samples from Gaussian distributions on each of the elemental inputs, rather than single values, for every sample at each iteration (epoch). These Gaussian distributions are chosen to specifically represent the unique statistical uncertainty of the dry-weight elements in the logging measurements. Sampling from Gaussian distributions during training reduces the potential for overfitting, provides robustness for log interpretations, and further enables a calibrated estimate of uncertainty on the mineral and reconstructed elemental outputs, all of which are lacking in prior methods. The framework of the algorithm is purposefully generalizable that it can be adapted across geochemical spectroscopy tools. The algorithm reasonably approximates a ‘global-average’ model that requires neither different calibrations nor expert parameterization or intervention for interpreting common oilfield sedimentary formations, although the framework is again purposefully generalizable so it can be optimized for local environments where desirable. The paper showcases field application of the method for estimating mineral type and abundance in oil
本文描述了一种基于变分自编码器框架的创新机器学习应用程序,该应用程序使用地球化学光谱测井中原子元素浓度的测量作为输入,量化沉积地层中常见矿物的浓度和相关不确定性。该算法包括输入、编码器、解码器、输出和一个新的成本函数,用于在训练过程中优化模型系数。该算法的输入是一组原子元素的干重浓度及其相关的不确定度。第一个输出是一组14种矿物的干重分数,第二个输出是一组原始元素的重建干重浓度。两组产出都包括对其预测的不确定性的估计。编码器和解码器是多层前馈人工神经网络(ANN),其系数(权值)在标定(训练)过程中优化。成本函数同时最小化矿物和重建元素输出的误差(精度度量)和方差(精度或鲁棒性度量)。权重的训练是使用一组数千个独立的、高保真的元素和矿物(石英、钾长石、斜长石、伊利石、蒙脱石、高岭石、绿泥石、云母、方解石、白云石、铁云石、菱铁矿、黄铁矿和硬石膏)数据的岩心样本完成的。与现有的依赖简单线性、经验或最近邻函数估计地层岩性或矿物学的方法相比,该算法具有显著的优势。人工神经网络在数值上捕捉了元素和矿物之间的多维和非线性地球化学关系(映射),这是以前的方法无法充分描述的。训练是通过反向传播和对每个元素输入的高斯分布样本进行迭代的,而不是每次迭代(epoch)的每个样本的单个值。选择这些高斯分布是为了具体表示测井测量中干重元素的独特统计不确定性。在训练过程中从高斯分布中采样减少了过拟合的可能性,为对数解释提供了鲁棒性,并进一步实现了对矿物和重建元素输出的不确定性的校准估计,所有这些都是以前的方法所缺乏的。该算法的框架是有目的的推广,它可以适应各种地球化学光谱工具。该算法合理地近似于“全球平均”模型,既不需要不同的校准,也不需要专家参数化或干预来解释常见的油田沉积地层,尽管该框架也是有目的地一般化的,因此可以针对需要的局部环境进行优化。本文介绍了利用测井资料估算油田地层矿物类型和丰度的方法在现场的应用。
{"title":"ENHANCED MINERAL QUANTIFICATION AND UNCERTAINTY ANALYSIS FROM DOWNHOLE SPECTROSCOPY LOGS USING VARIATIONAL AUTOENCODERS","authors":"P. Craddock, Prakhar Srivastava, H. Datir, D. Rose, T. Zhou, L. Mosse, Lalitha Venkataramanan","doi":"10.30632/SPWLA-2021-0069","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0069","url":null,"abstract":"This paper describes an innovative machine learning application, based on variational autoencoder frameworks, to quantify the concentrations and associated uncertainties of common minerals in sedimentary formations using the measurement of atomic element concentrations from geochemical spectroscopy logs as inputs. The algorithm comprises an input(s), encoder, decoder, output(s), and a novel cost function to optimize the model coefficients during training. The input to the algorithm is a set of dry-weight concentrations of atomic elements with their associated uncertainty. The first output is a set of dry-weight fractions of fourteen minerals, and the second output is a set of reconstructed dry-weight concentrations of the original elements. Both sets of outputs include estimates of uncertainty on their predictions. The encoder and decoder are multilayer feed-forward artificial neural networks (ANN), with their coefficients (weights) optimized during calibration (training). The cost function simultaneously minimizes error (the accuracy metric) and variance (the precision or robustness metric) on the mineral and reconstructed elemental outputs. Training of the weights is done using a set of several-thousand core samples with independent, high-fidelity elemental and mineral (quartz, potassium-feldspar, plagioclase-feldspar, illite, smectite, kaolinite, chlorite, mica, calcite, dolomite, ankerite, siderite, pyrite, and anhydrite) data. The algorithm provides notable advantages over existing methods to estimate formation lithology or mineralogy relying on simple linear, empirical, or nearest-neighbor functions. The ANN numerically capture the multi-dimensional and nonlinear geochemical relationship (mapping) between elements and minerals that is insufficiently described by prior methods. Training is iterative via backpropagation and samples from Gaussian distributions on each of the elemental inputs, rather than single values, for every sample at each iteration (epoch). These Gaussian distributions are chosen to specifically represent the unique statistical uncertainty of the dry-weight elements in the logging measurements. Sampling from Gaussian distributions during training reduces the potential for overfitting, provides robustness for log interpretations, and further enables a calibrated estimate of uncertainty on the mineral and reconstructed elemental outputs, all of which are lacking in prior methods. The framework of the algorithm is purposefully generalizable that it can be adapted across geochemical spectroscopy tools. The algorithm reasonably approximates a ‘global-average’ model that requires neither different calibrations nor expert parameterization or intervention for interpreting common oilfield sedimentary formations, although the framework is again purposefully generalizable so it can be optimized for local environments where desirable. The paper showcases field application of the method for estimating mineral type and abundance in oil","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"20 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125972028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
ULTRASONIC ANGLE-DEPENDENT REFLECTIVITY IN COMPLEX ROCKS FOR IMPROVED INTERPRETATION OF SONIC AND ULTRASONIC LOGS 复杂岩石的超声角度相关反射率,以改善声波和超声测井的解释
Pub Date : 2021-05-17 DOI: 10.30632/spwla-2021-0058
Daria Olszowska, Gabriel Gallardo-Giozza, C. Torres‐Verdín
Porous rocks are rarely homogeneous. Significant spatial variations in elastic properties are often observed in rocks due to depositional, diagenetic, and structural processes. In laminated sandstones, complex carbonates, or unconventional formations, elastic properties can vary on scales from millimeters to tens of meters. Detection of inhomogeneities and their size in rocks is crucial for fracture propagation design, height containment assessment, and for improving well/reservoir productivity. Most laboratory techniques used to measure rock elastic properties fail to distinguish mid-scale anisotropy; results are subject to spatial averaging effects. We introduce a new experimental method to measure continuous compressional- and shear-wave logs of core samples based on measurements of angle-dependent ultrasonic reflection coefficients. Simultaneously with reflected waves, we detect and interpret refracted waves as an independent way to estimate acoustic wave velocities to support the analysis. Our laboratory system is equipped with an array of receivers to continuously collect measurements. At each core location, we acquire acoustic waveforms at multiple transmitter-receiver angles using a pitch-catch acquisition mode (similar to standard sonic tools). This acquisition mode uses multiple receivers, allowing us to obtain measurements at different incidence angles without moving the sample and keeping the distance traveled by reflected waves constant, thereby eliminating the need for geometrical spreading corrections in reflection-coefficient calculations. Reflectivity-vs.-angle measurements are then matched with numerical simulations to estimate rock elastic properties. Ultrasonic reflection-coefficient measurements are successfully used to estimate dynamic elastic rock properties of homogeneous and layered rock samples. For homogenous samples, values are within a 5% range when compared to those obtained with the standard acoustic transmission method. Measurements acquired on natural and artificially constructed samples show significant departures from homogeneous behavior caused by layering. Laboratory reflection-coefficient measurements enable detection of inch-scale anisotropy within the rock, leading to improved assessment of formation elastic properties. Furthermore, continuous core measurements provide high-resolution reflection-coefficient information which is complementary to open-hole ultrasonic logs.
多孔岩石很少是均匀的。由于沉积作用、成岩作用和构造作用,在岩石中经常观察到弹性性质的显著空间变化。在层状砂岩、复杂碳酸盐或非常规地层中,弹性特性可以在毫米到几十米的尺度上变化。检测岩石中的非均质性及其尺寸对于裂缝扩展设计、高度控制评估以及提高油井/油藏产能至关重要。大多数用于测量岩石弹性特性的实验室技术无法区分中等尺度的各向异性;结果受空间平均效应影响。本文介绍了一种测量岩心样品连续纵波和横波测井曲线的实验方法,该方法基于测量角度相关的超声反射系数。与反射波同时,我们检测和解释折射波作为一种独立的方法来估计声波速度,以支持分析。我们的实验室系统配备了一系列接收器来连续收集测量数据。在每个核心位置,我们使用音调捕获采集模式(类似于标准声波工具)获取多个发射器-接收器角度的声波波形。这种采集模式使用多个接收器,使我们能够在不移动样品的情况下获得不同入射角的测量结果,并保持反射波传播的距离恒定,从而消除了在反射系数计算中对几何扩展修正的需要。Reflectivity-vs。然后将-角测量值与数值模拟相匹配,以估计岩石的弹性特性。超声反射系数测量成功地用于估计均匀和层状岩石样品的动态弹性岩石性质。对于均匀样品,与标准声透射法获得的值相比,值在5%的范围内。在自然和人工构造的样品上获得的测量结果表明,分层造成的均匀行为明显偏离。实验室反射系数测量可以检测岩石内部英寸尺度的各向异性,从而改进对地层弹性特性的评估。此外,连续岩心测量可以提供高分辨率的反射系数信息,这是裸眼超声测井的补充。
{"title":"ULTRASONIC ANGLE-DEPENDENT REFLECTIVITY IN COMPLEX ROCKS FOR IMPROVED INTERPRETATION OF SONIC AND ULTRASONIC LOGS","authors":"Daria Olszowska, Gabriel Gallardo-Giozza, C. Torres‐Verdín","doi":"10.30632/spwla-2021-0058","DOIUrl":"https://doi.org/10.30632/spwla-2021-0058","url":null,"abstract":"Porous rocks are rarely homogeneous. Significant spatial variations in elastic properties are often observed in rocks due to depositional, diagenetic, and structural processes. In laminated sandstones, complex carbonates, or unconventional formations, elastic properties can vary on scales from millimeters to tens of meters. Detection of inhomogeneities and their size in rocks is crucial for fracture propagation design, height containment assessment, and for improving well/reservoir productivity. Most laboratory techniques used to measure rock elastic properties fail to distinguish mid-scale anisotropy; results are subject to spatial averaging effects. We introduce a new experimental method to measure continuous compressional- and shear-wave logs of core samples based on measurements of angle-dependent ultrasonic reflection coefficients. Simultaneously with reflected waves, we detect and interpret refracted waves as an independent way to estimate acoustic wave velocities to support the analysis. Our laboratory system is equipped with an array of receivers to continuously collect measurements. At each core location, we acquire acoustic waveforms at multiple transmitter-receiver angles using a pitch-catch acquisition mode (similar to standard sonic tools). This acquisition mode uses multiple receivers, allowing us to obtain measurements at different incidence angles without moving the sample and keeping the distance traveled by reflected waves constant, thereby eliminating the need for geometrical spreading corrections in reflection-coefficient calculations. Reflectivity-vs.-angle measurements are then matched with numerical simulations to estimate rock elastic properties. Ultrasonic reflection-coefficient measurements are successfully used to estimate dynamic elastic rock properties of homogeneous and layered rock samples. For homogenous samples, values are within a 5% range when compared to those obtained with the standard acoustic transmission method. Measurements acquired on natural and artificially constructed samples show significant departures from homogeneous behavior caused by layering. Laboratory reflection-coefficient measurements enable detection of inch-scale anisotropy within the rock, leading to improved assessment of formation elastic properties. Furthermore, continuous core measurements provide high-resolution reflection-coefficient information which is complementary to open-hole ultrasonic logs.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"105 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132223512","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
RAPID CROSS-PLOT DISCRIMINATION OF COMMERCIAL POTASH MINERALIZATION – CASE HISTORIES 商业钾肥矿化的快速交叉图判别——案例历史
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0109
D. Hill, E. Crain
Potash minerals are a source of potassium, which is used for the manufacture of gunpowder and fertilizer. Commercial potash mineralization is often discovered when petroleum wells are drilled through evaporite sequences and the Gamma Ray log “goes off scale”. This is because potassium is one of the naturally occurring radioactive elements, emitting gamma rays from the 40K isotope, in its decay to 40Ar. However, not all potash minerals may be commercial sources of potassium via underground mechanical or solution mining techniques and Potassium is not the only radioactive element. For example, the mineralogy of the McNutt “Potash” member of the Salado Formation in SE New Mexico, is extremely complex, consisting of multiple thin (i.e., less than 10 ft thick) beds of six low-grade (radioactive) potash minerals, only two of which are commercial. There are also four non-radioactive evaporite minerals, one of which interferes with potash milling chemistry, and numerous claystones and Marker Beds (shales), with GR count rates comparable to the low-grade potash. Because of this complexity, traditional wireline and Logging While Drilling Potash Assay techniques, such as Gamma Ray log-to-core assay transforms, may not be sufficient to identify potentially commercial potash mineralization, for underground mining. Crain and Anderson (1966) and Hill (2019) developed linear programming, and multi-mineral analyses, respectively, to estimate Potash mineralogy and grades. However, both of these approaches require complete sets of multiple log measurements. In SE New Mexico, petroleum wells are drilled through the McNutt “Potash” member of the Salado Formation, with air, cased and drilled out to TD in the underlying sediments, with water based mud. Complete log suites are then run from TD to the casing shoe, with only the GR and neutron logs recorded through the cased evaporite sequence for stratigraphic and structural correlation. As a result, numerous oil and gas wells, in SE New Mexico, have cased hole gamma ray and neutron logs, through the Salado Evaporite. Logs, from these wells could provide a rapid Potash screening database, if used properly. A simple screening cross-plot technique, the Potash Identification (PID) plot, utilizing only Gamma Ray and Neutron Porosity, is proposed and successfully demonstrated, as a potential screening tool. This technique can be used with both open and cased-hole petroleum well logs, as well as core hole wire-line logs, and provides discrimination of commercial potash mineralization from non-commercial (potash and non-potash) radioactive mineralization. Case histories of the use of PID cross plots in the evaporite basins of Michigan, Nova Scotia, Saskatchewan, and SE New Mexico are described. The technique may also be useful in screening potential potash deposits in China, Europe, North Africa, and South America.
钾盐矿物是钾的来源,钾可用于制造火药和化肥。商业钾肥矿化通常是在油井钻穿蒸发岩层,而伽马射线测井曲线“偏离规模”时发现的。这是因为钾是一种天然存在的放射性元素,在其衰变为40Ar时,会从40K同位素释放出伽马射线。然而,并不是所有的钾肥矿物都可以通过地下机械或溶液开采技术成为钾的商业来源,钾也不是唯一的放射性元素。例如,新墨西哥州东南部Salado组的McNutt“Potash”成员的矿物学非常复杂,由6种低品位(放射性)钾盐矿物组成的多个薄层(即厚度小于10英尺)组成,其中只有两种具有商业价值。还有四种非放射性蒸发岩矿物,其中一种会干扰钾肥的研磨化学,还有许多粘土和标志层(页岩),其GR计数率与低品位钾肥相当。由于这种复杂性,传统的电缆和随钻随测钾盐分析技术,如伽马射线测井-岩心分析转换,可能不足以识别潜在的商业钾盐矿化,用于地下开采。Crain和Anderson(1966)以及Hill(2019)分别开发了线性规划和多矿物分析,以估计钾肥矿物学和品位。然而,这两种方法都需要完整的多个日志测量集。在新墨西哥州东南部,通过Salado地层的McNutt“Potash”成员钻取油井,使用空气、套管和水基泥浆,在下部沉积物中钻至深度。然后从TD到套管鞋进行完整的测井,仅通过套管蒸发岩层序记录GR和中子测井,以进行地层和构造对比。因此,新墨西哥州东南部的许多油气井都通过Salado蒸发岩进行了套管井伽马射线和中子测井。如果使用得当,这些井的测井数据可以提供一个快速的碳酸钾筛选数据库。提出了一种简单的筛选交叉图技术,即仅利用伽马射线和中子孔隙度的钾鉴别(PID)图,并成功地证明了它是一种潜在的筛选工具。该技术可用于裸眼和套管井测井,以及岩心井电缆测井,并可区分商业钾盐矿化与非商业(钾盐和非钾盐)放射性矿化。描述了在密歇根州、新斯科舍省、萨斯喀彻温省和新墨西哥东南部的蒸发岩盆地使用PID交叉图的案例历史。该技术也可用于筛选中国、欧洲、北非和南美的潜在钾肥矿床。
{"title":"RAPID CROSS-PLOT DISCRIMINATION OF COMMERCIAL POTASH MINERALIZATION – CASE HISTORIES","authors":"D. Hill, E. Crain","doi":"10.30632/SPWLA-2021-0109","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0109","url":null,"abstract":"Potash minerals are a source of potassium, which is used for the manufacture of gunpowder and fertilizer. Commercial potash mineralization is often discovered when petroleum wells are drilled through evaporite sequences and the Gamma Ray log “goes off scale”. This is because potassium is one of the naturally occurring radioactive elements, emitting gamma rays from the 40K isotope, in its decay to 40Ar. However, not all potash minerals may be commercial sources of potassium via underground mechanical or solution mining techniques and Potassium is not the only radioactive element. For example, the mineralogy of the McNutt “Potash” member of the Salado Formation in SE New Mexico, is extremely complex, consisting of multiple thin (i.e., less than 10 ft thick) beds of six low-grade (radioactive) potash minerals, only two of which are commercial. There are also four non-radioactive evaporite minerals, one of which interferes with potash milling chemistry, and numerous claystones and Marker Beds (shales), with GR count rates comparable to the low-grade potash. Because of this complexity, traditional wireline and Logging While Drilling Potash Assay techniques, such as Gamma Ray log-to-core assay transforms, may not be sufficient to identify potentially commercial potash mineralization, for underground mining. Crain and Anderson (1966) and Hill (2019) developed linear programming, and multi-mineral analyses, respectively, to estimate Potash mineralogy and grades. However, both of these approaches require complete sets of multiple log measurements. In SE New Mexico, petroleum wells are drilled through the McNutt “Potash” member of the Salado Formation, with air, cased and drilled out to TD in the underlying sediments, with water based mud. Complete log suites are then run from TD to the casing shoe, with only the GR and neutron logs recorded through the cased evaporite sequence for stratigraphic and structural correlation. As a result, numerous oil and gas wells, in SE New Mexico, have cased hole gamma ray and neutron logs, through the Salado Evaporite. Logs, from these wells could provide a rapid Potash screening database, if used properly. A simple screening cross-plot technique, the Potash Identification (PID) plot, utilizing only Gamma Ray and Neutron Porosity, is proposed and successfully demonstrated, as a potential screening tool. This technique can be used with both open and cased-hole petroleum well logs, as well as core hole wire-line logs, and provides discrimination of commercial potash mineralization from non-commercial (potash and non-potash) radioactive mineralization. Case histories of the use of PID cross plots in the evaporite basins of Michigan, Nova Scotia, Saskatchewan, and SE New Mexico are described. The technique may also be useful in screening potential potash deposits in China, Europe, North Africa, and South America.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130371682","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
USING PROXY SIMULATOR FOR RESERVOIR ZONE SELECTION AND REDUCING THE FORMATION TESTER CLEANUP OPERATIONAL TIME 使用代理模拟器进行储层选择,减少地层测试器清理作业时间
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0114
A. Bertolini, V. Simoes, Marianna Dantas, P. Machado
The filtrate contamination cleanup time on a complex carbonate well using a traditional wireline formation tester (WFT) tool can vary from a couple of hours to over half a day. The method proposed aims at reducing operational time to collect a low-contamination formation fluid sample by determining regions with a smaller depth of invasion using a forward model simulation that considers static and dynamic formation properties to predict the radial profile of invasion. The mud filtrate invasion process was modeled considering the static and dynamic properties of the near-wellbore region in an industry reference reservoir simulator, and it integrates three mechanisms for fluid flow: Darcy’s law, material balance, and capillary pressure. The physical robustness of the reservoir simulator was united to a data-driven model to reduce the computational cost. This proxy model is based on a trained neural network with a broad range of scenarios to predict the numerical simulation results with high accuracy. The invasion estimation from the model is then used to predict the filtrate cleanup time using an industry consolidated numerical modeling. One of the variables influencing most of the cleanup time is the depth of mud filtrate invasion. Thus, reducing this time is a determinant for the WFT operational efficiency. The model for mud invasion has been successfully tested on a complex carbonate well, and the results for the depth of mud invasion were comparable to the results obtained with a commercial data-driven inversion using multiple resistivity channels. The estimated cleanup time using the results of depth of invasion predicted by the forward model has been compared and matched with real carbonate sampling stations, and there was a high correlation indicating that zones with lower depth of invasion required less cleanup time. Besides, using the history-matched cases, different WFT technologies such as single and radial probes, focused, unfocused, and dual-packer WFT inlets were evaluated, showing a high potential for reduction of operational time when properly planned and selected for the specific type of reservoir. The proposed methodology is a viable method for understanding the clean-up behavior in different reservoir scenarios using different WFT technologies. The innovation of this method relies on the data calibration using basic and advanced petrophysical properties through a data-driven model based on a trained neural network to reduce the uncertainty in the predicted invasion radial profile and the WFT cleanup time. The reliability on the theoretical results was increased using real data calibration, and this calibrated theoretical model has been used to guide the sampling depth selection, saving operational time.
在复杂的碳酸盐井中,使用传统的电缆地层测试(WFT)工具进行滤液污染清理的时间从几个小时到半天以上不等。该方法旨在通过考虑静态和动态地层特性的正演模型模拟来预测侵入的径向剖面,从而确定侵入深度较小的区域,从而减少收集低污染地层流体样本的操作时间。在一个行业参考油藏模拟器中,考虑了近井区域的静态和动态特性,对泥浆滤液侵入过程进行了建模,并集成了流体流动的三种机制:达西定律、物质平衡和毛细压力。油藏模拟器的物理鲁棒性与数据驱动模型相结合,降低了计算成本。该代理模型是基于经过训练的具有广泛场景的神经网络,对数值模拟结果进行高精度预测。然后,使用行业统一的数值模拟,将模型中的入侵估计用于预测滤液清理时间。影响大部分清理时间的变量之一是泥浆滤液侵入的深度。因此,减少这一时间是WFT操作效率的决定因素。泥浆侵入模型已在一口复杂的碳酸盐岩井中成功进行了测试,泥浆侵入深度的计算结果与利用多电阻率通道进行的商业数据驱动反演结果相当。利用正演模型预测的侵入深度估算的清理时间与实际碳酸盐采样站进行了比较和匹配,相关性较高,表明侵入深度较低的区域需要较少的清理时间。此外,使用历史匹配的案例,评估了不同的WFT技术,如单探头和径向探头,聚焦、非聚焦和双封隔器WFT入口,显示了在适当规划和选择特定类型油藏时,减少作业时间的巨大潜力。本文提出的方法是一种可行的方法,可用于了解不同WFT技术在不同油藏场景下的清理行为。该方法的创新之处在于,通过基于训练神经网络的数据驱动模型,利用基本和高级岩石物理特性进行数据校准,以减少预测侵入径向剖面的不确定性和WFT清理时间。通过实际数据标定,提高了理论结果的可靠性,并将标定后的理论模型用于指导采样深度的选择,节省了操作时间。
{"title":"USING PROXY SIMULATOR FOR RESERVOIR ZONE SELECTION AND REDUCING THE FORMATION TESTER CLEANUP OPERATIONAL TIME","authors":"A. Bertolini, V. Simoes, Marianna Dantas, P. Machado","doi":"10.30632/SPWLA-2021-0114","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0114","url":null,"abstract":"The filtrate contamination cleanup time on a complex carbonate well using a traditional wireline formation tester (WFT) tool can vary from a couple of hours to over half a day. The method proposed aims at reducing operational time to collect a low-contamination formation fluid sample by determining regions with a smaller depth of invasion using a forward model simulation that considers static and dynamic formation properties to predict the radial profile of invasion. The mud filtrate invasion process was modeled considering the static and dynamic properties of the near-wellbore region in an industry reference reservoir simulator, and it integrates three mechanisms for fluid flow: Darcy’s law, material balance, and capillary pressure. The physical robustness of the reservoir simulator was united to a data-driven model to reduce the computational cost. This proxy model is based on a trained neural network with a broad range of scenarios to predict the numerical simulation results with high accuracy. The invasion estimation from the model is then used to predict the filtrate cleanup time using an industry consolidated numerical modeling. One of the variables influencing most of the cleanup time is the depth of mud filtrate invasion. Thus, reducing this time is a determinant for the WFT operational efficiency. The model for mud invasion has been successfully tested on a complex carbonate well, and the results for the depth of mud invasion were comparable to the results obtained with a commercial data-driven inversion using multiple resistivity channels. The estimated cleanup time using the results of depth of invasion predicted by the forward model has been compared and matched with real carbonate sampling stations, and there was a high correlation indicating that zones with lower depth of invasion required less cleanup time. Besides, using the history-matched cases, different WFT technologies such as single and radial probes, focused, unfocused, and dual-packer WFT inlets were evaluated, showing a high potential for reduction of operational time when properly planned and selected for the specific type of reservoir. The proposed methodology is a viable method for understanding the clean-up behavior in different reservoir scenarios using different WFT technologies. The innovation of this method relies on the data calibration using basic and advanced petrophysical properties through a data-driven model based on a trained neural network to reduce the uncertainty in the predicted invasion radial profile and the WFT cleanup time. The reliability on the theoretical results was increased using real data calibration, and this calibrated theoretical model has been used to guide the sampling depth selection, saving operational time.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"35 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"117265886","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
DEEP DIELECTRIC-BASED WATER SATURATION IN FRESHWATER AND MIXED SALINITY ENVIRONMENTS 淡水和混合盐度环境中深层介电基水饱和度
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0033
Ping Zhang, Wael Abdallah, G. Wang, S. Ma
It is desirable to evaluate the possibility of developing a deeper dielectric permittivity based Sw measurement for various petrophysical applications. The low frequency, (< MHz), resistivity-based method for water saturation (Sw) evaluation is the desired method in the industry due to its deepest depth of investigation (DOI, up to 8 ft). However, the method suffers from higher uncertainty when formation water is very fresh or has mixed salinity. Dielectric permittivity and conductivity dispersion have been used to estimate Sw and salinity. The current dielectric dispersion tools, however, have very shallow DOI due to their high measurement frequency up to GHz, which most likely confines the measurements within the near wellbore mud-filtrate invaded zones. In this study, effective medium-model simulations were conducted to study different electromagnetic (EM) induced-polarization effects and their relationships to rock petrophysical properties. Special attention is placed on the complex conductivity at 2 MHz due to its availability in current logging tools. It is known that the complex dielectric saturation interpretation at the MHz range is quite difficult due to lack of fully understood of physics principles on complex dielectric responses, especially when only single frequency signal is used. Therefore, our study is focused on selected key parameters: water-filled porosity, salinity, and grain shape, and their effects on the modeled formation conductivity and permittivity. To simulate field logs, some of the petrophysical parameters mentioned above are generated randomly within expected ranges. Formation conductivity and permittivity are then calculated using our petrophysical model. The calculated results are then mixed with random noises of 10% to make them more realistic like downhole logs. The synthetic conductivity and permittivity logs are used as inputs in a neural network application to explore possible correlations with water-filled porosity. It is found that while the conductivity and permittivity logs are generated from randomly selected petrophysical parameters, they are highly correlated with water-filled porosity. Furthermore, if new conductivity and permittivity logs are generated with different petrophysical parameters, the correlations defined before can be used to predict water-filled porosity in the new datasets. We also found that for freshwater environments, the conductivity has much lower correlation with water-filled porosity than the one derived from the permittivity. However, the correlations are always improved when both conductivity and permittivity were used. This exercise serves as proof of concept, which opens an opportunity for field data applications. Field logs confirm the findings in the model simulations. Two propagation resistivity logs measured at 2 MHz are processed to calculate formation conductivity and permittivity. Using independently estimated water-filled porosity, a model was trained using a neural
在各种岩石物理应用中,有必要评估开发更深层次介电常数的Sw测量方法的可能性。低频(< MHz)、基于电阻率的水饱和度(Sw)评估方法由于其最深的调查深度(DOI,高达8英尺)而成为业内最理想的方法。然而,当地层水非常新鲜或含盐量混合时,该方法存在较高的不确定性。介质介电常数和电导率色散被用来估计Sw和盐度。然而,目前的介质色散工具由于测量频率高达GHz, DOI非常浅,这很可能将测量限制在近井段的泥浆滤液侵入区域。为了研究不同的电磁诱导极化效应及其与岩石物性的关系,本文进行了有效的介质模型模拟。由于目前的测井工具中有2 MHz的复杂电导率,因此需要特别注意。由于对复杂介电响应的物理原理缺乏充分的了解,特别是当仅使用单频信号时,在MHz范围内的复杂介电饱和解释是相当困难的。因此,我们的研究重点是选择关键参数:充水孔隙度、矿化度和颗粒形状,以及它们对模拟地层电导率和介电常数的影响。为了模拟现场测井,上面提到的一些岩石物理参数是在预期范围内随机生成的。然后使用我们的岩石物理模型计算地层电导率和介电常数。然后将计算结果与10%的随机噪声混合,使其更接近井下测井。合成电导率和介电常数测井被用作神经网络应用的输入,以探索与充水孔隙度的可能相关性。研究发现,虽然电导率和介电常数测井曲线是由随机选择的岩石物理参数生成的,但它们与充水孔隙度高度相关。此外,如果使用不同的岩石物理参数生成新的电导率和介电常数测井曲线,则可以使用之前定义的相关性来预测新数据集中的充水孔隙度。我们还发现,在淡水环境中,电导率与充水孔隙度的相关性远低于介电常数。然而,当同时使用电导率和介电常数时,相关性总是得到改善。该练习可以作为概念验证,为现场数据应用提供了机会。现场测井证实了模型模拟的结果。在2 MHz下测量的两个传播电阻率测井数据被处理以计算地层电导率和介电常数。利用独立估计的含水孔隙度,利用神经网络对其中一条测井曲线进行模型训练。在训练模型中,观察到地层电导率、介电常数和充水孔隙度之间具有良好的相关性。该神经网络生成的模型可用于预测从不同井收集的其他测井资料的含水量,相关系数高达96%。给出了利用电导率和介电常数预测含水孔隙度的最佳实践。其中包括如何有效地训练神经网络相关模型,训练后的模型在不同领域日志中的一般应用。利用已建立的方法,可以在淡水和混合盐度环境中获得基于电介质的深层水饱和度,从而提高地层评价、井位和储层饱和度监测。
{"title":"DEEP DIELECTRIC-BASED WATER SATURATION IN FRESHWATER AND MIXED SALINITY ENVIRONMENTS","authors":"Ping Zhang, Wael Abdallah, G. Wang, S. Ma","doi":"10.30632/SPWLA-2021-0033","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0033","url":null,"abstract":"It is desirable to evaluate the possibility of developing a deeper dielectric permittivity based Sw measurement for various petrophysical applications. The low frequency, (< MHz), resistivity-based method for water saturation (Sw) evaluation is the desired method in the industry due to its deepest depth of investigation (DOI, up to 8 ft). However, the method suffers from higher uncertainty when formation water is very fresh or has mixed salinity. Dielectric permittivity and conductivity dispersion have been used to estimate Sw and salinity. The current dielectric dispersion tools, however, have very shallow DOI due to their high measurement frequency up to GHz, which most likely confines the measurements within the near wellbore mud-filtrate invaded zones. In this study, effective medium-model simulations were conducted to study different electromagnetic (EM) induced-polarization effects and their relationships to rock petrophysical properties. Special attention is placed on the complex conductivity at 2 MHz due to its availability in current logging tools. It is known that the complex dielectric saturation interpretation at the MHz range is quite difficult due to lack of fully understood of physics principles on complex dielectric responses, especially when only single frequency signal is used. Therefore, our study is focused on selected key parameters: water-filled porosity, salinity, and grain shape, and their effects on the modeled formation conductivity and permittivity. To simulate field logs, some of the petrophysical parameters mentioned above are generated randomly within expected ranges. Formation conductivity and permittivity are then calculated using our petrophysical model. The calculated results are then mixed with random noises of 10% to make them more realistic like downhole logs. The synthetic conductivity and permittivity logs are used as inputs in a neural network application to explore possible correlations with water-filled porosity. It is found that while the conductivity and permittivity logs are generated from randomly selected petrophysical parameters, they are highly correlated with water-filled porosity. Furthermore, if new conductivity and permittivity logs are generated with different petrophysical parameters, the correlations defined before can be used to predict water-filled porosity in the new datasets. We also found that for freshwater environments, the conductivity has much lower correlation with water-filled porosity than the one derived from the permittivity. However, the correlations are always improved when both conductivity and permittivity were used. This exercise serves as proof of concept, which opens an opportunity for field data applications. Field logs confirm the findings in the model simulations. Two propagation resistivity logs measured at 2 MHz are processed to calculate formation conductivity and permittivity. Using independently estimated water-filled porosity, a model was trained using a neural ","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"141 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"123759490","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
FORMATION PRESSURE ESTIMATION IN ULTRA-LOW PERMEABLE RESERVOIRS EMPLOYING FORMATION RATE ANALYSIS (FRA) AND ARTIFICIAL INTELLIGENCE CONTROLLED TOOLS 利用地层速率分析(fra)和人工智能控制工具估算超低渗透储层地层压力
Pub Date : 2021-05-17 DOI: 10.30632/SPWLA-2021-0086
Y. Askoul, G. Sibbald, A. Hooker, J. Banks, Total E P Uk Ltd
The necessity of knowing formation pressure is crucial to classifying pressure regimes for better understanding in well planning and to de-risk potential abnormal pressure conditions before any future field development wells are drilled, consequently minimizing operational cost. Wireline formation pressure testing has been a useful and reliable technology, that has evolved to confront the challenge of ultra-low permeable reservoir conditions by innovating and improving pump capability, accuracy in pressure measurements, automated control and the implantation of Formation Rate Analysis (FRA) intertwined with an Artificial Intelligent tool. In any pressure testing, the key factor is to be able to withdraw volume from the formation to generate a disturbance on formation pore pressure that a pressure gauge can measure. In the past this has been a difficult task in ultra-low permeable zones. The new generation of wireline tools are capable of withdrawing volume from ultra-low permeable reservoirs, with mobilities lower than 0.01mD/cP. This has been made possible by utilizing control of the pump speed down to 0.0003cc/s which then gives the operator the ability to test ultra-tight formations without the need for inflatable packers. By pulling down the pressure at an extremely low rate and using Artificial Intelligence to control the rate by knowing the formation rate, a proportional amount of volume can be extracted without calling it a tight test. During the operation by observing the rate, and making sure the pump is not oscillating, which indicates the formation rate is lower than the lowest rate the pump can withdraw, the test can be validated for formation flow and the pressure transient of the build – up can be analysed to confirm that at least spherical flow is observed. Once reservoir communication has been confirmed and by analysing drawdown and build-up pressure versus volume withdrawn and implementing the FRA equation, the reservoir pressure can be back calculated by considering isothermal compressibility and FRA slope. This paper highlights the best technical approach to quality check and quality control these tests, showing examples of various wells, where the technique has been used to predict a formation pressure, which can be used for further use for field development, drilling optimisation and production profiles. These pressures would never have been possible using static rates and volume.
了解地层压力对于分类压力状态至关重要,有助于更好地了解井计划,并在未来的任何油田开发井钻探之前消除潜在异常压力状况的风险,从而最大限度地降低运营成本。电缆地层压力测试是一项有用且可靠的技术,通过创新和提高泵的性能、压力测量的准确性、自动化控制以及地层速率分析(FRA)的植入与人工智能工具的结合,电缆地层压力测试已经发展成为一项应对超低渗透油藏条件挑战的技术。在任何压力测试中,关键因素是能够从地层中提取体积,从而对压力表可以测量的地层孔隙压力产生扰动。在过去,这在超低渗透层是一项艰巨的任务。新一代电缆工具能够从超低渗透油藏中提取体积,移动速度低于0.01mD/cP。通过将泵速控制到0.0003cc/s,这使得作业人员能够在不需要膨胀封隔器的情况下测试超密地层。通过以极低的速率降低压力,并利用人工智能通过了解地层速率来控制速率,可以在不称之为紧密测试的情况下提取出一定比例的体积。在作业过程中,通过观察速率,并确保泵没有振荡,这表明地层速率低于泵可以提取的最低速率,可以验证地层流动的测试,并且可以分析堆积的压力瞬态,以确认至少观察到球形流动。一旦确认了储层的连通,通过分析降压和积聚压力与抽采体积的关系,并执行FRA方程,就可以通过考虑等温压缩率和FRA斜率来计算储层压力。本文重点介绍了质量检查和质量控制这些测试的最佳技术方法,并展示了各种井的示例,其中该技术已用于预测地层压力,可用于进一步的油田开发,钻井优化和生产剖面。如果使用静态速率和体积,这些压力是不可能的。
{"title":"FORMATION PRESSURE ESTIMATION IN ULTRA-LOW PERMEABLE RESERVOIRS EMPLOYING FORMATION RATE ANALYSIS (FRA) AND ARTIFICIAL INTELLIGENCE CONTROLLED TOOLS","authors":"Y. Askoul, G. Sibbald, A. Hooker, J. Banks, Total E P Uk Ltd","doi":"10.30632/SPWLA-2021-0086","DOIUrl":"https://doi.org/10.30632/SPWLA-2021-0086","url":null,"abstract":"The necessity of knowing formation pressure is crucial to classifying pressure regimes for better understanding in well planning and to de-risk potential abnormal pressure conditions before any future field development wells are drilled, consequently minimizing operational cost. Wireline formation pressure testing has been a useful and reliable technology, that has evolved to confront the challenge of ultra-low permeable reservoir conditions by innovating and improving pump capability, accuracy in pressure measurements, automated control and the implantation of Formation Rate Analysis (FRA) intertwined with an Artificial Intelligent tool. In any pressure testing, the key factor is to be able to withdraw volume from the formation to generate a disturbance on formation pore pressure that a pressure gauge can measure. In the past this has been a difficult task in ultra-low permeable zones. The new generation of wireline tools are capable of withdrawing volume from ultra-low permeable reservoirs, with mobilities lower than 0.01mD/cP. This has been made possible by utilizing control of the pump speed down to 0.0003cc/s which then gives the operator the ability to test ultra-tight formations without the need for inflatable packers. By pulling down the pressure at an extremely low rate and using Artificial Intelligence to control the rate by knowing the formation rate, a proportional amount of volume can be extracted without calling it a tight test. During the operation by observing the rate, and making sure the pump is not oscillating, which indicates the formation rate is lower than the lowest rate the pump can withdraw, the test can be validated for formation flow and the pressure transient of the build – up can be analysed to confirm that at least spherical flow is observed. Once reservoir communication has been confirmed and by analysing drawdown and build-up pressure versus volume withdrawn and implementing the FRA equation, the reservoir pressure can be back calculated by considering isothermal compressibility and FRA slope. This paper highlights the best technical approach to quality check and quality control these tests, showing examples of various wells, where the technique has been used to predict a formation pressure, which can be used for further use for field development, drilling optimisation and production profiles. These pressures would never have been possible using static rates and volume.","PeriodicalId":153712,"journal":{"name":"SPWLA 62nd Annual Online Symposium Transactions","volume":"97 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2021-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"131535044","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
SPWLA 62nd Annual Online Symposium Transactions
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1