CN101261235B - 原油中含气率和含水率的双能χ射线测量方法 - Google Patents
原油中含气率和含水率的双能χ射线测量方法 Download PDFInfo
- Publication number
- CN101261235B CN101261235B CN2008100972031A CN200810097203A CN101261235B CN 101261235 B CN101261235 B CN 101261235B CN 2008100972031 A CN2008100972031 A CN 2008100972031A CN 200810097203 A CN200810097203 A CN 200810097203A CN 101261235 B CN101261235 B CN 101261235B
- Authority
- CN
- China
- Prior art keywords
- detector
- ray
- ray machine
- energy
- crude oil
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000010779 crude oil Substances 0.000 title claims abstract description 81
- 230000009977 dual effect Effects 0.000 title claims abstract description 32
- 238000000691 measurement method Methods 0.000 title description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 81
- 238000000034 method Methods 0.000 claims abstract description 52
- 239000003921 oil Substances 0.000 claims abstract description 39
- 238000001514 detection method Methods 0.000 claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 27
- 239000007789 gas Substances 0.000 claims abstract description 21
- 239000000463 material Substances 0.000 claims abstract description 17
- 238000007493 shaping process Methods 0.000 claims description 40
- 238000012937 correction Methods 0.000 claims description 36
- 230000003321 amplification Effects 0.000 claims description 34
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 34
- 239000011800 void material Substances 0.000 claims description 34
- 238000005259 measurement Methods 0.000 claims description 21
- 238000001228 spectrum Methods 0.000 claims description 21
- 238000012360 testing method Methods 0.000 claims description 19
- 239000008358 core component Substances 0.000 claims description 15
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000000151 deposition Methods 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 5
- 239000003345 natural gas Substances 0.000 claims description 5
- 239000011435 rock Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 239000012467 final product Substances 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 230000002285 radioactive effect Effects 0.000 abstract description 6
- 230000003993 interaction Effects 0.000 abstract description 5
- 239000000523 sample Substances 0.000 description 34
- 235000019198 oils Nutrition 0.000 description 24
- 230000000694 effects Effects 0.000 description 9
- 239000008186 active pharmaceutical agent Substances 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 238000013016 damping Methods 0.000 description 6
- 239000000047 product Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 230000009897 systematic effect Effects 0.000 description 5
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 229910052750 molybdenum Inorganic materials 0.000 description 4
- 239000011733 molybdenum Substances 0.000 description 4
- 230000005251 gamma ray Effects 0.000 description 3
- 238000002156 mixing Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 235000019476 oil-water mixture Nutrition 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000003723 Smelting Methods 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000005250 beta ray Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000005658 nuclear physics Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/26—Oils; Viscous liquids; Paints; Inks
- G01N33/28—Oils, i.e. hydrocarbon liquids
- G01N33/2823—Raw oil, drilling fluid or polyphasic mixtures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/083—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/12—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the material being a flowing fluid or a flowing granular solid
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- General Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Toxicology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种原油中含气率和含水率的双能χ射线测量方法,系利用χ光机产生的χ射线与物质作用原理,在油田生产中,油水气三相介质并存的条件下,测量输油管道中含水率和含气率指标。该测量方法依托于主要由三个大分系统组成的测量设备和一套专用软件,包括二种能量χ射线的产生分系统,一个或二套探测器构成的探测器分系统,以及一个总控和数据处理分系统组成,还包括准直器。专用软件中采用了特殊算法来求解含水率ω1和含气率ω3。本发明建立的理论模型精度较高,各种参数的物理意义比较明确,使用简单,还考虑了温度、压力等因素的影响。该系统还克服了放射源带来的重大安全隐患,特别适用于油田生产中自动在线的计量系统。
Description
技术领域
本发明涉及石油测量技术领域,具体涉及原油含气含水率的测量系统,具体涉及一种利用χ光机产生的χ射线与物质作用原理,油田生产中、油水气三相介质并存的条件下测量输油管道中含水率和含气率指标的双能χ射线测量技术,属于国际专利分类G01N技术领域。
背景技术
原油作为最重要的能源之一,从油井开采出的原油,是油、水、天然气等多相介质组成的混合物。处理这种混合物首先要进行气液分离,剩下的油水混合液体经脱水处理后得到含水率很低的成品油进行外输或销售。在原油脱水处理等一系列生产活动中,需要及时准确地掌握原油含水率和含气率等情况,以便控制生产过程,保证生产出合格的成品原油。因此,原油含水率和含气率等指标是石化行业石油采集、冶炼及运输过程中一个重要参数。特别是许多老油田,目前主要采用注水采油工艺,采出原油的含水率普遍偏高。因此,对原油进行含水率和含气率的准确检测在原油生产、贸易中有着重要作用。
目前测量原油含水率主要有以下几种方法:人工蒸馏化验法、微波法(或射频法)、电容法、短波法、导热法、振动密度计法和γ射线法。1、微波法(如CN1112677)是根据电磁波与介电物质相互作用,其耗散与物质的大小及相对介电常数有关,油和水的介电常数不同导致被测对象所呈现的射频阻抗特性不同,当射频信号传到以油水混合物为介质的电容式射频传感器时,其负载阻抗随着混合介质的不同油水比而变化,即当原油含水率变化时,波参量随之变化,从而实现含水率测量;2、电容法(如:CN1186236)是根据油水的介电常数不同,反映到由极板构成的电容器的电容量不同,测量电容量的变化,就可以测量含水率的变化;3、短波法(如:CN2349574)是利用一个(后)探头向原油中发射3.579MHz的短波信号,把当前原油状态查清,间隔几秒钟后,在通过另一个(前)探头向原油中发射3.579MHz的短波信号,又取出油中含水的信号,然后取二次测得的差值,经处理后可得出瞬时含水率;4、导热法(如:CN1259671)利用液-液两相流体的热物理性质的差异,如导热、比热、粘度等,同时测量原油的含水率和油水流量;5、振动密度计法(如:CN1789969,CN2359692)利用液位测量元件测量储油罐(或分离器)内原油的液位,压力测量仪表测量储油罐(或分离器)中无原油部分的压力和底部承受的压力,最后通过经验公式算出;6、射线法(如:CN86105543A,CN2359692Y,CN1086602A,CN2383068Y)是根据γ射线穿过不同介质时,其衰减不同的原理工作的。除γ射线法以外的其它各种测量方法,都属接触式测量,由于原油腐蚀性较强,结垢、结蜡严重,致使仪表长期运行的可靠性差,尤其是这些仪表都无法消除含气对含水率测量带来的影响,而导致了比较大的测量误差。对此,专利CN2452022Y、CN2646704Y和CN2646705Y专门设计了不同的擦除器来传感器外面累积的杂质。另外,电容法、射频法和微波法测量的含水率变化与被测量之间是非线性关系,在某一含水率范围内有拐点,而原油是油水气混合体,其物理化学性质多变,所以除γ射线法外的上面几种测量方法在实际应用中,都不能很好地满足生产要求。
根据γ射线与物质相互作用规律而工作的原油含水分析仪与混合流体的宏观流态和化学性质无关,能够对复杂的原油进行含水率和含气率的测量,深受油田的欢迎。
发明专利CN86105543A公开了一种利用放射源(如109Cd,或243Am等)发出的单能γ射线,对二相油水混合体的体积含水率的测量原理。实用新型专利CN2359692Y公布了一种利用238Pu放射源对二相油水混合体的含水率进行测量的装置。发明专利CN1086602A公开了一种在三相油水气混合体中,测量原油中含气、含水率的自动测量仪;在测量管道的侧壁上沿径向中心线对称位置两侧分别固定有γ射线源和透射探测器;在与γ射线源和透射探测器所在中心线成夹角且沿测量管道轴向与之相距一定距离的中心线侧壁上固定有散射探测器;最后根据测量的结果,经过数据处理而得到体积含气率和体积含水率。文献[1]对其测量原理,从理论做了进一步的探讨。实用新型专利CN2383068Y对依据上述原理设计的装置,做了改进,增加了一个搅拌装置,使实际上从油井里出来的油水气混合均匀,以便更进一步满足理论假设条件,以便提高测量的精度。不管怎么说,这个理论模型取近似的地方太多,各个参数物理意义不明确,对压力、温度等变化参数的影响,没有考虑修正,这些最终还是影响了其使用的方法和测量精度。
另外,使用放射源产生的γ射线测量法,还有一个比较大的弱点,就是存在放射性安全问题,特别是在目前反恐形势比较严峻的时期,这个弱点更加突出。
目前市场上尚未发现应用双能χ射线法测量原油中含气率和含水率的原理和装置。
发明内容
本发明的目的,在于针对三相油、水、气混合体中,实时高精度测量原油中含气、含水率的需求,提供一种采用双能技术测量含气、含水率的原理和方法。
本发明的目的是通过以下技术方案实现的。
一种原油中含气率和含水率的双能χ射线测量方法,系利用χ光机产生的χ射线与物质作用原理,在油田生产中,油水气三相介质并存的条件下,测量输油管道中含水率和含气率指标。该测量方法依托于主要由三个大分系统组成的测量设备和一套专用软件,即:二种能量χ射线的产生分系统,一个或二套探测器构成的探测器分系统,以及一个总控和数据处理分系统组成;其它装置还包括准直器;所述的探测器分系统包括探测器、前置放大器或光电倍增管、信号成形、放大、采样保持、AD转换等单元;总的控制和数据处理分系统包括数据的传输、同步、显示、控制和报警等部分;所述的专用软件中采用了特殊的算法来求解含水率ω1和含气率ω3;ω1表示水所占的重量百分比即含水率,ω3表示天然气所占的重量百分比即含气率,ω2表示油所占的重量百分比ω2=1-ω1-ω3。
在油水气三相状态下,该软件采用了如下的两个方程来求解含水率ω1和含气率ω3,
上述方程中:EH *、EL *分别代表χ光机的高能和低能χ射线所对应的等效能量;ρ代表油管中油水气三相状态下的实际密度,ρ1代表实际油管里所对应的温度、压力等条件下纯水的密度,ρ2代表实际油管里所对应的温度、压力等条件下纯原油的密度,ρ3代表实际油管里所对应的温度、压力等条件下纯天然气的密度;μ1、μ2、μ3分别代表纯水、纯原油、纯天然气在对应的等效射线能量下的线性衰减系数;x代表油管里,测试系统测量空间的线性厚度;N0(E*)代表在对应的等效射线能量条件下,油管里没有任何物质存在时,测试系统所测量的计数;N0(x,E*)代表在对应的测量厚度、等效射线能量条件下,测试系统所测量的计数;k、c分别为修正系数,同μ1、μ2、μ3一起,可以通过预先测量指数衰减曲线求得。
在只考虑油水二相状态时,可以把测量系统的χ光机能量简化为单能,这时所述软件中采用了如下公式来求解含水率ω1,
公式中:ρ代表油管中油水二相状态下的实际密度,ρ1代表实际油管里所对应的温度、压力等条件下纯水的密度,ρ2代表实际油管里所对应的温度、压力等条件下纯原油的密度;μ1、μ2分别代表纯水、纯原油在对应的等效射线能量下的线性衰减系数;x代表油管里,测试系统测量空间的线性厚度;N0代表在对应的等效射线能量条件下,油管里没有任何物质存在时,测试系统所测量的计数;N0(x)代表在对应的测量厚度、等效射线能量条件下,测试系统所测量的计数;k、c分别为修正系数,同μ1、μ2一起,可以通过预先测量指数衰减曲线求得。
设计χ光机系统时,将高、低能量之间满足的一定条件考虑进去,高能EH *与低能EL *的差别越大,测量精度越高;例如: 简单一点, 所述高能χ光机的能量范围可取在10keV~1MeV之间。
该测量方法依托于真双能χ光机测量系统,其测量设备核心部件的安装方式,是在原油管道的水平直径一端装设χ射线产生分系统,包括真双能χ光机控制系统,连接准直器和准直器,准直器屏蔽室中设有位置基本重合的高能χ射线的靶点和低能χ射线的靶点;原油管道水平直径的另一端装设探测器分系统,包括探测器和屏蔽管,连接总控和数据处理分系统,包括相互连接的第一路信号成形、放大和采样保持单元、第一组探测器的高压电源、第一路AD转换单元、第一路控制单元和数据处理计算机。
该测量方法依托于二台单能χ光机产生高低能χ射线,其测量设备核心部件的安装方式,是由并行排列的高能χ光机的控制系统和低能χ光机的控制系统,分别连接高能χ光机的靶点和高能χ光机准直器和屏蔽室、以及低能χ光机靶点和低能χ光机的准直器和屏蔽室,所述两台单能χ光机,分别对应两组探测器分系统——第一组探测器和屏蔽管,第二组探测器和屏蔽管,形成两条平行的探测分路,两组探测器分系统分别连接各自的总控和数据处理分系统,即:一组连接第一路信号成形、放大和采样保持单元、第一组探测器的高压电源、第一路AD转换单元、第一路控制单元,另一组连接第二路信号成形、放大和采样保持单元、第二组探测器的高压电源、第二路AD转换单元、第二路控制单元;两组数据均分别传至数据处理计算机。
该测量方法依托于伪双能χ光机测量系统,其测量设备核心部件的安装方式,包括一台单能χ光机控制系统,对应旋转机构,能谱预硬化装置同旋转机构安装在一起,外罩准直器,与之对应安装第一组探测器和屏蔽管;χ光机的靶点、准直器、探测器和屏蔽管组成一套χ射线探测通路,连接第一路信号成形、放大和采样保持单元、第一组探测器的高压电源、第一路AD转换单元、第一路控制单元和数据处理计算机。
该测量方法依托于二台单能χ光机测量系统,其测量设备核心部件的安装方式,采用了二台单能χ光机和探测器的横向安装方式,是在原油管道的上下左右四个端点部位进行安装,在其左端和下端分别装载高能χ光机的控制系统和低能χ光机的控制系统,分别配合有高能χ光机的靶点、高能χ光机的准直器和屏蔽室,以及低能χ光机的靶点和低能χ光机的准直器和屏蔽室,高能χ光机的靶点、准直器、与原油管道右端安装的探测器和屏蔽管组成一套高能χ射线探测通路,低能χ光机的靶点、准直器与原油管道右端安装的探测器和屏蔽管组成一套低能χ射线探测通路;两套探测通路处于原油管道的同一个横截面上,以减少测量设备的长度;两套探测通路的夹角可以变化,理想状态是,夹角越小越好,在原油管道横截面上均匀性比较好时,其夹角可以大一些。
该测量方法依托于伪双能探测器,其测量设备核心部件的安装方式,包括一台单能χ光机控制系统,外罩准直器,与之对应在探测器前面安装能谱过滤片,第一组探测器和屏蔽管;χ光机的靶点、准直器、能谱过滤片、探测器和屏蔽管组成一套高能χ射线探测通路,连接第一路信号成形、放大和采样保持单元、第一组探测器的高压电源、第一路AD转换单元、第一路控制单元和数据处理计算机。同时χ光机的靶点、准直器、第二组探测器和第二组探测器的屏蔽管组成一套低能χ射线探测通路,连接第二路信号成形、放大和采样保持单元、第二组探测器的高压电源、第二路AD转换单元、第二路控制单元和数据处理计算机
如果χ光机的束流不太稳定,为了提高测量精度,可以加入一路亮度校正探测器。其特征在于:在每个χ光机的出口安装一路亮度校正探测器或;亮度校正探测器的位置在χ光机的出口,但最好不要遮挡用于测量的主束流。以亮度校正探测器为例说明其安装关系:χ光机的靶点、准直器和亮度校正探测器组成一套χ射线束流大小探测通路,连接第三路信号成形、放大和采样保持单元、第三路探测器的高压电源、第三路AD转换单元、第三路控制单元和计算机。
本发明双能χ射线测量技术,即利用χ光机产生的χ射线与物质作用原理,在油田生产中,油水气三相介质并存的条件下,测量输油管道中含水率和含气率指标。该系统克服了放射源带来的重大安全隐患,特别适用于油田生产中,自动在线计量系统。本发明理论模型的精度比较高,各种参数的物理意义比较明确,使用简单,还能考虑温度、压力等因素的影响,特别适用于油田生产中,自动在线计量系统。当利用X光机作为射线源时,可以免除丢失放射源的困扰,提高了辐射防护的安全系数,从根本上杜绝了恐怖分子获取脏弹原料的机会,对国家安全有着特别重要的意义。
附图说明
图1-1为20keV钼靶χ光机的能谱图;
图1-2为30keV钼靶χ光机的能谱图;
图1-3为35keV钼靶χ光机的能谱图;
图1-4为50keV钼靶χ光机的能谱图;
图2为真双能χ光机和探测器安装方式示意图;
图3为二台单能χ光机和探测器纵向安装方式示意图;
图4为伪双能χ光机和探测器安装方式示意图;
图5为二台单能χ光机和探测器横向安装方式示意图;
图6为单能X光机和伪双能探测器安装方式示意图(二组探测器也可在其它位置,如在同一横截面上)。
图中:1-测量设备;2-高能χ光机的靶点;3-高能χ光机的准直器和屏蔽室;4-(第一组)探测器;5-(第一组)探测器的屏蔽管;6-低能χ光机的靶点;7-低能χ光机的准直器和屏蔽室;8-第二组探测器;9-第二组探测器的屏蔽管;10-原油管道;12-能谱预硬化装置;13-旋转机构;20-(第一路)信号成形、放大和采样保持单元;21-第二路信号成形、放大和采样保持单元;22-(第一组)探测器的高压电源;23-第二组探测器的高压电源;24-(第一路)AD转换单元;25-第二路AD转换单元;26-(第一路)控制单元;27-第二路控制单元;28-计算机;30-真双能χ光机的控制系统;31-高能χ光机的控制系统;32-低能χ光机的控制系统;33-(第一路)亮度校正探测器;34-第三路信号成形、放大和采样保持单元;35-第三路探测器的高压电源;36-第三路AD转换单元;37-第三路控制单元;38-第二路亮度校正探测器;39-第四路信号成形、放大和采样保持单元;40-第四路探测器的高压电源;41-第四路AD转换单元;42-第四路控制单元;50-能谱过滤片。
具体实施方式
本发明技术方案以如下方式实现:
对单能γ射线来讲,与物质的相互作用遵循指数衰减规律,即公式(1)成立。
其中:N0 ----射线穿过空气后,测量的计数。
N(xm)----射线穿过质量厚度为xm的物质后,测量的计数。
xm ----射线穿过物质的质量厚度。
μm ----射线穿过物质的质量衰减系数。
由于χ光机产生的χ射线的能谱是连续的,文献[3]给出了几种能量χ光机产生的能谱,如图1所示。具有连续能谱的的χ射线,其与物质的相互作用是否还遵循指数衰减规律呢?对这个问题的理论研究,参见文献[4]。
参见说明书末尾引述的文献[4],推导了下列命题:
这里,将[c,d]m等分,分点为
在分点xj处有
记:
R=max{|Rj| j=0,1,2,Λ,m}
经过一番推导后,得出:
在本课题的实际应用中,相当于实测数据;ke-αx相当于计算拟合数据;f(x)相当于在x点对应的绝对误差;R相当于所有绝对误差中最大误差的绝对值;G相当于实测数据中最大的;α相当于等效质量吸收系数μm;αi相当于射线能量为Ei所对应的质量吸收系数μm(Ei);x相当于样品的质量厚度xm;Δx相当于选取实验数据的间隔Δxm(自变量)。理论上,Δx可以取得很小,这时|f(x)|≤R。R可以通过设计实验获得,从而得到整个函数的误差估计。如果所得到的最大误差可以接受的话,其衰减规律就可以用指数衰减规律来近似。
根据理论分析和文献[3]的实验数据发现,在吸收体不是太厚的情况下,其与物质的相互作用近似遵循指数衰减规律,如果太厚,误差就大了一点。对具体的χ光机,首先应该用实验测试一下,找出近似遵循指数衰减规律的条件。由此再依据测试样品的最大厚度,确定χ光机的最低能量EL *(对应于低能χ光机的管高压VL)。根据经验,建议依据 确定EH *(对应于高能χ光机的管高压VH)的值,。
因此,不妨假设所讨论的χ光机的χ射线,在我们所讨论的测试样品的厚度范围内,仍然近似遵循指数衰减规律,即公式(1)近似成立。为了理论模型与实验数据更好的符合,特加入二个拟合系数k和c,如公式(5)所示。
注意:
(1)这时在公式(5)中,应采用等效能量E*的概念。
(2)k和c的理论值为:k=1;c=0。在缺乏实验值时,可直接引用理论值。
(3)μm(E*)、k和c,可以事先通过实验测得。
(4)为了书写方便,下面中的μm(E*),均简写为μm。
式中ω1,ω2,...,ωi...,ωN分别为组成元素的重量百分比。
注意:元素的质量衰减系数表示为μm,线性厚度为x,质量厚度为xm。即:
μm=μ/ρ xm=x·ρ
1、考虑油水二相状态(无下标-原油+水混合状态,下标1-纯水状态,下标2-纯原油状态)
μm=ω1μm1+ω2μm2
=ω1μm1+(1-ω1)μm2
=ω1(μm1-μm2)+μm2 (7)
Θμm=μ/ρ xm=x·ρ
∴
把(8)代入(5),并化简得:
即:
在这种情况下,一个未知量含水率ω1,一个方程,故采用一台单能χ光机即可解决测量问题。
2、考虑油水气三相状态(无下标-原油+水+气混合状态,下标1-纯水状态,下标2-纯原油状态,下标3-纯天然气状态)
μm=ω1μm1+ω2μm2+ω3μm3
=ω1μm1+(1-ω1-ω3)μm2+ω3μm3
=ω1(μm1-μm2)+ω3(μm3-μm2)+μm2 (10)
Θμm=μ/ρ xm=x·ρ
∴
把(11)代入(5),并化简得:
为了求得含水率ω1和含气率ω3,需要列出类似(12)的二个方程。从核物理上,可以通过二种不同能量χ光机的χ射线来测量求得。
本文中,采用双能测量模式讨论如下:
设EH *、EL *分别代表χ光机的高能和低能χ射线所对应的等效能量,则(12)可以表示为:
理论上,可根据方程(13)和(14)求得ω1和ω3,这就是油水气三相状态下的双能χ射线测量含水率和含气率的原理。
注意:
1)在实验测量ρ1、ρ2、ρ3和ρ时,需要同时检测样品的温度、压力等参数的影响。
2)因为气体的状态与温度、压力密切相关,应用中,要测量与实际条件相一致的ρ3和μm3。
3)求解方程时,要利用查表法,采用与实际条件相对应的ρ值,可以通过实时测量得到。
4)高能EH *与低能EL *的差别越大,测量精度要越好。例如: 简单一点,
5)公式中的各种系数,μ1、μ2、μ3、k和c等,可以在实验室中,分别用高、低能χ光机,照射不同质量厚度的标定介质(纯原油、纯水、纯天然气),用衰减测量法获得的实验数据,再按最小二乘法拟合求得。注意:k和c可以用原油所对应的值来近似,也可以针对各种情况,在实验室测出数据,建成一个数据库,在现场使用时,用查表法获取该数据。最后根据本发明推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
下面结合附图和具体的实施方式对本发明作进一步的描述。
参见各附图,本发明的总体设计思想:本测量装置1主要由三个大分系统组成。二种能量的χ光机构成χ射线产生分系统;一个或二套探测器构成探测器分系统;和一个总的控制和数据处理分系统组成。其它分系统部分还包括准直器,或其它机构等。探测器分系统包括探测器、前置放大器或光电倍增管、信号成形、放大、采样保持、AD转换等单元。数据的传输、同步、显示、各部分控制和报警等工作由总的控制和数据处理分系统统一完成。
图中标号含义:1-测量设备;2-高能χ光机的靶点;3-高能χ光机的准直器和屏蔽室;4-(第一组)探测器;5-(第一组)探测器的屏蔽管;6-低能χ光机的靶点;7-低能χ光机的准直器和屏蔽室;8-第二组探测器;9-第二组探测器的屏蔽管;10-原油管道;12-能谱预硬化装置;13-旋转机构;20-(第一路)信号成形、放大和采样保持单元;21-第二路信号成形、放大和采样保持单元;22-(第一组)探测器的高压电源;23-第二组探测器的高压电源;24-(第一路)AD转换单元;25-第二路AD转换单元;26-(第一路)控制单元;27-第二路控制单元;28-计算机;30-真双能χ光机的控制系统;31-高能χ光机的控制系统;32-低能χ光机的控制系统;33-(第一路)亮度校正探测器;34-第三路信号成形、放大和采样保持单元;35-第三路探测器的高压电源;36-第三路AD转换单元;37-第三路控制单元;38-第二路亮度校正探测器;39-第四路信号成形、放大和采样保持单元;40-第四路探测器的高压电源;41-第四路AD转换单元;42-第四路控制单元;50-能谱过滤片。
根据本发明的原理,给出以下四种应用实例:
1、真双能χ光机测量系统的工作原理
测量设备1的核心部件的安装位置如图2所示。
真双能χ光机的特征是高能χ射线的靶点2和低能χ射线的靶点6的位置基本重合,真双能χ光机的控制系统30控制高、低能χ射线的分时交替输出,输出脉冲的频率与流体的速度有关。理论上希望高、低能χ射线能同时打在介质的同一个位置上,实践中,可根据介质的均匀度、流速、要求监测数据的间隔等,调节其参数,以保证测试条件尽可能满足理论模型和误差要求。
由靶点2和6发出的χ射线,经过准直器3和7(二者已合为一体,为了和下文其它实例中的编号一致,特保留两个标号)后,穿过原油管道10中的介质,被探测器4转换成电信号。探测器的屏蔽管5的作用是保护探测器4和减少探测本底、散射信号对探测器4的影响。
探测器的高压电源22给探测器4提供工作电压,探测器4的信号输出到信号成形、放大和采样保持单元20,信号经过放大、处理后,送到AD转换单元24转换成数字信号,最后送到计算机28进行分析处理。控制单元26用来同步、协调各个单元或分系统的工作。
在本实例中省了第二套探测器系统,高、低能信号的识别、同步是通过真双能χ光机的控制系统30和控制单元26的信号交互来实现的。
如果χ光机的束流随时间的变化较大,即束流不稳定,则在实际处理数据时,应该对其修正。为了获得χ光机的束流随时间的变化量,需要增加用于亮度校正的探测系统。即:在χ光机的出口安装一路亮度校正探测器33,第三路探测器的高压电源35给探测器33提供高压,探测器33的信号输出到第三路信号成形、放大和采样保持单元34,信号经过放大、处理后,送到第三路AD转换单元36转换成数字信号,最后送到计算机28进行分析处理。控制单元37用来同步、协调各个单元或分系统的工作。
如果χ光机的束流稳定,对系统测量带来的误差,可以忽略,则有关亮度校正的探测系统可省略。
计算机28上的专用软件,先把探测到的数据分成高能组数据系列和低能组数据系列,应用亮度探测器获得的高、低能数据,分别对其对应时刻的所有数据进行修正,首先消除χ光机的束流随时间的变化的影响。然后,应用本发明中推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
2、二台单能χ光机和探测器纵向安装方式测量系统的工作原理
测量设备1的核心部件的安装位置如图3所示。
本应用实例的特点是,利用目前现有的二台单能χ光机来产生高、低能χ射线。高能χ射线的靶点2、准直器3、探测器4和屏蔽管5等组成一套高能χ射线探测通路,同理低能χ射线的靶点6、准直器7、探测器8和屏蔽管9等组成一套低能χ射线探测通路。两套探测通路越靠近,越能满足理论上希望高低能χ射线能同时打在介质的同一个位置上的理想要求。实践中,可根据介质的均匀度、流速、要求监测数据的间隔等,调节其参数,以保证测试条件尽可能满足理论模型和误差要求。
高能χ光机的控制系统31控制高能χ光机的靶点2发出高能χ射线,经过准直器和屏蔽室3后,穿过原油管道10中的介质,被第一组探测器4转换成电信号。第一组探测器的屏蔽管5的作用是保护第一组探测器4,同时减少本底、散射信号对第一组探测器4的影响。第一组探测器的高压电源22给第一组探测器4提供工作电压,第一组探测器4的信号输出到第一路信号成形、放大和采样保持单元20,信号经过放大、处理后,送到第一路AD转换单元24转换成数字信号,最后送到计算机28进行分析处理。第一路控制单元26用来同步、协调第一路各个单元和与其相关的分系统的工作。
同理低能χ光机的控制系统32控制低能χ光机的靶点6发出低能χ射线,经过准直器和屏蔽室7后,穿过原油管道10中的介质,被第二组探测器8转换成电信号。第二组探测器的屏蔽管9的作用是保护第二组探测器8,同时减少本底、散射信号对第二组探测器8的影响。第二组探测器的高压电源23给第二组探测器8提供工作电压,第二组探测器8的信号输出到第二路信号成形、放大和采样保持单元21,信号经过放大、处理后,送到第二路AD转换单元25转换成数字信号,最后送到计算机28进行分析处理。第二路控制单元27用来同步、协调第二路各个单元和与其相关的分系统的工作。
在本实例中,降低了对χ光机的设计要求,只要使用目前市场上的产品即可。同一位置上的高、低信号的可以用流体的速度和探测时间来加以同步。
如果χ光机的束流随时间的变化较大,实际处理数据时,应该对其修正。为了获得高能χ光机的束流随时间的变化量,需要增加用于亮度校正的探测系统。即:在高能χ光机的出口安装第一路亮度校正探测器33,第三路探测器的高压电源35给探测器33提供高压,探测器33的信号输出到第三路信号成形、放大和采样保持单元34,信号经过放大、处理后,送到第三路AD转换单元36转换成数字信号,最后送到计算机28进行分析处理。控制单元37用来同步、协调各个单元或分系统的工作。
同理,为了获得低能χ光机的束流随时间的变化量,在低能χ光机的出口安装第二路亮度校正探测器38,第四路探测器的高压电源40给探测器38提供高压,探测器38的信号输出到第四路信号成形、放大和采样保持单元39,信号经过放大、处理后,送到第四路AD转换单元41转换成数字信号,最后送到计算机28进行分析处理。第四路控制单元42用来同步、协调各个单元或分系统的工作。
如果χ光机的束流稳定,对系统测量带来的误差,可以忽略,则有关亮度校正的探测系统可省略。
计算机28上的专用软件,先把探测到高能组数据和低能组数据,应用亮度探测器获得的高、低能数据,分别对其对应时刻的所有数据进行修正,首先消除χ光机的束流随时间的变化的影响。然后,应用本发明中推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
3、伪双能χ光机测量系统的工作原理
测量设备1的核心部件的安装位置如图4所示。
本应用实例的特点是,利用目前现有的一台单能χ光机通过分时预硬化技术来产生高、低能χ射线。χ光机的靶点2、准直器3、探测器4和屏蔽管5组成一套χ射线探测通路。能谱预硬化装置12同旋转机构13安装在一起,旋转机构13的旋转中心不能遮挡χ光机的靶点2。当旋转机构13旋转的位置使能谱预硬化装置12挡住靶点2发出的χ射线时,这时从准直器3射出来的χ射线,就是高能χ射线,系统当作高能测量系统。否则,系统当作低能测量系统使用。实践中,可根据介质的均匀度、流速、要求监测数据的间隔等,调节旋转机构13旋转速度,以保证测试条件尽可能满足理论模型和误差要求。
由靶点2发出的χ射线,经过准直器3后,穿过原油管道10中的介质,被探测器4转换成电信号。探测器的屏蔽管5的作用是保护探测器4,同时减少探测本底、散射信号对探测器4的影响。探测器的高压电源22给探测器4提供工作电压,探测器4的信号输出到信号成形、放大和采样保持单元20,信号经过放大、处理后,送到AD转换单元24转换成数字信号,最后送到计算机28进行分析处理。控制单元26用来同步、协调各个单元或分系统的工作,特别是高、低能的识别是通过χ光机的控制系统30和控制单元26的信号交互来实现的。
在本实例中,降低了对χ光机的设计要求,只要使用目前市场上的产品即可。
如果χ光机的束流随时间的变化较大,实际处理数据时,应该对其修正。为了获得χ光机的束流随时间的变化量,需要增加用于亮度校正的探测系统。即:在χ光机的出口安装一路亮度校正探测器33,第三路探测器的高压电源35给探测器33提供高压,探测器33的信号输出到第三路信号成形、放大和采样保持单元34,信号经过放大、处理后,送到第三路AD转换单元36转换成数字信号,最后送到计算机28进行分析处理。控制单元37用来同步、协调各个单元或分系统的工作。
如果χ光机的束流稳定,对系统测量带来的误差,可以忽略,则有关亮度校正的探测系统可省略。
计算机28上的专用软件,先把探测到的数据分成高能组数据系列和低能组数据系列,应用亮度探测器获得的高、低能数据,分别对其对应时刻的所有数据进行修正,首先消除χ光机的束流随时间的变化的影响。然后,应用本发明中推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
4、二台单能χ光机和探测器横向安装方式测量系统的工作原理
测量设备1的核心部件的安装位置如图5所示。
本应用实例的特点是,利用目前现有的二台单能χ光机来产生高、低能χ射线。高能χ光机的靶点2、准直器3、探测器4和屏蔽管5组成一套高能χ射线探测通路,同理低能χ光机的靶点6、准直器7、探测器8和屏蔽管9组成一套低能χ射线探测通路。两套探测通路在原油管道10的同一个横截面上,这样减少了测量设备1的长度。注意图5中所示的两套探测通路的夹角不一定要求90°,只要能完整安装好测试部件即可。本实例要求流体在原油管道10的同一个横截面分布近似相同,这样才能满足理论上希望高低能χ射线能同时打在介质的同一个位置上的理想要求。实际应用中,可在介质流入测试设备前,采取措施对流体加以搅拌,使之混合均匀即可。
高能χ光机的控制系统31控制高能χ光机的靶点2发出高能χ射线,经过准直器和屏蔽室3后,穿过原油管道10中的介质,被第一组探测器4转换成电信号。第一组探测器的屏蔽管5的作用是保护第一组探测器4,同时减少本底、散射信号对第一组探测器4的影响。第一组探测器的高压电源22给第一组探测器4提供工作电压,第一组探测器4的信号输出到第一路信号成形、放大和采样保持单元20,信号经过放大、处理后,送到第一路AD转换单元24转换成数字信号,最后送到计算机28进行分析处理。第一路控制单元26用来同步、协调第一路各个单元和与其相关的分系统的工作。
同理低能χ光机的控制系统32控制低能χ光机的靶点6发出低能χ射线,经过准直器和屏蔽室7后,穿过原油管道10中的介质,被第二组探测器8转换成电信号。第二组探测器的屏蔽管9的作用是保护第二组探测器8,同时减少本底、散射信号对第二组探测器8的影响。第二组探测器的高压电源23给第二组探测器8提供工作电压,第二组探测器8的信号输出到第二路信号成形、放大和采样保持单元21,信号经过放大、处理后,送到第二路AD转换单元25转换成数字信号,最后送到计算机28进行分析处理。第二路控制单元27用来同步、协调第二路各个单元和与其相关的分系统的工作。
在本实例中,降低了对χ光机的设计要求,只要使用目前市场上的产品即可。同一位置上的高、低信号的可以采用混合流体的措施,使流体的横截面保持均匀,以此逼近理论模型的条件。
如果χ光机的束流随时间的变化较大,实际处理数据时,应该对其修正。为了获得高能χ光机的束流随时间的变化量,需要增加用于亮度校正的探测系统。即:在高能χ光机的出口安装第一路亮度校正探测器33,第三路探测器的高压电源35给探测器33提供高压,探测器33的信号输出到第三路信号成形、放大和采样保持单元34,信号经过放大、处理后,送到第三路AD转换单元36转换成数字信号,最后送到计算机28进行分析处理。控制单元37用来同步、协调各个单元或分系统的工作。
同理,为了获得低能χ光机的束流随时间的变化量,在低能χ光机的出口安装第二路亮度校正探测器38,第四路探测器的高压电源40给探测器38提供高压,探测器38的信号输出到第四路信号成形、放大和采样保持单元39,信号经过放大、处理后,送到第四路AD转换单元41转换成数字信号,最后送到计算机28进行分析处理。第四路控制单元42用来同步、协调各个单元或分系统的工作。
如果χ光机的束流稳定,对系统测量带来的误差,可以忽略,则有关亮度校正的探测系统可省略。
计算机28上的专用软件,先把探测到的高能组数据系列和低能组数据系列,应用亮度探测器获得的高、低能数据,分别对其对应时刻的所有数据进行修正,首先消除χ光机的束流随时间的变化的影响。然后,应用本发明中推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
5、单能X光机和伪双能探测器安装方式
测量设备1的核心部件的安装位置如图6所示。注意:二路探测器也可安装在其它位置,如在原油管道10同一横截面上。
本应用实例的特点是,利用目前现有的一台单能χ光机,通过对一路探测器采用预硬化技术来产生高能χ射线探测通路,另一路探测器则为低能χ射线探测通路。能谱过滤片50安装在探测器4的前面,χ光机的靶点2、准直器3、能谱过滤片50、探测器4和屏蔽管5组成一套高能χ射线探测通路。这时这一路系统当作高能测量系统。χ光机的靶点2、准直器3、第二组探测器8和第二组探测器的屏蔽管9组成另一套低能χ射线探测通路,这一路系统当作低能测量系统使用。
第一路:由靶点2发出的χ射线,经过准直器3后,穿过原油管道10中的介质,经过能谱过滤片50把能谱预硬化后,变成高能能谱,被探测器4转换成电信号。探测器的屏蔽管5的作用是保护探测器4,同时减少探测本底、散射信号对探测器4的影响。探测器的高压电源22给探测器4提供工作电压,探测器4的信号输出到信号成形、放大和采样保持单元20,信号经过放大、处理后,送到AD转换单元24转换成数字信号,最后送到计算机28进行分析处理。控制单元26用来同步、协调各个单元或分系统的工作。
第二路:由靶点2发出的χ射线,经过准直器3后,穿过原油管道10中的介质,被第二组探测器8转换成电信号。第二组探测器的屏蔽管9的作用是保护探测器8,同时减少探测本底、散射信号对探测器8的影响。第二组探测器的高压电源23给探测器8提供工作电压,探测器8的信号输出到第二路信号成形、放大和采样保持单元21,信号经过放大、处理后,送到第二路AD转换单元25转换成数字信号,最后送到计算机28进行分析处理。第二路控制单元27用来同步、协调各个单元或分系统的工作。
在本实例中,降低了对χ光机的设计要求,只要使用目前市场上的产品即可。
如果χ光机的束流随时间的变化较大,实际处理数据时,应该对其修正。为了获得χ光机的束流随时间的变化量,需要增加用于亮度校正的探测系统。即:在χ光机的出口安装一路亮度校正探测器33,第三路探测器的高压电源35给探测器33提供高压,探测器33的信号输出到第三路信号成形、放大和采样保持单元34,信号经过放大、处理后,送到第三路AD转换单元36转换成数字信号,最后送到计算机28进行分析处理。控制单元37用来同步、协调各个单元或分系统的工作。
如果χ光机的束流稳定,对系统测量带来的误差,可以忽略,则有关亮度校正的探测系统可省略。
计算机28上的专用软件,先把探测到的数据分成高能组数据系列和低能组数据系列,应用亮度探测器获得的高、低能数据,分别对其对应时刻的所有数据进行修正,首先消除χ光机的束流随时间的变化的影响。然后,应用本发明中推导的模型(也可以采用其它合适的模型),算出原油中的含水率、含气率等指标。
参考文献:
[1]任晓峰,李建等,“射线体积型原油含气、含水率分析仪的原理及应用”,自动化仪表,Vol.28,No.10,2007。
[2]李星洪,《辐射防护基础》,原子能出版社,1982年。
[3]何承发,巴维真,吾勤之等,“衰减测量法确定X光机能谱”,第九届全国核电子学与核探测技术学术年会论文集,1998年。
[4]罗平安,缪卫东,朱中梅等,“β射线的吸收规律研究”,北京核学会‘99青年核科技论文报告会文集,1999。(Luo Pingan,Miao Weidong,Zhu Zhongmei etc.“The Research on β Rays AbsorbingLaw”,’99 Youth Nuclear Technological Conference ofthe Beijing Nuclear Academy,1999.)
Claims (9)
1.一种原油中含气率和含水率的双能x射线测量方法,系利用x光机产生的x射线与物质作用原理,在油田生产中,油水气三相介质并存的条件下,测量输油管道中含水率和含气率指标,其特征在于:该测量方法依托于主要由三个大分系统组成的测量设备(1)和一套软件,即:二种能量x射线的产生分系统,一套或二套探测器分系统,以及一个总控和数据处理分系统组成;其它装置还包括准直器;所述的探测器分系统包括探测器、前置放大器或光电倍增管、以及信号成形、放大、采样保持单元、AD转换单元;总控和数据处理分系统包括数据的传输、同步、显示、控制和报警部分;所述的软件中采用了本发明提供的模型算法来求解含水率ω1和含气率ω3;ω1表示水所占的重量百分比即含水率,ω3表示天然气所占的重量百分比即含气率,ω2表示油所占的重量百分比ω2=1-ω1-ω3;
本发明提供的模型算法,其特征在于:在油水气三相状态下,采用了如下的两个方程来求解含水率ω1和含气率ω3,
上述方程中:分别代表x光机的高能和低能x射线所对应的等效能量;ρ代表油管中油水气三相状态下的实际密度,ρ1代表实际油管里所对应的温度、压力条件下纯水的密度,ρ2代表实际油管里所对应的温度、压力条件下纯原油的密度,ρ3代表实际油管里所对应的温度、压力条件下纯天然气的密度;μ1、μ2、μ3分别代表纯水、纯原油、纯天然气在对应的等效射线能量下的线性衰减系数;x代表油管里,测试系统测量空间的线性厚度;N0(E*)代表在对应的等效射线能量条件下,油管里没有任何物质存在时,测试系统所测量的计数;N(x,E*)代表在对应的测量厚度、等效射线能量条件下,测试系统所测量的计数;k、c分别为修正系数,同μ1、μ2、μ3一起,通过预先测量指数衰减曲线求得。
2.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:在只考虑油水二相状态时,把测量系统的x光机能量简化为单能,这时所述软件中采用了如下公式来求解含水率ω1,
公式中:ρ代表油管中油水二相状态下的实际密度,ρ1代表实际油管里所对应的温度、压力条件下纯水的密度,ρ2代表实际油管里所对应的温度、压力条件下纯原油的密度;μ1、μ2分别代表纯水、纯原油在对应的等效射线能量下的线性衰减系数;x代表油管里,测试系统测量空间的线性厚度;N0代表在对应的等效射线能量条件下,油管里没有任何物质存在时,测试系统所测量的计数;N(x)代表在对应的测量厚度、等效射线能量条件下,测试系统所测量的计数;k、c分别为修正系数,同μ1、μ2一起,通过预先测量指数衰减曲线求得。
4.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:该测量方法依托于真双能x光机测量系统,其测量设备(1)核心部件的安装方式,是在原油管道(10)的水平直径一端装设x射线产生分系统,包括真双能x光机控制系统(30),连接准直器(3),准直器屏蔽室中设有位置基本重合的高能x射线的靶点(2)和低能x射线的靶点(6);原油管道(10)水平直径的另一端装设探测器分系统,包括第一组探测器(4)和屏蔽管(5),连接第一路信号成形、放大和采样保持单元(20)、第一组探测器的高压电源(22)、第一路AD转换单元(24)、第一路控制单元(26)和数据处理计算机(28)。
5.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:该测量方法依托于二台单能x光机产生高低能x射线,其测量设备(1)核心部件的安装方式,是由并行排列的高能x光机的控制系统(31)和低能x光机的控制系统(32),分别连接高能x光机的靶点(2)和高能x光机准直器和屏蔽室(3)、以及低能x光机靶点(6)和低能x光机的准直器和屏蔽室(7),所述两台单能x光机,分别对应两套探测器分系统——第一组探测器(4)和屏蔽管(5),第二组探测器(8)和屏蔽管(9),形成两条平行的探测分路,两套探测器分系统分别连接各自的总控和数据处理分系统,即:一套探测器分系统的第一组探测器(4)连接第一路信号成形、放大和采样保持单元(20)、第一组探测器的高压电源(22)、第一路AD转换单元(24)、第一路控制单元(26),另一套探测器分系统的第二组探测器(8)连接第二路信号成形、放大和采样保持单元(21)、第二组探测器的高压电源(23)、第二路AD转换单元(25)、第二路控制单元(27);两组数据均分别传至数据处理计算机(28)。
6.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:该测量方法依托于伪双能x光机测量系统,其测量设备(1)核心部件的安装方式,包括一台单能x光机控制系统(30),对应旋转机构(13),能谱预硬化装置(12)同旋转机构(13)安装在一起,外罩准直器(3),与准直器(3)对应安装的第一组探测器(4)和屏蔽管(5);x光机的靶点(2)、准直器(3)、第一组探测器(4)和屏蔽管(5)组成一套x射线探测通路,连接第一路信号成形、放大和采样保持单元(20)、第一组探测器的高压电源(22)、第一路AD转换单元(24)、第一路控制单元(26)和数据处理计算机(28)。
7.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:该测量方法依托于二台单能x光机测量系统,其测量设备(1)核心部件的安装方式,采用了二台单能x光机和探测器的横向安装方式,是在原油管道(10)的上下左右四个端点部位进行安装,在其左端和下端分别装载高能x光机的控制系统(31)和低能x光机的控制系统(32),分别配合有高能x光机的靶点(2)、高能x光机的准直器和屏蔽室(3),以及低能x光机的靶点(6)和低能x光机的准直器和屏蔽室(7),高能x光机的靶点(2)、准直器(3)、与原油管道(10)右端安装的第一组探测器(4)和屏蔽管(5)组成一套高能x射线探测通路,低能x光机的靶点(6)、准直器(7)与原油管道(10)右端安装的第二组探测器(8)和屏蔽管(9)组成一套低能x射线探测通路;两套探测通路处于原油管道(10)的同一个横截面上,以减少测量设备(1)的长度;两套探测通路的夹角可以变化,理想状态是,夹角越小越好,在原油管道(10)横截面上均匀性比较好时,其夹角可以大一些。
8.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,其特征在于:该测量方法依托于伪双能探测器,其测量设备(1)核心部件的安装方式,包括一台单能x光机控制系统(30),外罩准直器(3),第一组探测器(4)和屏蔽管(5),且在第一组探测器(4)前面安装能谱过滤片(50);x光机的靶点(2)、准直器(3)、能谱过滤片(50)、第一组探测器(4)和屏蔽管(5)组成一套高能x射线探测通路,连接第一路信号成形、放大和采样保持单元(20)、第一组探测器的高压电源(22)、第一路AD转换单元(24)、第一路控制单元(26)和数据处理计算机(28),同时x光机的靶点(2)、准直器(3)、第二组探测器(8)和第二组探测器的屏蔽管(9)组成一套低能x射线探测通路,连接第二路信号成形、放大和采样保持单元(21)、第二组探测器的高压电源(23)、第二路AD转换单元(25)、第二路控制单元(27)和数据处理计算机(28)。
9.根据权利要求1所述的原油中含气率和含水率的双能x射线测量方法,如果x光机的束流不太稳定,为了提高测量精度,加入一路亮度校正探测器;其特征在于:在以二台单能x光机来产生高能和低能的x光机出口各安装一路亮度校正探测器(33)、(38),在真双能和伪双能x光机上只需要安装一路亮度校正探测器(33)即可;亮度校正探测器的位置在x光机的出口,但最好不要遮挡用于测量的主束流,如果遮挡了主束流,则会减少实际的输出剂量,并改变能谱成分的构成,对引起的影响要注意修正;以安装一路亮度校正探测器(33)为例说明其安装关系:x光机的靶点、准直器和亮度校正探测器(33)组成一套x射线束流大小探测通路,连接第三路信号成形、放大和采样保持单元(34)、第三组探测器的高压电源(35)、第三路AD转换单元(36)、第三路控制单元(37)和计算机(28)。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100972031A CN101261235B (zh) | 2008-05-06 | 2008-05-06 | 原油中含气率和含水率的双能χ射线测量方法 |
PCT/CN2009/070253 WO2009135390A1 (zh) | 2008-05-06 | 2009-01-21 | 多相流中物质含量的测量方法和系统 |
US12/991,170 US8718230B2 (en) | 2008-05-06 | 2009-01-21 | Method and system for determining the constituent content of a multiphase fluid |
RU2010146142/28A RU2466383C2 (ru) | 2008-05-06 | 2009-01-21 | Способ и система для определения содержания компонентов в многофазном флюиде |
EP09741665.5A EP2275804B1 (en) | 2008-05-06 | 2009-01-21 | A method and system for determining the constituent content of a multiphase fluid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100972031A CN101261235B (zh) | 2008-05-06 | 2008-05-06 | 原油中含气率和含水率的双能χ射线测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101261235A CN101261235A (zh) | 2008-09-10 |
CN101261235B true CN101261235B (zh) | 2010-12-08 |
Family
ID=39961829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100972031A Expired - Fee Related CN101261235B (zh) | 2008-05-06 | 2008-05-06 | 原油中含气率和含水率的双能χ射线测量方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8718230B2 (zh) |
EP (1) | EP2275804B1 (zh) |
CN (1) | CN101261235B (zh) |
RU (1) | RU2466383C2 (zh) |
WO (1) | WO2009135390A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107085089A (zh) * | 2017-05-21 | 2017-08-22 | 魏建军 | 一种原油含水率实时测量系统 |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101261235B (zh) | 2008-05-06 | 2010-12-08 | 罗平安 | 原油中含气率和含水率的双能χ射线测量方法 |
CN101852746A (zh) * | 2010-06-29 | 2010-10-06 | 北京中油联自动化技术开发有限公司 | 多相流中含水率与含气率的x射线测量装置及其测量方法 |
RU2569909C2 (ru) * | 2011-09-20 | 2015-12-10 | Сименс Акциенгезелльшафт | Устройство для измерения состава потока многофазной смеси |
EP2574919B1 (en) * | 2011-09-29 | 2014-05-07 | Service Pétroliers Schlumberger | Apparatus and method for fluid phase fraction determination using X-rays |
CN102749111B (zh) * | 2012-07-24 | 2013-12-25 | 兰州海默科技股份有限公司 | 湿气流量测量方法及其装置 |
CN102879410B (zh) * | 2012-08-22 | 2014-07-30 | 北京乾达源科技有限公司 | 一种在线测量油水气多相流含率的方法及装置 |
US20150226589A1 (en) * | 2012-08-27 | 2015-08-13 | Siemens Aktiengesellschaft | X-Ray Based Multiphase Flow Meter with Energy Resolving Matrix Detector |
CN102967609B (zh) * | 2012-11-28 | 2015-05-13 | 中国石油大学(北京) | 多相流体饱和度分布测量油藏渗流实验模型和测量方法 |
RU2558001C2 (ru) * | 2013-06-19 | 2015-07-27 | Закрытое акционерное общество "Научно-исследовательский институт интроскопии МНПО "СПЕКТР" | Способ двухэнергетической оценки средней плотности и эффективного атомного номера многокомпонентных материалов |
RU2530460C1 (ru) * | 2013-07-05 | 2014-10-10 | Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт автоматики им. Н.Л. Духова" (ФГУП "ВНИИА") | Анализатор многофазной жидкости |
RU2632249C1 (ru) | 2013-11-08 | 2017-10-03 | Шлюмбергер Текнолоджи Б.В. | Определение режима течения для адаптации модели потока |
WO2015070008A1 (en) | 2013-11-08 | 2015-05-14 | Schlumberger Canada Limited | Spectral analysis with spectrum deconvolution |
CN104515781B (zh) * | 2014-12-31 | 2018-01-23 | 清华大学 | X射线检查系统 |
CN104516010B (zh) * | 2014-12-31 | 2018-12-11 | 清华大学 | X射线束流强度监控装置和x射线检查系统 |
CN104535592A (zh) * | 2014-12-31 | 2015-04-22 | 同方威视技术股份有限公司 | 一种滤波装置和方法及一种物质探测装置和方法 |
RU2695099C2 (ru) | 2014-12-31 | 2019-07-19 | Тсинхуа Юниверсити | Устройство контроля интенсивности рентгеновского пучка и система контроля рентгеновского излучения |
CN110487992A (zh) * | 2018-05-14 | 2019-11-22 | 中国石油天然气股份有限公司 | 原油含水率的修正装置和方法 |
CN108535284B (zh) * | 2018-06-26 | 2023-08-22 | 西安石油大学 | 一种基于射频法的双天线结构原油含水率测量传感器 |
CN108843315B (zh) * | 2018-06-28 | 2022-05-10 | 梁安民 | 一种传感式综合自动计量装置及油液质量的计算方法 |
CN110887859A (zh) * | 2018-09-07 | 2020-03-17 | 中石化石油工程技术服务有限公司 | 钻井液含气量和含油量一体化在线检测仪器 |
NO20181382A1 (en) * | 2018-10-26 | 2020-04-27 | Roxar Flow Measurement As | Flow measuring system |
US10890544B1 (en) * | 2019-12-18 | 2021-01-12 | Field Service Solutions LLC | Nuclear densitometer assemblies for hydraulic fracturing |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4663711A (en) * | 1984-06-22 | 1987-05-05 | Shell Oil Company | Method of analyzing fluid saturation using computerized axial tomography |
CN1086602A (zh) * | 1993-06-16 | 1994-05-11 | 中国科学院近代物理研究所 | 原油含气、含水率自动测量仪 |
CN2629032Y (zh) * | 2003-08-08 | 2004-07-28 | 兰州科庆仪器仪表有限责任公司 | 原油低含水率在线监测仪 |
Family Cites Families (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2088050A (en) * | 1980-11-25 | 1982-06-03 | Kendall Ernest John Michael | Gamma Ray Analysis of Multi- component Material |
US4618975A (en) * | 1984-12-21 | 1986-10-21 | At&T Technologies, Inc. | Method and apparatus for analyzing a porous nonhomogeneous cylindrical object |
CA1257712A (en) * | 1985-11-27 | 1989-07-18 | Toshimasa Tomoda | Metering choke |
CN86105543A (zh) | 1986-07-24 | 1988-02-24 | 中国科学院近代物理研究所 | 原油含水率自动监测仪 |
US5014288A (en) * | 1989-04-20 | 1991-05-07 | Measurex Corporation | X-ray coating weight controller and sensor |
MY123677A (en) * | 1993-04-26 | 2006-05-31 | Shell Int Research | Fluid composition meter |
CN1112677A (zh) | 1994-05-25 | 1995-11-29 | 刘瑞复 | 阻抗法原油含水率高精度测试仪 |
AU719913B2 (en) * | 1996-05-02 | 2000-05-18 | Shell Internationale Research Maatschappij B.V. | Method and meter for measuring the composition of a multiphase fluid |
CN1186236A (zh) | 1996-12-27 | 1998-07-01 | 石油大学(华东) | 电容传感器非均匀探头 |
US6097786A (en) * | 1998-05-18 | 2000-08-01 | Schlumberger Technology Corporation | Method and apparatus for measuring multiphase flows |
CN2349574Y (zh) | 1998-10-15 | 1999-11-17 | 长春市宏宇电子节能设备开发有限责任公司 | 一种原油含水率测定仪双探头传感器 |
CN2359692Y (zh) | 1998-12-02 | 2000-01-19 | 大庆石油管理局第五采油厂 | 原油低含水率自动监测仪 |
CN2383068Y (zh) | 1999-03-16 | 2000-06-14 | 徐谦 | 原油含水监测仪 |
NZ518315A (en) * | 1999-10-21 | 2003-10-31 | Foss Electric As | Dual X-ray absorptiometry, method and apparatus for determination of properties of food or feed |
CN1129005C (zh) | 1999-12-29 | 2003-11-26 | 西安交通大学 | 测量原油含水率的装置及其方法 |
CN2452022Y (zh) | 2000-11-29 | 2001-10-03 | 长春市宏宇电子节能设备开发有限责任公司 | 一种高精度原油含水率测定仪 |
CN2646705Y (zh) | 2003-08-13 | 2004-10-06 | 长春市宏宇电子节能设备开发有限责任公司 | 一种有纵向擦除器的原油含水率测定仪 |
CN2646704Y (zh) | 2003-08-13 | 2004-10-06 | 长春市宏宇电子节能设备开发有限责任公司 | 一种有横向旋转擦除器的原油含水率测定仪 |
NO321752B1 (no) | 2004-06-07 | 2006-06-26 | Roxar Flow Measurement As | System for sanddeteksjon ved innsnevringer eller stromhindringer i ror |
CN1789969A (zh) | 2005-12-02 | 2006-06-21 | 仇小琦 | 测量储油罐中原油含水率和纯原油含量的装置 |
JP4725350B2 (ja) | 2006-02-20 | 2011-07-13 | 株式会社島津製作所 | 透過x線測定方法 |
JP4614001B2 (ja) * | 2006-04-13 | 2011-01-19 | 株式会社島津製作所 | 透過x線を用いた三次元定量方法 |
US7507952B2 (en) * | 2006-05-25 | 2009-03-24 | Schlumberger Technology Corporation | Apparatus and method for fluid density determination |
US7684540B2 (en) * | 2006-06-20 | 2010-03-23 | Schlumberger Technology Corporation | Apparatus and method for fluid phase fraction determination using x-rays |
RU2307328C1 (ru) * | 2006-06-23 | 2007-09-27 | Открытое акционерное общество "Техприбор" | Система измерения покомпонентного расхода трехкомпонентного газожидкостного потока нефтяных скважин |
CN101261235B (zh) | 2008-05-06 | 2010-12-08 | 罗平安 | 原油中含气率和含水率的双能χ射线测量方法 |
-
2008
- 2008-05-06 CN CN2008100972031A patent/CN101261235B/zh not_active Expired - Fee Related
-
2009
- 2009-01-21 US US12/991,170 patent/US8718230B2/en not_active Expired - Fee Related
- 2009-01-21 WO PCT/CN2009/070253 patent/WO2009135390A1/zh active Application Filing
- 2009-01-21 EP EP09741665.5A patent/EP2275804B1/en not_active Not-in-force
- 2009-01-21 RU RU2010146142/28A patent/RU2466383C2/ru not_active IP Right Cessation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4663711A (en) * | 1984-06-22 | 1987-05-05 | Shell Oil Company | Method of analyzing fluid saturation using computerized axial tomography |
CN1086602A (zh) * | 1993-06-16 | 1994-05-11 | 中国科学院近代物理研究所 | 原油含气、含水率自动测量仪 |
CN2629032Y (zh) * | 2003-08-08 | 2004-07-28 | 兰州科庆仪器仪表有限责任公司 | 原油低含水率在线监测仪 |
Non-Patent Citations (6)
Title |
---|
任晓峰等.射线型原油含气含水率分析仪的原理及应用.自动化仪表28 10.2007,28(10),11-13,17. |
任晓峰等.射线型原油含气含水率分析仪的原理及应用.自动化仪表28 10.2007,28(10),11-13,17. * |
田晶京.γ射线原油含水率检测系统的研究.中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑 6.2006,(6),12-24. |
田晶京.γ射线原油含水率检测系统的研究.中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑 6.2006,(6),12-24. * |
白秋果等.双能γ射线透射法原油含水含气的分析与计量系统.核电子学与探测技术20 3.2000,20(3),187-190. |
白秋果等.双能γ射线透射法原油含水含气的分析与计量系统.核电子学与探测技术20 3.2000,20(3),187-190. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107085089A (zh) * | 2017-05-21 | 2017-08-22 | 魏建军 | 一种原油含水率实时测量系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101261235A (zh) | 2008-09-10 |
EP2275804A4 (en) | 2012-02-01 |
US8718230B2 (en) | 2014-05-06 |
WO2009135390A1 (zh) | 2009-11-12 |
EP2275804B1 (en) | 2018-03-07 |
US20120020458A1 (en) | 2012-01-26 |
EP2275804A1 (en) | 2011-01-19 |
RU2466383C2 (ru) | 2012-11-10 |
RU2010146142A (ru) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101261235B (zh) | 原油中含气率和含水率的双能χ射线测量方法 | |
CN101261236A (zh) | 原油中含气率和含水率的双能γ射线测量方法 | |
Roshani et al. | Intelligent densitometry of petroleum products in stratified regime of two phase flows using gamma ray and neural network | |
El Abd | Intercomparison of gamma ray scattering and transmission techniques for gas volume fraction measurements in two phase pipe flow | |
CN101183065B (zh) | 利用伽玛背散射的密度测量 | |
CN105890689B (zh) | 一种测量湿气中气油水三相质量流量的测量装置及测量方法 | |
US2953681A (en) | System for measuring mass flow rate by radiation | |
CN103852475A (zh) | 一种基于伽马射线的多道测钾仪 | |
Salgado et al. | Monitoring system of oil by-products interface in pipelines using the gamma radiation attenuation | |
Abdul-Majid | Determination of wax deposition and corrosion in pipelines by neutron back diffusion collimation and neutron capture gamma rays | |
CN103399025A (zh) | 一种在线测量多相流中质量含砂率的方法和装置 | |
CN102879410B (zh) | 一种在线测量油水气多相流含率的方法及装置 | |
CN115266658A (zh) | 一种湿气油气水三相流体的含液率测量方法及装置 | |
CN201335814Y (zh) | 油气田中多相流指标的真双能χ射线测量装置 | |
US2957986A (en) | Means of studying oil consumption in an engine | |
US3508047A (en) | Method and apparatus for the accurate analysis of hydrocarbon material using a multiple of radiation sources | |
CN101852746A (zh) | 多相流中含水率与含气率的x射线测量装置及其测量方法 | |
US12025481B2 (en) | Device for measuring mass flow rate of multiphase flow based on ray coincidence measurement | |
CN103076057B (zh) | 一种多相流流量计 | |
CN113945248A (zh) | 一种四相混相质量流量的在线计量方法及装置 | |
Osman et al. | Measurement of void fraction in pipes by nuclear transmission based techniques | |
CN201352205Y (zh) | 油气田中多相流指标的伪双能χ射线测量装置 | |
CN205785377U (zh) | 一种测量湿气中气油水三相质量流量的测量装置 | |
Kopteva et al. | Noncontact measurement method for assessing oil quality and quantity in main pipeline | |
CN201416444Y (zh) | 油气田中多相流指标的双能χ射线测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101208 Termination date: 20210506 |
|
CF01 | Termination of patent right due to non-payment of annual fee |