CN102809526B  Method for measuring diffusion coefficient of carbon dioxide in saturated oil core  Google Patents
Method for measuring diffusion coefficient of carbon dioxide in saturated oil core Download PDFInfo
 Publication number
 CN102809526B CN102809526B CN201210307962.2A CN201210307962A CN102809526B CN 102809526 B CN102809526 B CN 102809526B CN 201210307962 A CN201210307962 A CN 201210307962A CN 102809526 B CN102809526 B CN 102809526B
 Authority
 CN
 China
 Prior art keywords
 diffusion
 pressure
 curve
 carbon dioxide
 eff
 Prior art date
Links
 238000009792 diffusion process Methods 0.000 title claims abstract description 105
 CURLTUGMZLYLDIUHFFFAOYSAN carbon dioxide Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSAyMzUuMjMxLDEzNy40NTQgTCAyMTIuMjI3LDEzNy40NTQnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gMjEyLjIyNywxMzcuNDU0IEwgMTg5LjIyMywxMzcuNDU0JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0wJyBkPSdNIDIzNS4yMzEsMTYyLjU0NiBMIDIxMi4yMjcsMTYyLjU0Nicgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSAyMTIuMjI3LDE2Mi41NDYgTCAxODkuMjIzLDE2Mi41NDYnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiMzQjQxNDM7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTEnIGQ9J00gMTA5Ljc3MiwxMzcuNDU0IEwgODYuNzY3MywxMzcuNDU0JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDg2Ljc2NzMsMTM3LjQ1NCBMIDYzLjc2MywxMzcuNDU0JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDEwOS43NzIsMTYyLjU0NiBMIDg2Ljc2NzMsMTYyLjU0Nicgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSA4Ni43NjczLDE2Mi41NDYgTCA2My43NjMsMTYyLjU0Nicgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHRleHQgeD0nMjYzLjQ2JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8dGV4dCB4PScxMzgnIHk9JzE3MCcgY2xhc3M9J2F0b20tMScgc3R5bGU9J2ZvbnQtc2l6ZTo0MHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPkM8L3RleHQ+Cjx0ZXh0IHg9JzEyLjU0MDMnIHk9JzE3MCcgY2xhc3M9J2F0b20tMicgc3R5bGU9J2ZvbnQtc2l6ZTo0MHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjwvc3ZnPgo= data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSA3MS40OTQ3LDM4LjU1MDMgTCA1OS4xODU5LDM4LjU1MDMnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gNTkuMTg1OSwzOC41NTAzIEwgNDYuODc3LDM4LjU1MDMnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiMzQjQxNDM7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gNzEuNDk0Nyw0NS40NDk3IEwgNTkuMTg1OSw0NS40NDk3JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0wJyBkPSdNIDU5LjE4NTksNDUuNDQ5NyBMIDQ2Ljg3Nyw0NS40NDk3JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDM2Ljk5OCwzOC41NTAzIEwgMjQuNjg5MSwzOC41NTAzJyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDI0LjY4OTEsMzguNTUwMyBMIDEyLjM4MDIsMzguNTUwMycgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAzNi45OTgsNDUuNDQ5NyBMIDI0LjY4OTEsNDUuNDQ5Nycgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAyNC42ODkxLDQ1LjQ0OTcgTCAxMi4zODAyLDQ1LjQ0OTcnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+Cjx0ZXh0IHg9JzcyLjM1NzEnIHk9JzQ4Ljg5OTQnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTNweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8dGV4dCB4PSczNy44NjA0JyB5PSc0OC44OTk0JyBjbGFzcz0nYXRvbS0xJyBzdHlsZT0nZm9udC1zaXplOjEzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+QzwvdGV4dD4KPHRleHQgeD0nMy4zNjM2NCcgeT0nNDguODk5NCcgY2xhc3M9J2F0b20tMicgc3R5bGU9J2ZvbnQtc2l6ZToxM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjwvc3ZnPgo= O=C=O CURLTUGMZLYLDIUHFFFAOYSAN 0.000 title claims abstract description 65
 229910002092 carbon dioxide Inorganic materials 0.000 title claims abstract description 64
 239000001569 carbon dioxide Substances 0.000 title claims abstract description 64
 239000010779 crude oil Substances 0.000 claims abstract description 38
 239000003921 oil Substances 0.000 claims abstract description 26
 239000011435 rock Substances 0.000 claims description 47
 230000015572 biosynthetic process Effects 0.000 claims description 18
 238000005755 formation reaction Methods 0.000 claims description 18
 230000000875 corresponding Effects 0.000 claims description 14
 241001212149 Cathetus Species 0.000 claims description 12
 239000012530 fluid Substances 0.000 claims description 8
 238000007789 sealing Methods 0.000 claims description 6
 230000000694 effects Effects 0.000 claims description 5
 239000011148 porous material Substances 0.000 claims description 3
 239000000565 sealant Substances 0.000 claims description 3
 238000004088 simulation Methods 0.000 claims description 3
 238000000034 method Methods 0.000 abstract description 7
 238000004090 dissolution Methods 0.000 abstract 1
 229960004424 Carbon Dioxide Drugs 0.000 description 51
 239000000243 solution Substances 0.000 description 7
 XLYOFNOQVPJJNPUHFFFAOYSAN water Substances data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTAwLjUwMScgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+SDwvdGV4dD4KPHRleHQgeD0nMTI2LjExNCcgeT0nMTg2JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+MjwvdGV4dD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTMuMzEyOCcgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPkg8L3RleHQ+Cjx0ZXh0IHg9JzI4LjE1NjknIHk9JzYyLjg2MzYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTVweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4yPC90ZXh0Pgo8dGV4dCB4PSczNS4wNDU1JyB5PSc1My41OTA5JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== O XLYOFNOQVPJJNPUHFFFAOYSAN 0.000 description 4
 238000005516 engineering process Methods 0.000 description 3
 238000005259 measurement Methods 0.000 description 3
 238000011084 recovery Methods 0.000 description 3
 241001506308 Potato virus T Species 0.000 description 2
 238000002715 modification method Methods 0.000 description 2
 229920002102 polyvinyl toluene Polymers 0.000 description 2
 230000002123 temporal effect Effects 0.000 description 2
 238000010792 warming Methods 0.000 description 2
 208000008425 Protein Deficiency Diseases 0.000 description 1
 238000004364 calculation method Methods 0.000 description 1
 235000011089 carbon dioxide Nutrition 0.000 description 1
 238000006073 displacement reaction Methods 0.000 description 1
 238000002474 experimental method Methods 0.000 description 1
 238000002347 injection Methods 0.000 description 1
 239000007924 injection Substances 0.000 description 1
 230000005012 migration Effects 0.000 description 1
 230000004048 modification Effects 0.000 description 1
 238000006011 modification reaction Methods 0.000 description 1
 239000000126 substance Substances 0.000 description 1
 238000006467 substitution reaction Methods 0.000 description 1
 230000001629 suppression Effects 0.000 description 1
Abstract
The invention relates to a method for measuring a diffusion coefficient of carbon dioxide in a saturated oil core. The method comprises the following steps of: performing nondimensionalization by using a convectiondiffusion mathematical model, and performing numerical solution, wherein the diffusion coefficient Deff of the carbon dioxide in the saturated oil core is finally obtained by using the method and required to be used in the numerical solution process; and continuously correcting the value of D'eff, and repeatedly executing the step (13), so a curve c of delta Co and t<1/2> is completely coincided with a curve b of delta PTh and the t<1/2>, wherein the D'eff is the Deff and represents the diffusion coefficient of the carbon dioxide in the saturated oil core. The method has the advantages that the influence of the volumetric expansion of crude oil caused by dissolution of the carbon dioxide in the crude oil on the diffusion process and the influence of a porous medium on the diffusion process are taken into consideration, and the obtained effective diffusion coefficient of the carbon dioxide in the saturated oil core can accurately reflect the true diffusion process.
Description
Technical field
The present invention relates to a kind of method of measuring carbon dioxide coefficient of diffusion in saturated oil rock core, belong to the technical field of petrochemical complex.
Background technology
The global warming problem causing due to a large amount of discharges of carbon dioxide is increasingly severe.And collecting carbonic anhydride with bury (CO
_{2}capture and Storage) technology thought a kind of potential, alternative carbon dioxide dischargereduction scheme widely, can reduce the carbon dioxide content in atmosphere, alleviates climate warming problem.At carbon dioxide geological, bury in scheme, improve oil recovery factor (CO with carbon dioxide
_{2}EOR), not only can reach the effect of burying of carbon dioxide, and can improve the recovery ratio of Reservoir Crude Oil, therefore carbon dioxide raising oil recovery factor technology is widely used in oil field.In the process of employing carbondioxide flooding, the quality transfer law of carbon dioxide in saturated crude oil rock core seems particularly important for the prediction of the migration characteristics of injecting carbon dioxide.Therefore the mensuration of the coefficient of diffusion of carbon dioxide in saturated crude oil rock core has great importance for the development of carbon dioxide flooding oil tech.
Measure at present carbon dioxide method of coefficient of diffusion in porous medium and be mainly the PVT method that nonexpansion type spreads, the method has only been considered the diffusion of carbon dioxide in porous medium, do not have to consider, because the dissolving of carbon dioxide causes the impact of crude oil volumetric expansion on diffusion, can not well simulate the true diffusion process of carbon dioxide in saturated crude oil rock core.Meanwhile, the PVT of employing methods are measured the coefficient of diffusion of carbon dioxide in crude oil more at present, and the measurement of this coefficient of diffusion does not reflect the impact of porous medium on diffusion process, thereby can not be used for measuring the coefficient of diffusion of carbon dioxide at saturated crude oil rock core.
Summary of the invention
Summary of the invention
For above technical deficiency, the invention provides a kind of method of measuring carbon dioxide coefficient of diffusion in saturated oil rock core.The present invention, under the combined influence effect that considers crude oil volumetric expansion and porous medium, calculates the mathematical model of carbon dioxide coefficient of diffusion in saturated crude oil rock core by having derived in conjunction with one dimension advection diffusion equation and real gas state equation.By measure the Pressure Drop of carbon dioxide in conjunction with its corresponding calculated with mathematical model the coefficient of diffusion of carbon dioxide in saturated crude oil rock core.
The diffusion model of carbon dioxide in the rock core of saturation water and the difference of the diffusion model of carbon dioxide of the present invention in the rock core of saturated oil are: carbon dioxide solubility in water after, water expands hardly, be that the expansion of water is on not impact of diffusion, and the dissolving of carbon dioxide in oil can cause significantly expansion of oil in the present invention, in rock core, the expansion of oil has obvious influence to diffusion, so it is convection currentdiffusion mathematical model that the present invention will be spread mathematics model modification, but because the Analytical Solution of convection current diffusion mathematical model is more difficult, so the present invention is directed to convection currentdiffusion mathematical model, first carry out nondimensionalization, then carry out numerical solution, and in numerical solution process, need to use the diffusion coefficient D of the carbon dioxide of finally trying to achieve in saturated oil rock core
_{eff}, described D
_{eff}that the method for the invention solves and obtains.
Terminological interpretation:
Saturated fluid processing: refer to by vacuumizing the whole suckingoffs of air of loading in diffusion barrel in sample rock core interior detail fine porosity, then utilize constantflux pump that Experimental Flowing Object displacement is entered in diffusion barrel, injection experiments fluid generation build the pressure, it is logical when pressure is suppressed 10MPa, Experimental Flowing Object in diffusion barrel is under high pressure pressed in the hole of filling sample rock core, has realized whole saturated upper Experimental Flowing Objects in the hole that loads sample rock core.
Detailed Description Of The Invention
Technical scheme of the present invention is as follows:
A kind of carbon dioxide diffusion coefficient D in saturated oil rock core of measuring
_{eff}method, comprise that step is as follows:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as filling sample rock core, described filling sample rock core will be dried, after will upper and lower filling sample rock core both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, standby;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel
_{0}, gather n time point t
_{1}～t
_{n}pressure P
_{1}～P
_{n}, n is more than or equal to 2;
(10) Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}; By the pressure drop Δ P calculating
_{ex, 1}~ Δ P
_{ex, n}respectively with time point t
_{1}~ t
_{n}the corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records
_{ex}with t
^{1/2}curve a, calculate the slope k of described curve a cathetus section part;
(11) to Δ P
_{ex}with t
^{1/2}curve a revise: because carbon dioxide and the saturated crude oil at filling sample core surface place rigidly connect while touching, need to set up a stable liquidgas interface boundary condition, thereby the pressure drop of experiment incipient stage is very fast, the Pressure Drop of experiment starting stage belongs to the pressure drop Δ P that the unstable stage of diffusion causes measurement
_{ex}with t
^{1/2}curve map be upwards offset a segment distance, for the ease of calculating pressure falloff curve b with theory below, contrast, need to be to Δ P
_{ex}with t
^{1/2}curve a revise.Modification method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point
_{s}; Then by formula, 1. obtain average error
Formula 1. in, Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}～P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}, kPa; Kbe the slope of the curve a cathetus section part of trying to achieve in step (10),
t
_{i}for testing i the time point recording, min; t
_{s}be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P
_{ex, i}revise, obtain revising pressure drop values Δ P
_{co, i}:
Make Δ P
_{co}with t
^{1/2}curve c;
(12) bring the slope k in step (10) into formula 3.,
Formula 3. in, M
_{∞}the time levels off to carbon dioxide when infinite and diffuses into the amount in filling sample rock core, mol; Zcompressibility factor is determined by the pressure and temperature of testing; Runiversal gas constant, 8.314Pam
^{3}k
^{1}mol
^{1}; Texperimental temperature, K; D '
_{eff}carbon dioxide effective diffusion cofficient inunexpansive system, m
^{2}/ s; r
_{0}the xsect radius offilling sample rock core, m; Annular volume between Vfilling sample rock core and diffusion barrel, m
^{3}; Ttime point, comprises t
_{1}~ t
_{n}; KΔ P
_{ex}with t
^{1/2}the slope of the straightline segment part that forms;
By formula, 3. obtained,
obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system
_{eff};
(13) by D '
_{eff}bring in convection currentdiffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute
with crude oil volumetric expansion velocity distribution
(14) gas concentration lwevel is distributed
with crude oil volumetric expansion velocity distribution
there is bonding state equation after dimension, obtain corresponding theory and calculate pressure drop Δ P
_{th}, utilize described theory to calculate pressure drop Δ P
_{th}with t
^{1/2}curve plotting b, curve b and the Δ P that has just started experiment and record
_{co}with t
^{1/2}curve a have error, the source of error is exactly D '
_{eff};
(15) constantly revise D '
_{eff}value, repeating step (13), makes Δ P
_{co}with t
^{1/2}curve c and Δ P
_{th}with t
^{1/2}curve b overlap completely, D ' now
_{eff}be D
_{eff}, wherein said D
_{eff}for the coefficient of diffusion of carbon dioxide in saturated oil rock core.
The invention has the advantages that:
Present invention is directed at the diffusion process of carbon dioxide in saturated crude oil rock core, considered the crude oil volumetric expansion that causes due to the dissolving of carbon dioxide in the crude oil impact on diffusion process, considered again the impact of porous medium on diffusion process, the effective diffusion cofficient of the carbon dioxide of trying to achieve in saturated crude oil rock core can reflect real diffusion process comparatively accurately simultaneously.
Accompanying drawing explanation
Fig. 1 is that embodiment passes through the P in hisandhers watches 1
_{1}P
_{n}corresponding time t
_{1}~ t
_{n}mapping, i.e. pressure temporal evolution figure;
Fig. 2 is that embodiment passes through the pressure differential deltap P in hisandhers watches 1
_{1}Δ P
_{n}to time t
^{1/2}mapping, tests pressure differential deltap P
_{ex}with t
^{1/2}variation a, and through revised Pressure Drop Δ P
_{co}with t
^{1/2}curve c, by Convectiondiffusion model iterative computation Δ P out
_{th}with t
^{1/2}curve b.
Embodiment
According to embodiment and Figure of description, the present invention is described in detail below, but is not limited to this.
Embodiment,
As shown in Figure 12.
A kind of carbon dioxide diffusion coefficient D in saturated oil rock core of measuring
_{eff}method, comprise that step is as follows:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as filling sample rock core, described filling sample rock core will be dried, after will upper and lower filling sample rock core both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, standby;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel
_{0}, gather n time point t
_{1}~ t
_{n}pressure P
_{1}～P
_{n}, (n is more than or equal to 2);
(10) Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}; By the pressure drop Δ P calculating
_{ex, 1}~ Δ P
_{ex, n}respectively with time point t
_{1}~ t
_{n}the corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records
_{ex}with t
^{1/2}curve a, calculate the slope of described curve a cathetus section part
(11) to Δ P
_{ex}with t
^{1/2}curve a revise: because carbon dioxide and the saturated crude oil at filling sample core surface place rigidly connect while touching, need to set up a stable liquidgas interface boundary condition, thereby the pressure drop of experiment incipient stage is very fast, the Pressure Drop of experiment starting stage belongs to the pressure drop Δ P that the unstable stage of diffusion causes measurement
_{ex}with t
^{1/2}curve map be upwards offset a segment distance, for the ease of calculating pressure falloff curve b with theory below, contrast, need to be to Δ P
_{ex}with t
^{1/2}curve a revise.Modification method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point
_{s}=19.27min; Then by formula, 1. obtain average error
,
$\stackrel{\‾}{\mathrm{\δ}}=245.31\left(\mathrm{kPa}\right)$
Formula 1. in, Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}～P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}, kPa; Kis the slope of the curve a cathetus section part of trying to achieve in step (10),
t
_{i}for testing i the time point recording, min; t
_{s}be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P
_{ex, i}revise, obtain revising pressure drop values Δ P
_{co, i}:
Make Δ P
_{co}with t
^{1/2}curve c;
Table 1: utilize pressure transducer to gather the pressure P of each time point
_{1}～P
_{n}with the pressure differential deltap P calculating
_{ex}with correction pressure differential deltap P
_{co}
By the pressure in hisandhers watches 1, the time is mapped, obtain pressure temporal evolution figure, as shown in Figure 1.
(12) bring the slope k in step (10) into formula 3.,
Formula 3. in, M
_{∞}the time levels off to carbon dioxide when infinite and diffuses into the amount in filling sample rock core, mol; Zcompressibility factor is determined by the pressure and temperature of testing; Runiversal gas constant, 8.314Pam
^{3}k
^{1}mol
^{1}; Texperimental temperature, K; D '
_{eff}carbon dioxide effective diffusion cofficient inunexpansive system, m
^{2}/ s; r
_{0}the xsect radius offilling sample rock core, m; Annular volume between Vfilling sample rock core and diffusion barrel, m
^{3}; Ttime point, comprises t
_{1}~ t
_{n}; KΔ P
_{ex}with t
^{1/2}the slope of the straightline segment part that forms;
By formula, 3. obtained,
$k=\frac{4{M}_{\∞}\mathrm{ZRT}\sqrt{{D}_{\mathrm{eff}}^{\′}}}{{r}_{0}V\sqrt{\mathrm{\π}}}=24.3033\mathrm{kPa}/\sqrt{\mathrm{min}},$ Obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system
_{eff}=2.95 × 10
^{10}m
^{2}/ s;
(13) by D '
_{eff}bring in convection currentdiffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute
with crude oil volumetric expansion velocity distribution
Formula 4. in, dimensionless time:
dimensionless radius:
nondimensional velocity:
dimensionless concentration:
$\stackrel{\‾}{c}=\frac{c}{{c}_{0}};$ $\mathrm{\λ}=\stackrel{\‾}{u}\frac{1}{\stackrel{\‾}{r}};$
Consider that dimensionless diffusionconvection equation fully implicit solution finite difference method that crude oil expands carries out solution procedure as follows, in departure process, the single order of speed single order differential and concentration and secondorder differential central difference schemes, the single order differential of time adopts forward difference form.Abovementioned equation is carried out to discretize, discretize form:
Applying finite difference scheme partial differential equation is above write as:
Wherein:
${a}_{i}=\frac{\mathrm{\Δ\τ}}{\mathrm{\Δ}{\stackrel{\‾}{r}}^{2}}\frac{\mathrm{\λ\Δ\τ}}{2\mathrm{\Δ}\stackrel{\‾}{r}},$ ${b}_{i}=1+\frac{\mathrm{\Δ\τ}}{2\mathrm{\Δ}\stackrel{\‾}{r}}({\stackrel{\‾}{u}}_{i+1}^{n+1}{\stackrel{\‾}{u}}_{i1}^{n+1})+\frac{\mathrm{\Δ\τ}}{{\stackrel{\‾}{r}}_{i}}{\stackrel{\‾}{u}}_{i}^{n+1}+\frac{2\mathrm{\Δ\τ}}{\mathrm{\Δ}{\stackrel{\‾}{r}}^{2}},$ ${e}_{i}=\frac{\mathrm{\λ\Δ\τ}}{2\mathrm{\Δ}\stackrel{\‾}{r}}\frac{\mathrm{\Δ\τ}}{\mathrm{\Δ}{\stackrel{\‾}{r}}^{2}},$
Boundary condition discretize can obtain:
The system of equations of substitution boundary condition discretize can be write as:
Your alternative manner of application GaussSaden calculates concentration and the velocity distribution of each time step.To each time step, by the concentration of previous time step as the velocity distribution in this time step of calculation of initial value.Adopt velocity distribution in this time step to calculate the CONCENTRATION DISTRIBUTION of next time step.Then by new CONCENTRATION DISTRIBUTION, calculate the velocity distribution in next time step.The condition that abovementioned iterative process finishes is that the concentration maximum relative error of each radial direction Nodes is less than limits of error e (e=10 in the present invention
^{4}).
(14) gas concentration lwevel is distributed
with crude oil volumetric expansion velocity distribution
there is bonding state equation after dimension, obtain corresponding theory and calculate pressure drop Δ P
_{th}, utilize described theory to calculate pressure drop Δ P
_{th}with t
^{1/2}curve plotting b, curve b and the Δ P that has just started experiment and record
_{co}with t
^{1/2}curve a have error, the source of error is exactly D '
_{eff}:
The dimensionless gas concentration lwevel of trying to achieve in step (13) is distributed
and crude oil volumetric expansion velocity distribution
Δ P falls in the pressure carbon dioxide that has dimension and obtain intumescent system
_{th}t in time
^{1/2}variation relation
Formula 5. in, Zcompressibility factor by the pressure and temperature of testing determine; Runiversal gas constant, 8.314Pam
^{3}k
^{1}mol
^{1}; Texperimental temperature, K; P
_{0}carbon dioxide original pressure indiffusion barrel, MPa; Annular volume between Vrock sample and diffusion barrel, m
^{3}; Δ Vobtained by 6. formula, crude oil volumetric expansion does not cause the reduction of diffusion barrel inner volume, m in the same time
^{3}; Qdo not diffuse into the in the same time CO in saturated oil rock core
_{2}be dissolved in the CO that expands in crude oil
_{2}amount of substance sum, mol.
Formula 6. in, r
_{0}the xsect radius offilling sample rock core, m; The length of hfilling sample rock core, m;
nondimensional velocity; τdimensionless time.
The Pressure Drop Δ P of the Convectiondiffusion model of trying to achieve
_{th}with time relationship as following table 2,
Table 2: Δ P falls in Convectiondiffusion model calculating pressure
_{th}with time relationship table
(15) constantly revise D '
_{eff}value, repeating step (13), makes Δ P
_{co}with t
^{1/2}curve a and Δ P
_{th}with t
^{1/2}curve b overlap completely, D ' now
_{eff}be D
_{eff}, wherein said D
_{eff}for the coefficient of diffusion of carbon dioxide in saturated oil rock core:
By the Δ P doing in step (10)
_{co}with t
^{1/2}curve c and step (14) in do Δ P
_{th}with t
^{1/2}curve b contrast, as Fig. 2.If two curve differs larger, revise effective diffusion cofficient D '
_{eff}and get back in step (13).If two curve relative error is less than limits of error e (e=10 in this method
^{4}), stop iterative process, effective diffusion cofficient D now
_{eff}==3.19 × 10
^{10}m
^{2}/ s is the effective diffusion cofficient of carbon dioxide in saturated crude oil rock core.
Claims (1)
1. measure carbon dioxide diffusion coefficient D in saturated oil rock core for one kind
_{eff}method, it is characterized in that, it is as follows that the method comprising the steps of:
(1) diffusion barrel is dried;
(2) will sample cylinder rock core or artificial cylinder rock core as loading core sample, described filling core sample will be dried, after will upper and lower fillings core sample both ends of the surface sealing with fluid sealant, pack in diffusion barrel;
(3) diffusion barrel is vacuumized;
(4) utilize well heater to heat diffusion barrel, constant temperature is to wanting simulated formation temperature, standby;
(5) in diffusion barrel, pump into Experimental Flowing Object, and be forced into 10MPa, to loading core sample, carry out saturated fluid processing; Described Experimental Flowing Object is simulated formation crude oil, and described simulated formation crude oil is the Simulation of Crude Oil identical with formation pore oil property;
(6) continue diffusion barrel inside Experimental Flowing Object to carry out pressurized operation, be forced into and want simulated formation pressure;
(7) in the bottom of described diffusion barrel, be provided with check valve, adjust the set pressure of check valve, described set pressure is compared with the little 0.1MPa of diffusion barrel internal pressure;
(8) to diffusion barrel, pump into carbon dioxide, the pressure of described carbon dioxide with want simulated formation pressure identical, due to the effect of check valve, the Experimental Flowing Object in emptying diffusion barrel;
(9) sealing diffusion barrel, utilizes pressure transducer to gather the original pressure P in diffusion barrel
_{0}, gather n time point t
_{1}～t
_{n}pressure P
_{1}～P
_{n}, n is more than or equal to 2;
(10) Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}; By the pressure drop Δ P calculating
_{ex, 1}～Δ P
_{ex, n}respectively with time point t
_{1}～t
_{n}the corresponding mapping of 1/2 power, obtain the Pressure Drop Δ P that experiment records
_{ex}with t
^{1/2}curve a, calculate the slope k of described curve a cathetus section part;
(11) to Δ P
_{ex}with t
^{1/2}curve a revise, method is as follows: first find the flex point of described curve a, i.e. the corresponding time point t of curve a cathetus section part starting point
_{s}; Then by formula, 1. obtain average error
:
Formula 1. in, Δ P
_{ex, i}for testing the pressure P of n the time point recording
_{1}～P
_{n}respectively with original pressure P
_{0}pressure differential, as, Δ P
_{ex, 1}=P
_{1}P
_{0}, Δ P
_{ex, 2}=P
_{2}P
_{0}..., Δ P
_{ex, n}=P
_{n}P
_{0}, kPa; Kbe the slope of the curve a cathetus section part of trying to achieve in step (10),
t
_{i}for testing i the time point recording, min; t
_{s}be the corresponding time point of curve a cathetus section part starting point of trying to achieve in step (10), min;
According to 2. formula to experimental pressure drop Δ P
_{ex, i}revise, obtain revising pressure drop values Δ P
_{co, i}:
Make Δ P
_{co}with t
^{1/2}curve c;
(12) bring the slope k in step (10) into formula 3.,
Formula 3. in, M
_{∞}the time levels off to carbon dioxide when infinite and diffuses into the amount of loading in core sample, mol; Zcompressibility factor, is determined by the pressure and temperature of testing; Runiversal gas constant, 8.314Pam
^{3}k
^{1}mol
^{1}; Texperimental temperature, K; D '
_{eff}carbon dioxide effective diffusion cofficient inunexpansive system, m
^{2}/ s; r
_{0}the xsect radius offilling core sample, m; Annular volume between Vfilling core sample and diffusion barrel, m
^{3}; Ttime point, comprises t
_{1}～t
_{n}; KΔ P
_{ex}with t
^{1/2}the slope of the straightline segment part that forms;
By formula, 3. obtained,
obtain the effective diffusion cofficient D ' of carbon dioxide in unexpansive system
_{eff};
(13) by D '
_{eff}bring in convection currentdiffusion mathematical model of nondimensionalization and carry out numerical solution, obtain nondimensional gas concentration lwevel and distribute
with crude oil volumetric expansion velocity distribution
;
(14) gas concentration lwevel is distributed
with crude oil volumetric expansion velocity distribution
have after dimension in conjunction with real gas state equation, obtain corresponding theory and calculate pressure drop Δ P
_{th}, utilize described theory to calculate pressure drop Δ P
_{th}with t
^{1/2}curve plotting b, curve b and the Δ P that has just started experiment and record
_{co}with t
^{1/2}curve a have error, the source of error is exactly D '
_{eff};
(15) constantly revise D '
_{eff}value, repeating step (13), makes Δ P
_{co}with t
^{1/2}curve c and Δ P
_{th}with t
^{1/2}curve b overlap completely, D ' now
_{eff}be D
_{eff}, wherein said D
_{eff}for the coefficient of diffusion of carbon dioxide in saturated oil rock core.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201210307962.2A CN102809526B (en)  20120827  20120827  Method for measuring diffusion coefficient of carbon dioxide in saturated oil core 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201210307962.2A CN102809526B (en)  20120827  20120827  Method for measuring diffusion coefficient of carbon dioxide in saturated oil core 
Publications (2)
Publication Number  Publication Date 

CN102809526A CN102809526A (en)  20121205 
CN102809526B true CN102809526B (en)  20140416 
Family
ID=47233287
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201210307962.2A CN102809526B (en)  20120827  20120827  Method for measuring diffusion coefficient of carbon dioxide in saturated oil core 
Country Status (1)
Country  Link 

CN (1)  CN102809526B (en) 
Families Citing this family (14)
Publication number  Priority date  Publication date  Assignee  Title 

CN104975827B (en) *  20140403  20180119  中国石油化工股份有限公司  Predict the material balance method of carbon dioxide flooding oil reservoir index 
CN104237079B (en) *  20141010  20170118  中国石油大学(华东)  Experimental method and device for measuring diffusion coefficient of carbon dioxide in water 
CN104502236B (en) *  20141215  20170322  中国石油大学(华东)  Method for measuring diffusion coefficient and equilibrium concentration of CO2 in process of diffusion from water phase to oil phase 
WO2017181392A1 (en) *  20160421  20171026  深圳市樊溪电子有限公司  Method for measuring diffusion coefficient of carbon dioxide in crude oil 
CN107314950B (en) *  20160427  20191105  中国石油化工股份有限公司  A method of measurement carbon dioxide diffusion coefficient in porous media 
CN106596353B (en) *  20161219  20190226  中国科学院武汉岩土力学研究所  Gas diffusivity measuring device and method in coal cinder based on differential pressure method 
CN106840973B (en) *  20170216  20190604  中国石油大学(华东)  CO in a kind of test porous media2The device and its test method of diffusion concentration and diffusion coefficient 
CN106872594B (en) *  20170216  20181116  中国石油大学(华东)  CO in a kind of test saturated oils porous media2The method of concentration distribution and diffusion coefficient 
CN107063935B (en) *  20170421  20180717  中国石油大学(华东)  A kind of measurement CO2Diffusion process and the method for diffusion coefficient between wateroil phase 
CN108801860B (en) *  20170502  20210101  中国石油化工股份有限公司  Method for measuring diffusion coefficient of carbon dioxide in crude oil 
CN108362614B (en) *  20180119  20200327  中国石油大学(华东)  Device and method for measuring diffusion coefficient in huff and puff process of shale oil carbon dioxide 
CN109269943B (en) *  20181127  20210420  衡阳师范学院  Method for rapidly measuring effective diffusion coefficient of radon in film 
CN109283098B (en) *  20181127  20210420  衡阳师范学院  Method for measuring effective diffusion coefficient of radon in film 
CN109883891B (en) *  20190325  20200327  中国石油大学(北京)  Measuring device and method for injected gasliquid phase diffusion coefficient of fracturecavity oil reservoir 
Family Cites Families (5)
Publication number  Priority date  Publication date  Assignee  Title 

JP4915903B2 (en) *  20060131  20120411  国立大学法人横浜国立大学  Method and apparatus for measuring oxygen diffusion coefficient of porous material 
CN101929939B (en) *  20100806  20130313  中联煤层气国家工程研究中心有限责任公司  Method for detecting diffusion coefficient of coal bed gas in coal matrix 
CN101975718B (en) *  20100813  20121114  中国科学院山西煤炭化学研究所  Method for simultaneously measuring highpressure gas adsorption capacity and adsorption swell capacity of coal petrography and measuring equipment 
CN102121891A (en) *  20101216  20110713  西南石油大学  Device and method for testing molecular diffusion coefficients in multicomponent oilgas system under high temperature and high pressure 
CN202195998U (en) *  20110916  20120418  东北石油大学  Device for testing diffusion coefficient of natural gas 

2012
 20120827 CN CN201210307962.2A patent/CN102809526B/en active IP Right Grant
Also Published As
Publication number  Publication date 

CN102809526A (en)  20121205 
Similar Documents
Publication  Publication Date  Title 

CN104101564B (en)  A kind of method of unstable state High Temperature High Pressure test flow in low permeability core freeboundary problem  
Gasda et al.  Vertical equilibrium with subscale analytical methods for geological CO 2 sequestration  
Clarkson et al.  Production analysis of tightgas and shalegas reservoirs using the dynamicslippage concept  
CN104502224B (en)  Saturation water Coal Under rock isothermal desorption curve determination device and method  
Shook et al.  Determining reservoir properties and flood performance from tracer test analysis  
CN103323352B (en)  Natural gas hydrate deposit dynamic triaxial mechanicacousticelectrical synchronous test experimental device and method  
CN103246820B (en)  A kind of numerical simulation for oilgas reservoir computing method  
CN106021778B (en)  A kind of simulation CO2The determination method of displacement dynamic miscible pressure  
Li et al.  Diffusion coefficients of supercritical CO2 in oilsaturated cores under low permeability reservoir conditions  
CN103257151B (en)  In a kind of quantitative evaluation oil and gas secondary migration process, pore throat employs the method for rule  
CN104018829B (en)  A kind of method that utilization coal bed gas well creation data measures air water phase percolation curve  
Liu et al.  Coal permeability evolution and gas migration under nonequilibrium state  
Li et al.  Production behavior of methane hydrate in porous media using huff and puff method in a novel threedimensional simulator  
Connell et al.  History matching of enhanced coal bed methane laboratory core flood tests  
CN103278615B (en)  Test method for geological storage of carbon dioxide coal seams  
CN202300351U (en)  Long rock core model displacement device  
CN202330233U (en)  Experiment test device for permeability of rock core under condition of formation pressure  
CN103279647A (en)  Method for calculating tight sandstone reservoir permeability based on hole feature parameters  
CN104568694B (en)  Method for testing gaswater relative permeability of dense core  
CN203630018U (en)  Unstablestate determination device for permeability of hyposmosis rock  
CN105840187A (en)  Method for calculating staged fracturing productivity of compact reservoir horizontal well  
CN104895550B (en)  A kind of tight gas pressure break horizontal well numerical well testing model establishes method for solving  
CN102914494B (en)  Device for measuring dynamic leakoff of foam fracturing fluid and working method thereof  
CN103527185A (en)  Horizontalwell physical simulation experiment device and experimental method thereof  
CN102353625B (en)  Method for measuring overburden porosity with water in permeation fluid mechanics experiment 
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 