CN112129707A - 油田集输管道内采出液ph值的计算方法 - Google Patents

油田集输管道内采出液ph值的计算方法 Download PDF

Info

Publication number
CN112129707A
CN112129707A CN201910549378.XA CN201910549378A CN112129707A CN 112129707 A CN112129707 A CN 112129707A CN 201910549378 A CN201910549378 A CN 201910549378A CN 112129707 A CN112129707 A CN 112129707A
Authority
CN
China
Prior art keywords
concentration
ion concentration
ion
hydrogen
reaction
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.)
Granted
Application number
CN201910549378.XA
Other languages
English (en)
Other versions
CN112129707B (zh
Inventor
唐德志
谷坛
陈宏健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201910549378.XA priority Critical patent/CN112129707B/zh
Publication of CN112129707A publication Critical patent/CN112129707A/zh
Application granted granted Critical
Publication of CN112129707B publication Critical patent/CN112129707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/75Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
    • G01N21/77Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
    • G01N21/78Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator producing a change of colour
    • G01N21/80Indicating pH value
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/26Oils; Viscous liquids; Paints; Inks
    • G01N33/28Oils, i.e. hydrocarbon liquids
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computing Systems (AREA)
  • General Chemical & Material Sciences (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Plasma & Fusion (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

本发明公开了一种油田集输管道内采出液PH值的计算方法,属于油田输送技术领域。本发明通过获取油田集输管线内采出液的Ca2+、Mg2+、HCO3 以及SO4 2‑的浓度,结合油田集输管线内气体浓度及分压,并采用离子逐步添加的方法,对油田集输管道内采出液的原位pH值进行计算。这不仅考虑了油田集输管道内的温度以及各个气相的分压,而且也基于溶液正负电荷平衡及物质的量守恒定律,采用离子逐步添加的方式,分步计算各步骤的离子浓度及平衡情况,充分考虑了溶质添加过程中的离子沉积反应,提高了油田集输管道原位pH值的计算精度。

Description

油田集输管道内采出液PH值的计算方法
技术领域
本发明涉及油田输送技术领域,特别涉及一种油田集输管道内采出液PH值的计算方法。
背景技术
在油田集输管道输送采出液的过程中,由于采出液中含有大量的酸性介质(例如二氧化碳、硫化氢),这会腐蚀油田集输管道,严重时会造成油气的泄漏。其中,油田集输管道发生腐蚀的程度与采出液的PH值密切相关。
目前,通常先对油田集输管道内的采出液进行降温、降压和脱气处理,然后再测量采出液的PH值。
发明人发现现有技术至少存在以下问题:
现有技术所提供的方法,没有考虑油田集输管道的工况条件,即带压带温带气,无法精确确定油田集输管道内采出液的PH值,也就难以准确反映油田集输管道内采出液的原位腐蚀性信息,进而制约集输管道内腐蚀失效分析以及科学有效的防控措施。
发明内容
本发明实施例提供了一种油田集输管道内采出液PH值的计算方法,可以解决上述问题。所述技术方案如下:
一种油田集输管道内采出液PH值的计算方法,所述计算方法包括:
获取油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量以及所述油田集输管道内采出液中每种离子浓度;
基于所述油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量,确定所述采出液的碳酸浓度、硫化氢浓度;
基于所述油田集输管道内的温度以及所述采出液的每种离子浓度、碳酸浓度和硫化氢浓度,确定第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度,其中,所述第一体系为氯化钠和氯化钾所组成的混合溶液;
向所述第一体系中加入碳酸氢根离子,形成第二体系,然后基于所述油田集输管道内的温度、所述第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度以及所述采出液的碳酸氢根离子、碳酸浓度,确定所述第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子浓度;
向所述第二体系中加入硫酸根离子,形成第三体系,然后基于所述油田集输管道内的温度、所述第二体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的硫酸氢根离子、硫化氢浓度、碳酸浓度,确定所述第三体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度;
向所述第三体系中加入钙离子和镁离子,形成第四体系,然后基于所述油田集输管道内的温度、所述第三体系的氢离子浓度以及所述采出液的碳酸浓度,确定所述第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度;
将所述钙离子沉积反应的临界浓度与所述采出液的钙离子浓度以及所述镁离子沉积反应的临界浓度值与所述采出液的镁离子浓度进行依次比较,并基于比较结果,确定所述第四体系的氢离子浓度;
根据所述第四体系的氢离子浓度,确定所述油田集输管道内采出液的PH值。
在一种可能的设计中,所述第四体系的氢离子浓度通过如下方法确定:
若所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度,将所述第四体系记作4a1体系,然后基于所述第三体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的钙离子浓度、碳酸浓度,确定所述4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度,
若所述述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度,将所述第四体系记作4b1体系,然后基于所述第三体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的钙离子浓度、硫化氢浓度、碳酸浓度,确定所述4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度;
若所述镁离子沉积反应的临界浓度值小于或等于所述采出液的镁离子浓度,将所述第四体系记作4a2体系,且当所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度时,基于所述4a1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4a2体系的氢离子浓度,或者,当所述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度时,基于所述4b1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4a2体系的氢离子浓度,
若所述镁离子沉积反应的临界浓度值大于所述采出液的镁离子浓度,将所述第四体系记作4b2体系,且当所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度时,基于所述4a1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4b2体系的氢离子浓度,或者,当所述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度时,基于所述4b1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4b2体系的氢离子浓度;
若所述镁离子沉积反应的临界浓度值小于或等于所述采出液的镁离子浓度,所述4a2体系的氢离子浓度看作所述第四系统的氢离子浓度,若所述镁离子沉积反应的临界浓度值大于所述采出液的镁离子浓度,所述4b2体系的氢离子浓度看作所述第四系统的氢离子浓度。
在一种可能的设计中,所述采出液的碳酸浓度、硫化氢浓度通过公式(1)至公式(8)计算得到:
Figure BDA0002105050840000031
Figure BDA0002105050840000032
Figure BDA0002105050840000033
Figure BDA0002105050840000034
Figure BDA0002105050840000035
Figure BDA0002105050840000036
Figure BDA0002105050840000037
Figure BDA0002105050840000038
式中:pCO2—所述油田集输管道内的CO2分压值,P—所述油田集输管道内的气体总压,xCO2—所述油田集输管道内的CO2百分含量,pH2S—所述油田集输管道内的H2S分压值,xH2S—所述油田集输管道内的H2S百分含量,
Figure BDA0002105050840000041
—所述采出液的H2CO3浓度,Khy—反应
Figure BDA0002105050840000042
的平衡常数,取值为2.58×10-3,Ksol—反应
Figure BDA0002105050840000043
的平衡常数,
Figure BDA0002105050840000044
—气体逸度系数,T—所述油田集输管道内的绝对温度,
Figure BDA0002105050840000045
—所述采出液的H2S浓度,
Figure BDA0002105050840000046
—反应
Figure BDA0002105050840000047
的平衡常数,Tf—所述油田集输管道内的华氏温度,I—所述采出液的离子强度,mi—所述采出液中第i种离子的浓度,zi—所述采出液中第i种离子所对应的离子价态。
在一种可能的设计中,所述第一体系的氢离子浓度、碳酸氢根离子浓度和碳酸根离子浓度通过公式(9)至公式(16)计算得到:
Figure BDA0002105050840000048
Figure BDA0002105050840000049
Figure BDA00021050508400000410
Figure BDA00021050508400000411
Figure BDA00021050508400000412
Figure BDA00021050508400000413
Figure BDA00021050508400000414
Figure BDA00021050508400000415
式中:
Figure BDA00021050508400000416
—所述第一体系的H+浓度,Kwa—反应
Figure BDA00021050508400000417
的平衡常数,Kca—反应
Figure BDA00021050508400000418
的平衡常数,Ka,1—反应
Figure BDA00021050508400000419
的平衡常数,Ka,2—反应
Figure BDA00021050508400000420
的平衡常数,Kbi—反应
Figure BDA00021050508400000421
的平衡常数,
Figure BDA00021050508400000422
—所述第一体系的HCO3 -浓度,
Figure BDA00021050508400000423
—所述第一体系的CO3 2-浓度。
在一种可能的设计中,所述第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子通过公式(17)至(21)计算得到:
Figure BDA00021050508400000424
Figure BDA00021050508400000425
Figure BDA00021050508400000426
Figure BDA0002105050840000051
Figure BDA0002105050840000052
式中,
Figure BDA0002105050840000053
—所述第二体系的H+浓度,x—加入HCO3 -后通过反应
Figure BDA0002105050840000054
消耗的HCO3 -浓度,y—加入HCO3 -后通过反应
Figure BDA0002105050840000055
Figure BDA0002105050840000056
消耗的HCO3 -浓度,
Figure BDA0002105050840000057
—所述采出液的HCO3 -浓度,
Figure BDA0002105050840000058
—所述第二体系的HCO3 -浓度,
Figure BDA0002105050840000059
—所述第二体系的HS-浓度。
在一种可能的设计中,所述第三体系的氢离子浓度、碳酸氢根离子、碳酸根离子浓度、硫化氢根离子浓度通过公式(22)至(29)计算得到:
Figure BDA00021050508400000510
Figure BDA00021050508400000511
Figure BDA00021050508400000512
Figure BDA00021050508400000513
Figure BDA00021050508400000514
Figure BDA00021050508400000515
Figure BDA00021050508400000516
Figure BDA00021050508400000517
式中,
Figure BDA00021050508400000518
—所述第三体系的H+浓度,d—反应
Figure BDA00021050508400000519
消耗的SO4 2-浓度,e—反应
Figure BDA00021050508400000520
生成的H2S浓度,f—反应
Figure BDA00021050508400000521
Figure BDA00021050508400000522
生成的H2CO3浓度,
Figure BDA00021050508400000523
—所述采出液的SO4 2-浓度,
Figure BDA00021050508400000524
—反应
Figure BDA00021050508400000525
的平衡常数,
Figure BDA00021050508400000526
—所述第三体系的HCO3 -浓度,
Figure BDA00021050508400000527
—所述第三体系的CO3 2-浓度,
Figure BDA00021050508400000528
—所述第三体系的HS-浓度。
在一种可能的设计中,所述第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度通过公式(30)至公式(33)计算得到:
Figure BDA00021050508400000529
Figure BDA00021050508400000530
Figure BDA00021050508400000531
Figure BDA00021050508400000532
式中,
Figure BDA0002105050840000061
—所述第四体系的钙离子沉积反应的临界浓度值,
Figure BDA0002105050840000062
—所述第四体系的镁离子沉积反应的临界浓度,Ksp1—反应
Figure BDA0002105050840000063
的平衡常数,Ksp2—反应
Figure BDA0002105050840000064
的平衡常数。
在一种可能的设计中,所述4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(34)至公式(40)计算得到:
Figure BDA0002105050840000065
Figure BDA0002105050840000066
Figure BDA0002105050840000067
Figure BDA0002105050840000068
Figure BDA0002105050840000069
Figure BDA00021050508400000610
Figure BDA00021050508400000611
式中:
Figure BDA00021050508400000612
—所述4a1体系的氢离子浓度,i—通过反应
Figure BDA00021050508400000613
生成的CO3 2-浓度,j—通过反应
Figure BDA00021050508400000614
生成的HCO3 -浓度,h—通过反应
Figure BDA00021050508400000615
消耗的Ca2+浓度,
Figure BDA00021050508400000616
—所述4a1体系的HCO3 -浓度,
Figure BDA00021050508400000617
—所述4a1体系的CO3 2-浓度,
Figure BDA00021050508400000618
—所述4a1体系的HS-浓度;
所述4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(41)至公式(48)计算得到:
Figure BDA00021050508400000619
Figure BDA00021050508400000620
Figure BDA00021050508400000621
Figure BDA00021050508400000622
Figure BDA00021050508400000623
Figure BDA00021050508400000624
Figure BDA00021050508400000625
Figure BDA00021050508400000626
式中:
Figure BDA00021050508400000627
—所述4b1体系的氢离子浓度,k—通过反应
Figure BDA00021050508400000628
Figure BDA0002105050840000071
消耗的Ca2+浓度,l—通过反应
Figure BDA0002105050840000072
生成的H2CO3浓度,m—通过反应
Figure BDA0002105050840000073
生成的H2S浓度,Ksp5—反应
Figure BDA0002105050840000074
的平衡常数,
Figure BDA0002105050840000075
—所述4b1体系的HCO3 -浓度,
Figure BDA0002105050840000076
—所述4b1体系的CO3 2-浓度,
Figure BDA0002105050840000077
—所述4b1体系的HS-浓度。
在一种可能的设计中,所述体系4a2的氢离子浓度通过公式(49)至公式(52)计算得到:
Figure BDA0002105050840000078
Figure BDA0002105050840000079
Figure BDA00021050508400000710
Figure BDA00021050508400000711
式中:
Figure BDA00021050508400000712
—所述体系4a2的氢离子浓度,
n=a或b,若
Figure BDA00021050508400000713
则n=a,否则,n=b,
s—通过HCO3 -电离生成的CO3 2-浓度,
t—通过H2CO3电离生成的HCO3 -浓度,
r—通过反应
Figure BDA00021050508400000714
消耗的Mg2+浓度;
所述体系4b2的氢离子浓度通过公式(53)至公式(57)计算得到:
Figure BDA00021050508400000715
Figure BDA00021050508400000716
Figure BDA00021050508400000723
Figure BDA00021050508400000717
Figure BDA00021050508400000718
式中:
Figure BDA00021050508400000719
—所述体系4b2的氢离子浓度,
n=a或b,若
Figure BDA00021050508400000720
则n=a,否则,n=b,
u—通过反应
Figure BDA00021050508400000721
消耗的Mg2+浓度,
v—通过反应
Figure BDA00021050508400000722
生成的H2CO3浓度,
w—通过反应
Figure BDA0002105050840000081
生成的H2S浓度,
Ksp6—反应
Figure BDA0002105050840000082
的平衡常数。
在一种可能的设计中,所述油田集输管道内采出液的pH值通过公式(58)计算得到:
Figure BDA0002105050840000083
式中:
n=a或b,若
Figure BDA0002105050840000084
则n=a;否则,n=b。
本发明实施例提供的技术方案带来的有益效果至少包括:
本发明实施提供的油田集输管道内采出液PH值的计算方法,通过获取油田集输管线内采出液的Ca2+、Mg2+、HCO3 -以及SO4 2-的浓度,结合油田集输管线内气体浓度及分压,并采用离子逐步添加的方法,对油田集输管道内采出液的原位pH值进行计算。这不仅考虑了油田集输管道内的温度以及各个气相的分压,而且也基于溶液正负电荷平衡及物质的量守恒定律,采用离子逐步添加的方式,分步计算各步骤的离子浓度及平衡情况,充分考虑了溶质添加过程中的离子沉积反应,提高了油田集输管道原位pH值的计算精度。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种油田集输管道内采出液PH值的计算方法,该方法包括如下步骤:
步骤S1、获取油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量以及所述油田集输管道内采出液中每种离子浓度。
步骤S2、基于油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量,确定采出液的碳酸浓度、硫化氢浓度。
步骤S3、基于油田集输管道内的温度以及所述采出液的每种离子浓度、碳酸浓度和硫化氢浓度,确定第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度,其中,第一体系为氯化钠和氯化钾所组成的混合溶液。
步骤S4、向第一体系中加入碳酸氢根离子,形成第二体系,然后基于油田集输管道内的温度、第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度以及采出液的碳酸氢根离子、碳酸浓度,确定第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子浓度。
步骤S5、向第二体系中加入硫酸根离子,形成第三体系,然后基于油田集输管道内的温度、第二体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及采出液的硫酸氢根离子、硫化氢浓度、碳酸浓度,确定第三体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度。
步骤S6、向第三体系中加入钙离子和镁离子,形成第四体系,然后基于油田集输管道内的温度、第三体系的氢离子浓度以及采出液的碳酸浓度,确定第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度。
步骤S7、将钙离子沉积反应的临界浓度与采出液的钙离子浓度以及镁离子沉积反应的临界浓度值与采出液的镁离子浓度进行依次比较,并基于比较结果,确定第四体系的氢离子浓度。
步骤S8、根据第四体系的氢离子浓度,确定油田集输管道内采出液的PH值。
本发明实施提供的油田集输管道内采出液PH值的计算方法,通过获取油田集输管线内采出液的Ca2+、Mg2+、HCO3 -以及SO4 2-的浓度,结合油田集输管线内气体浓度及分压,并采用离子逐步添加的方法,对油田集输管道内采出液的原位pH值进行计算。这不仅考虑了油田集输管道内的温度以及各个气相的分压,而且也基于溶液正负电荷平衡及物质的量守恒定律,采用离子逐步添加的方式,分步计算各步骤的离子浓度及平衡情况,充分考虑了溶质添加过程中的离子沉积反应,提高了油田集输管道原位pH值的计算精度。
下面就本发明实施例提供的计算方法的每个步骤进行描述:
针对于步骤S1而言,可利用光谱分析法测试油田集输管道内采出液离子的每种离子浓度,例如Ca2+浓度
Figure BDA0002105050840000091
Mg2+浓度
Figure BDA0002105050840000092
HCO3 -浓度
Figure BDA0002105050840000093
SO4 2-浓度
Figure BDA0002105050840000094
可利用色长度检测管法测定油田集输管道内二氧化碳气体百分含量(xCO2)和硫化氢气体百分含量(xH2S)。
可利用温度检测枪测定油田集输管道内的温度。
针对于步骤S2而言,利用气液平衡原理,可通过公式(1)至公式(8)计算得到采出液的碳酸浓度、硫化氢浓度:
Figure BDA0002105050840000101
Figure BDA0002105050840000102
Figure BDA0002105050840000103
Figure BDA0002105050840000104
Figure BDA0002105050840000105
Figure BDA0002105050840000106
Figure BDA0002105050840000107
Figure BDA0002105050840000108
式中:pCO2—油田集输管道内的CO2分压值,P—油田集输管道内的气体总压,xCO2—油田集输管道内的CO2百分含量,pH2S—油田集输管道内的H2S分压值,xH2S—油田集输管道内的H2S百分含量,
Figure BDA0002105050840000109
—采出液的H2CO3浓度,Khy—反应
Figure BDA00021050508400001010
的平衡常数,取值为2.58×10-3,Ksol—反应
Figure BDA00021050508400001011
的平衡常数,
Figure BDA00021050508400001012
—气体逸度系数,T—油田集输管道内的绝对温度,
Figure BDA00021050508400001013
—采出液的H2S浓度,
Figure BDA00021050508400001014
—反应
Figure BDA00021050508400001015
的平衡常数,Tf—油田集输管道内的华氏温度,I—采出液的离子强度,mi—采出液中第i种离子的浓度,zi—采出液中第i种离子所对应的离子价态。
其中,T的单位为K,Tf的单位为华氏温度,这两者之间的关系为Tf=9×(T-273.15)/5+32。
针对步骤S3至步骤S7而言,是采用离子逐步添加的方法来分步计算油田集输管道内采出液的原位pH值。
具体地,针对于步骤S3而言,第一体系仅存在NaCl和KCl,此时溶液中不存在沉积过程。本发明实施例根据溶液中阴阳离子电荷平衡及溶液中H2CO3(液体)和H2S(液体)的电离平衡可以得到关于第一体系氢离子浓度的一元三次方程(即公式9),求解该一元三次方程,可以得到第一体系的H+浓度,并结合阴阳离子电荷平衡可以获得HCO3 -、CO3 2-浓度。其中,第一体系的氢离子浓度、碳酸氢根离子浓度和碳酸根离子浓度通过公式(9)至公式(16)计算得到:
Figure BDA00021050508400001016
Figure BDA00021050508400001017
Figure BDA0002105050840000111
Figure BDA0002105050840000112
Figure BDA0002105050840000113
Figure BDA0002105050840000114
Figure BDA0002105050840000115
Figure BDA0002105050840000116
式中:
Figure BDA0002105050840000117
—第一体系的H+浓度,Kwa—反应
Figure BDA0002105050840000118
的平衡常数,Kca—反应
Figure BDA0002105050840000119
的平衡常数,Ka,1—反应
Figure BDA00021050508400001110
的平衡常数,Ka,2—反应
Figure BDA00021050508400001111
的平衡常数,Kbi—反应
Figure BDA00021050508400001112
的平衡常数,
Figure BDA00021050508400001113
—第一体系的HCO3 -浓度,
Figure BDA00021050508400001114
—第一体系的CO3 2-浓度。
具体地,针对于步骤S4而言,向第一体系中加入碳酸氢根离子,形成第二体系。本发明实施例根据相应的水化学平衡反应及物质的量守恒计算出新平衡下的H+、HCO3 -以及HS-离子浓度。其中,第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子通过公式(17)至(21)计算得到:
Figure BDA00021050508400001115
Figure BDA00021050508400001116
Figure BDA00021050508400001117
Figure BDA00021050508400001118
Figure BDA00021050508400001119
式中,
Figure BDA00021050508400001120
—第二体系的H+浓度,x—加入HCO3 -后通过反应
Figure BDA00021050508400001121
Figure BDA00021050508400001122
消耗的HCO3 -浓度,y—加入HCO3 -后通过反应
Figure BDA00021050508400001123
消耗的HCO3 -浓度,
Figure BDA00021050508400001124
—采出液的HCO3 -浓度,
Figure BDA00021050508400001125
—第二体系的HCO3 -浓度,
Figure BDA00021050508400001126
—第二体系的HS-浓度。
具体地,针对于步骤S5而言,向第二体系中加入SO4 2-离子,形成第三体系。此时由于SO4 2-离子与H+结合会导致溶液的pH升高,同时由于H+主要来源于H2S和H2CO3的电离过程,所以主要考虑这三个反应的平衡,由此可得到第三体系在平衡条件下H+、HCO3 -、CO3 2-以及HS-的离子浓度。其中,第三体系的氢离子浓度、碳酸氢根离子、碳酸根离子浓度、硫化氢根离子浓度通过公式(22)至(29)计算得到:
Figure BDA0002105050840000121
Figure BDA0002105050840000122
Figure BDA0002105050840000123
Figure BDA0002105050840000124
Figure BDA0002105050840000125
Figure BDA0002105050840000126
Figure BDA0002105050840000127
Figure BDA0002105050840000128
式中,
Figure BDA0002105050840000129
—第三体系的H+浓度,d—反应
Figure BDA00021050508400001210
消耗的SO4 2-浓度,e—反应
Figure BDA00021050508400001211
生成的H2S浓度,f—反应
Figure BDA00021050508400001212
生成的H2CO3浓度,
Figure BDA00021050508400001213
—采出液的SO4 2-浓度,
Figure BDA00021050508400001214
—反应
Figure BDA00021050508400001215
的平衡常数,
Figure BDA00021050508400001216
—第三体系的HCO3 -浓度,
Figure BDA00021050508400001217
—第三体系的CO3 2-浓度,
Figure BDA00021050508400001218
—第三体系的HS-浓度。
具体地,针对于步骤S6而言,向第三体系中加入Ca2+和Mg2+离子,形成第四体系。此时,利用第三体系的H+浓度计算得到Ca2+和Mg2+沉积反应的临界浓度值。其中,第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度通过公式(30)至公式(33)计算得到:
Figure BDA00021050508400001219
Figure BDA00021050508400001220
Figure BDA00021050508400001221
Figure BDA00021050508400001222
式中,
Figure BDA00021050508400001223
—第四体系的钙离子沉积反应的临界浓度值,
Figure BDA00021050508400001224
—第四体系的镁离子沉积反应的临界浓度,Ksp1—反应
Figure BDA00021050508400001225
的平衡常数,Ksp2—反应
Figure BDA00021050508400001226
的平衡常数。
具体地,针对于步骤S7而言,将第四体系的钙离子沉积反应的临界浓度值与采出液的钙离子浓度进行比较,以及将第四体系的镁离子沉积反应的临界浓度值与采出液的镁离子浓度进行比较,若
Figure BDA00021050508400001227
(或
Figure BDA00021050508400001228
),则会形成钙镁的沉积物;否则钙镁离子浓度为其水解平衡浓度;之后根据二者数值大小确定Ca2+和Mg2+在溶液中的平衡反应,进而确定第四体系的H+浓度。
其中,第四体系的氢离子浓度通过如下方法确定:
(1)先考虑加入Ca2+后溶液离子浓度发生的变化:
(1.1)若钙离子沉积反应的临界浓度值小于或等于采出液的钙离子浓度,将第四体系记作4a1体系,然后基于第三体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及采出液的钙离子浓度、碳酸浓度,确定4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度。
作为一种示例,4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(34)至公式(40)计算得到:
Figure BDA0002105050840000131
Figure BDA0002105050840000132
Figure BDA0002105050840000133
Figure BDA0002105050840000134
Figure BDA0002105050840000135
Figure BDA0002105050840000136
Figure BDA0002105050840000137
式中:
Figure BDA0002105050840000138
—4a1体系的氢离子浓度,i—通过反应
Figure BDA0002105050840000139
生成的CO3 2-浓度,j—通过反应
Figure BDA00021050508400001310
生成的HCO3 -浓度,h—通过反应
Figure BDA00021050508400001311
消耗的Ca2+浓度,
Figure BDA00021050508400001312
—4a1体系的HCO3 -浓度,
Figure BDA00021050508400001313
—4a1体系的CO3 2-浓度,
Figure BDA00021050508400001314
—4a1体系的HS-浓度。
(1.2)若述钙离子沉积反应的临界浓度值大于采出液的钙离子浓度,将第四体系记作4b1体系,然后基于第三体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及采出液的钙离子浓度、硫化氢浓度、碳酸浓度,确定4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度。
作为一种示例,4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(41)至公式(48)计算得到:
Figure BDA00021050508400001315
Figure BDA00021050508400001316
Figure BDA0002105050840000141
Figure BDA0002105050840000142
Figure BDA0002105050840000143
Figure BDA0002105050840000144
Figure BDA0002105050840000145
Figure BDA0002105050840000146
式中:
Figure BDA0002105050840000147
—4b1体系的氢离子浓度,k—通过反应
Figure BDA0002105050840000148
消耗的Ca2+浓度,l—通过反应
Figure BDA0002105050840000149
生成的H2CO3浓度,m—通过反应
Figure BDA00021050508400001410
生成的H2S浓度,Ksp5—反应
Figure BDA00021050508400001411
的平衡常数,
Figure BDA00021050508400001412
—4b1体系的HCO3 -浓度,
Figure BDA00021050508400001413
—4b1体系的CO3 2-浓度,
Figure BDA00021050508400001414
—4b1体系的HS-浓度。
(2)再考虑加入Mg2+后溶液离子浓度发生的变化:
(2.1)若镁离子沉积反应的临界浓度值小于或等于采出液的镁离子浓度,将第四体系记作4a2体系,且当钙离子沉积反应的临界浓度值小于或等于采出液的钙离子浓度时,基于4a1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及采出液的镁离子浓度、碳酸浓度,确定4a2体系的氢离子浓度,或者,当钙离子沉积反应的临界浓度值大于采出液的钙离子浓度时,基于4b1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及采出液的镁离子浓度、碳酸浓度,确定4a2体系的氢离子浓度。
作为一种示例,体系4a2的氢离子浓度通过公式(49)至公式(52)计算得到:
Figure BDA00021050508400001415
Figure BDA00021050508400001416
Figure BDA00021050508400001417
Figure BDA00021050508400001418
式中:
Figure BDA00021050508400001419
—体系4a2的氢离子浓度;
n=a或b,若
Figure BDA00021050508400001420
则n=a,否则,n=b;
s—通过HCO3 -电离生成的CO3 2-浓度;
t—通过H2CO3电离生成的HCO3 -浓度;
r—通过反应
Figure BDA0002105050840000151
消耗的Mg2+浓度。
(2.2)若镁离子沉积反应的临界浓度值大于采出液的镁离子浓度,将第四体系记作4b2体系,且当钙离子沉积反应的临界浓度值小于或等于采出液的钙离子浓度时,基于4a1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及采出液的镁离子浓度、碳酸浓度,确定4b2体系的氢离子浓度,或者,当钙离子沉积反应的临界浓度值大于采出液的钙离子浓度时,基于4b1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及采出液的镁离子浓度、碳酸浓度,确定4b2体系的氢离子浓度。
作为一种示例,体系4b2的氢离子浓度通过公式(53)至公式(57)计算得到:
Figure BDA0002105050840000152
Figure BDA0002105050840000153
Figure BDA0002105050840000154
Figure BDA0002105050840000155
Figure BDA0002105050840000156
式中:
Figure BDA0002105050840000157
—体系4b2的氢离子浓度;
n=a或b,若
Figure BDA0002105050840000158
则n=a,否则,n=b;
u—通过反应
Figure BDA0002105050840000159
消耗的Mg2+浓度;
v—通过反应
Figure BDA00021050508400001510
生成的H2CO3浓度;
w—通过反应
Figure BDA00021050508400001511
生成的H2S浓度;
Ksp6—反应
Figure BDA00021050508400001512
的平衡常数。
(3)若镁离子沉积反应的临界浓度值小于或等于采出液的镁离子浓度,4a2体系的氢离子浓度看作第四系统的氢离子浓度,若镁离子沉积反应的临界浓度值大于采出液的镁离子浓度,4b2体系的氢离子浓度看作第四系统的氢离子浓度。
具体地,针对于步骤S8而言,根据第四体系的氢离子浓度,确定油田集输管道内采出液的PH值。其中,
采出液的PH值通过公式(58)计算得到:
Figure BDA0002105050840000161
式中:n=a或b:若
Figure BDA0002105050840000162
则n=a;否则n=b。
下面以一个实例对本发明实施例提供的油田集输管道内采出液PH值的计算方法进行说明:
首先,从某油田集输管道取样口位置取得采出液,利用光谱分析获得其主要离子含量为:880mg/L Ca2+,1050mg/L Mg2+,600mg/L HCO3-以及3230mg/L SO4 2-
其次,利用着色长度检测管法测定管道内CO2及H2S气体含量,结合管道压力表读数,利用公式(1)和(2)计算得到管道内CO2及H2S分压值分别为:0.6MPa和20kPa;管道内温度读取为23.2℃。
其次,根据公式(5)~(16)计算得到第一体系的HCO3 -浓度、CO3 2-浓度、H+浓度分别为:0.00041mol/L、3.65×10-10mol/L、0.00042mol/L。。
其次,向第一体系中加入HCO3 -离子以形成第二体系,并利用公式(17)~(21)计算得到第二体系的HCO3 -浓度、H+浓度、HS-浓度分别为:0.0098mol/L、0.000036mol/L、4.97×10-5mol/L。
其次,向第二体系中加入SO4 2-离子以形成第三体系,并根据公式(22)~(29)计算第三体系的HCO3 -浓度、CO3 2-浓度、H+浓度、HS-浓度分别为:0.00992mol/L、1.01×10-7mol/L、0.000036mol/L、5.02×10-5mol/L;
其次,向第三体系中加入Ca2+和Mg2+以形成第四体系,并根据公式(28)~(31)计算确定
Figure BDA0002105050840000163
Figure BDA0002105050840000164
浓度分别为0.03257mol/L和0.1649mol/L。由于
Figure BDA0002105050840000165
Figure BDA0002105050840000166
判断出此时无沉淀反应发生。利用公式(41)~(48)以及(53)~(57)计算第四体系的H+浓度为:0.000036mol/L;
最后,利用公式(58)计算得到油田集输管道内采出液的原位pH值为4.44。
上述所有可选技术方案,可以采用任意结合形成本公开的可选实施例,在此不再一一赘述。
以上所述仅为本发明的说明性实施例,并不用以限制本发明的保护范围,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种油田集输管道内采出液PH值的计算方法,其特征在于,所述计算方法包括:
获取油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量以及所述油田集输管道内采出液中每种离子浓度;
基于所述油田集输管道内的温度、二氧化碳气体的百分含量、硫化氢气体的百分含量,确定所述采出液的碳酸浓度、硫化氢浓度;
基于所述油田集输管道内的温度以及所述采出液的每种离子浓度、碳酸浓度和硫化氢浓度,确定第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度,其中,所述第一体系为氯化钠和氯化钾所组成的混合溶液;
向所述第一体系中加入碳酸氢根离子,形成第二体系,然后基于所述油田集输管道内的温度、所述第一体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度以及所述采出液的碳酸氢根离子、碳酸浓度,确定所述第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子浓度;
向所述第二体系中加入硫酸根离子,形成第三体系,然后基于所述油田集输管道内的温度、所述第二体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的硫酸氢根离子、硫化氢浓度、碳酸浓度,确定所述第三体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度;
向所述第三体系中加入钙离子和镁离子,形成第四体系,然后基于所述油田集输管道内的温度、所述第三体系的氢离子浓度以及所述采出液的碳酸浓度,确定所述第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度;
将所述钙离子沉积反应的临界浓度与所述采出液的钙离子浓度以及所述镁离子沉积反应的临界浓度值与所述采出液的镁离子浓度进行依次比较,并基于比较结果,确定所述第四体系的氢离子浓度;
根据所述第四体系的氢离子浓度,确定所述油田集输管道内采出液的PH值。
2.根据权利要求1所述的计算方法,其特征在于,所述第四体系的氢离子浓度通过如下方法确定:
若所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度,将所述第四体系记作4a1体系,然后基于所述第三体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的钙离子浓度、碳酸浓度,确定所述4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度,
若所述述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度,将所述第四体系记作4b1体系,然后基于所述第三体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的钙离子浓度、硫化氢浓度、碳酸浓度,确定所述4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度和硫化氢根离子浓度;
若所述镁离子沉积反应的临界浓度值小于或等于所述采出液的镁离子浓度,将所述第四体系记作4a2体系,且当所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度时,基于所述4a1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4a2体系的氢离子浓度,或者,当所述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度时,基于所述4b1体系的氢离子浓度、碳酸根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4a2体系的氢离子浓度,
若所述镁离子沉积反应的临界浓度值大于所述采出液的镁离子浓度,将所述第四体系记作4b2体系,且当所述钙离子沉积反应的临界浓度值小于或等于所述采出液的钙离子浓度时,基于所述4a1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4b2体系的氢离子浓度,或者,当所述钙离子沉积反应的临界浓度值大于所述采出液的钙离子浓度时,基于所述4b1体系的氢离子浓度、硫化氢根离子浓度、碳酸氢根离子浓度以及所述采出液的镁离子浓度、碳酸浓度,确定所述4b2体系的氢离子浓度;
若所述镁离子沉积反应的临界浓度值小于或等于所述采出液的镁离子浓度,所述4a2体系的氢离子浓度看作所述第四系统的氢离子浓度,若所述镁离子沉积反应的临界浓度值大于所述采出液的镁离子浓度,所述4b2体系的氢离子浓度看作所述第四系统的氢离子浓度。
3.根据权利要求2所述的计算方法,其特征在于,所述采出液的碳酸浓度、硫化氢浓度通过公式(1)至公式(8)计算得到:
Figure FDA0002105050830000031
Figure FDA0002105050830000032
Figure FDA0002105050830000033
Figure FDA0002105050830000034
Figure FDA0002105050830000035
Figure FDA0002105050830000036
Figure FDA0002105050830000037
Figure FDA0002105050830000038
式中:pCO2—所述油田集输管道内的CO2分压值,P—所述油田集输管道内的气体总压,xCO2—所述油田集输管道内的CO2百分含量,pH2S—所述油田集输管道内的H2S分压值,xH2S—所述油田集输管道内的H2S百分含量,
Figure FDA0002105050830000039
—所述采出液的H2CO3浓度,Khy—反应
Figure FDA00021050508300000310
的平衡常数,取值为2.58×10-3,Ksol—反应
Figure FDA00021050508300000311
的平衡常数,
Figure FDA00021050508300000312
—气体逸度系数,T—所述油田集输管道内的绝对温度,
Figure FDA00021050508300000313
—所述采出液的H2S浓度,
Figure FDA00021050508300000314
—反应
Figure FDA00021050508300000315
的平衡常数,Tf—所述油田集输管道内的华氏温度,I—所述采出液的离子强度,mi—所述采出液中第i种离子的浓度,zi—所述采出液中第i种离子所对应的离子价态。
4.根据权利要求3所述的计算方法,其特征在于,所述第一体系的氢离子浓度、碳酸氢根离子浓度和碳酸根离子浓度通过公式(9)至公式(16)计算得到:
Figure FDA00021050508300000316
Figure FDA00021050508300000317
Figure FDA00021050508300000318
Figure FDA00021050508300000319
Figure FDA0002105050830000041
Figure FDA0002105050830000042
Figure FDA0002105050830000043
Figure FDA0002105050830000044
式中:
Figure FDA0002105050830000045
—所述第一体系的H+浓度,Kwa—反应
Figure FDA0002105050830000046
的平衡常数,Kca—反应
Figure FDA0002105050830000047
的平衡常数,Ka,1—反应
Figure FDA0002105050830000048
的平衡常数,Ka,2—反应
Figure FDA0002105050830000049
的平衡常数,Kbi—反应
Figure FDA00021050508300000410
的平衡常数,
Figure FDA00021050508300000411
—所述第一体系的HCO3 -浓度,
Figure FDA00021050508300000412
—所述第一体系的CO3 2-浓度。
5.根据权利要求4所述的计算方法,其特征在于,所述第二体系的氢离子浓度、碳酸氢根离子浓度、硫化氢根离子通过公式(17)至(21)计算得到:
Figure FDA00021050508300000413
Figure FDA00021050508300000414
Figure FDA00021050508300000415
Figure FDA00021050508300000416
Figure FDA00021050508300000417
式中,
Figure FDA00021050508300000418
—所述第二体系的H+浓度,x—加入HCO3 -后通过反应
Figure FDA00021050508300000419
消耗的HCO3 -浓度,y—加入HCO3 -后通过反应
Figure FDA00021050508300000420
Figure FDA00021050508300000421
消耗的HCO3 -浓度,
Figure FDA00021050508300000422
—所述采出液的HCO3 -浓度,
Figure FDA00021050508300000423
—所述第二体系的HCO3 -浓度,
Figure FDA00021050508300000424
—所述第二体系的HS-浓度。
6.根据权利要求5所述的计算方法,其特征在于,所述第三体系的氢离子浓度、碳酸氢根离子、碳酸根离子浓度、硫化氢根离子浓度通过公式(22)至(29)计算得到:
Figure FDA00021050508300000425
Figure FDA00021050508300000426
Figure FDA00021050508300000427
Figure FDA0002105050830000051
Figure FDA0002105050830000052
Figure FDA0002105050830000053
Figure FDA0002105050830000054
Figure FDA0002105050830000055
式中,
Figure FDA0002105050830000056
—所述第三体系的H+浓度,d—反应
Figure FDA0002105050830000057
消耗的SO4 2-浓度,e—反应
Figure FDA0002105050830000058
生成的H2S浓度,f—反应
Figure FDA0002105050830000059
Figure FDA00021050508300000510
生成的H2CO3浓度,
Figure FDA00021050508300000511
—所述采出液的SO4 2-浓度,
Figure FDA00021050508300000512
—反应
Figure FDA00021050508300000513
的平衡常数,
Figure FDA00021050508300000514
—所述第三体系的HCO3 -浓度,
Figure FDA00021050508300000515
—所述第三体系的CO3 2-浓度,
Figure FDA00021050508300000527
—所述第三体系的HS-浓度。
7.根据权利要求6所述的计算方法,其特征在于,所述第四体系的钙离子沉积反应的临界浓度值和镁离子沉积反应的临界浓度通过公式(30)至公式(33)计算得到:
Figure FDA00021050508300000516
Figure FDA00021050508300000517
Figure FDA00021050508300000518
Figure FDA00021050508300000519
式中,
Figure FDA00021050508300000520
—所述第四体系的钙离子沉积反应的临界浓度值,
Figure FDA00021050508300000521
—所述第四体系的镁离子沉积反应的临界浓度,Ksp1—反应
Figure FDA00021050508300000522
的平衡常数,Ksp2—反应
Figure FDA00021050508300000523
的平衡常数。
8.根据权利要求7所述的计算方法,其特征在于,所述4a1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(34)至公式(40)计算得到:
Figure FDA00021050508300000524
Figure FDA00021050508300000525
Figure FDA00021050508300000526
Figure FDA0002105050830000061
Figure FDA0002105050830000062
Figure FDA0002105050830000063
Figure FDA0002105050830000064
式中:
Figure FDA0002105050830000065
—所述4a1体系的氢离子浓度,i—通过反应
Figure FDA0002105050830000066
生成的CO3 2-浓度,j—通过反应
Figure FDA0002105050830000067
生成的HCO3 -浓度,h—通过反应
Figure FDA0002105050830000068
消耗的Ca2+浓度,
Figure FDA0002105050830000069
—所述4a1体系的HCO3 -浓度,
Figure FDA00021050508300000610
—所述4a1体系的CO3 2-浓度,
Figure FDA00021050508300000611
—所述4a1体系的HS-浓度;
所述4b1体系的氢离子浓度、碳酸氢根离子浓度、碳酸根离子浓度、硫化氢根离子浓度通过公式(41)至公式(48)计算得到:
Figure FDA00021050508300000612
Figure FDA00021050508300000613
Figure FDA00021050508300000614
Figure FDA00021050508300000615
Figure FDA00021050508300000616
Figure FDA00021050508300000617
Figure FDA00021050508300000618
Figure FDA00021050508300000619
式中:
Figure FDA00021050508300000620
—所述4b1体系的氢离子浓度,k—通过反应
Figure FDA00021050508300000621
Figure FDA00021050508300000622
消耗的Ca2+浓度,l—通过反应
Figure FDA00021050508300000623
生成的H2CO3浓度,m—通过反应
Figure FDA00021050508300000624
生成的H2S浓度,Ksp5—反应
Figure FDA00021050508300000625
的平衡常数,
Figure FDA00021050508300000626
—所述4b1体系的HCO3 -浓度,
Figure FDA00021050508300000627
—所述4b1体系的CO3 2-浓度,
Figure FDA00021050508300000628
—所述4b1体系的HS-浓度。
9.根据权利要求8所述的计算方法,其特征在于,所述体系4a2的氢离子浓度通过公式(49)至公式(52)计算得到:
Figure FDA00021050508300000629
Figure FDA0002105050830000071
Figure FDA0002105050830000072
Figure FDA0002105050830000073
式中:
Figure FDA0002105050830000074
—所述体系4a2的氢离子浓度,
n=a或b,若
Figure FDA0002105050830000075
则n=a,否则,n=b,
s—通过HCO3 -电离生成的CO3 2-浓度,
t—通过H2CO3电离生成的HCO3 -浓度,
r—通过反应
Figure FDA0002105050830000076
消耗的Mg2+浓度;
所述体系4b2的氢离子浓度通过公式(53)至公式(57)计算得到:
Figure FDA0002105050830000077
Figure FDA0002105050830000078
Figure FDA0002105050830000079
Figure FDA00021050508300000710
Figure FDA00021050508300000711
式中:
Figure FDA00021050508300000712
—所述体系4b2的氢离子浓度,
n=a或b,若
Figure FDA00021050508300000713
则n=a,否则,n=b,
u—通过反应
Figure FDA00021050508300000714
消耗的Mg2+浓度,
v—通过反应
Figure FDA00021050508300000715
生成的H2CO3浓度,
w—通过反应
Figure FDA00021050508300000716
生成的H2S浓度,
Ksp6—反应
Figure FDA00021050508300000717
的平衡常数。
10.根据权利要求9所述的计算方法,其特征在于,所述油田集输管道内采出液的pH值通过公式(58)计算得到:
Figure FDA00021050508300000718
式中:
n=a或b,若
Figure FDA00021050508300000719
则n=a;否则,n=b。
CN201910549378.XA 2019-06-24 2019-06-24 油田集输管道内采出液ph值的计算方法 Active CN112129707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910549378.XA CN112129707B (zh) 2019-06-24 2019-06-24 油田集输管道内采出液ph值的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910549378.XA CN112129707B (zh) 2019-06-24 2019-06-24 油田集输管道内采出液ph值的计算方法

Publications (2)

Publication Number Publication Date
CN112129707A true CN112129707A (zh) 2020-12-25
CN112129707B CN112129707B (zh) 2024-05-28

Family

ID=73849306

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910549378.XA Active CN112129707B (zh) 2019-06-24 2019-06-24 油田集输管道内采出液ph值的计算方法

Country Status (1)

Country Link
CN (1) CN112129707B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014155868A1 (ja) * 2013-03-29 2014-10-02 三菱電機株式会社 水質診断方法、水質診断装置および給湯器システム
CN106153415A (zh) * 2016-06-23 2016-11-23 郭子锐 一种二元弱酸酸式盐溶液离子浓度的处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014155868A1 (ja) * 2013-03-29 2014-10-02 三菱電機株式会社 水質診断方法、水質診断装置および給湯器システム
CN106153415A (zh) * 2016-06-23 2016-11-23 郭子锐 一种二元弱酸酸式盐溶液离子浓度的处理方法

Also Published As

Publication number Publication date
CN112129707B (zh) 2024-05-28

Similar Documents

Publication Publication Date Title
Clegg Activity coefficients in natural waters
Boyd et al. pH, carbon dioxide, and alkalinity
Sridhar et al. Effects of water and gas compositions on the internal corrosion of gas pipelines modeling and experimental studies
Christov Study of bromide salts solubility in the (m1KBr+ m2CaBr2)(aq) system at T= 323.15 K. Thermodynamic model of solution behaviour and (solid+ liquid) equilibria in the ternaries (m1KBr+ m2CaBr2)(aq), and (m1MgBr2+ m2CaBr2)(aq), and in the quinary (Na+ K+ Mg+ Ca+ Br+ H2O) systems to high concentration and temperature
Loewenthal et al. Control of corrosion and aggression in drinking water systems
Nešić et al. Implementation of a comprehensive mechanistic prediction model of mild steel corrosion in multiphase oil and gas pipelines
Christov Isopiestic investigation of the osmotic coefficients of aqueous CaBr2 and study of bromide salt solubility in the NaBr–CaBr2–H2O system at 50° C: Thermodynamic model of solution behavior and solid–liquid equilibria in the CaBr2–H2O, and NaBr–CaBr2–H2O systems to high concentration and temperature
Fosbøl et al. Review and recommended thermodynamic properties of FeCO3
Stiff Jr et al. A method for predicting the tendency of oil field waters to deposit calcium sulfate
CN104897521B (zh) 一种二氧化碳无水压裂液导流能力测试系统
Plennevaux et al. pH prediction in concentrated aqueous solutions under high pressure of acid gases and high temperature
Boyd et al. Carbon dioxide, pH, and alkalinity
Christov et al. Thermodynamic Modeling of Aqueous Aluminum Chemistry and Solid-Liquid Equilibria to High Solution Concentration and Temperature. I. The Acidic H-Al-Na-K-Cl-H 2 O System from 0 to 100 C
Okocha et al. Scale prediction for iron, zinc and lead sulphides and its relation to scale test design
CN112228019B (zh) 一种评价二氧化碳驱盐垢沉淀对采收率影响的方法
CN112129707A (zh) 油田集输管道内采出液ph值的计算方法
Tomson et al. Scale formation and prevention in the presence of hydrate inhibitors
Anderko et al. Non-Ideal Gases and Solutions, Complexes and Ion Pairs in Corrosion
CN115236118B (zh) 一种分析地热流体在管道内结垢情况的方法
JPH09318586A (ja) ステンレス鋼油井管の寿命予測方法
García-Antón et al. Use of differential pulse polarography to study corrosion of galvanized steel in aqueous lithium bromide solution
Ferguson Mineral scale prediction and control at extreme TDS
Bågander et al. In situ determination of the apparent solubility product of amorphous iron sulphide
Panthi Carbonate chemistry and calcium carbonate saturation state of rural water supply projects in Nepal
Norton et al. Ammonia in simulated Hanford double-shell tank wastes: Solubility and effects on surface tension

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant