CN104213899B  A kind of Logging Identification Method of formation rock skeleton  Google Patents
A kind of Logging Identification Method of formation rock skeleton Download PDFInfo
 Publication number
 CN104213899B CN104213899B CN201310218843.4A CN201310218843A CN104213899B CN 104213899 B CN104213899 B CN 104213899B CN 201310218843 A CN201310218843 A CN 201310218843A CN 104213899 B CN104213899 B CN 104213899B
 Authority
 CN
 China
 Prior art keywords
 value
 porosity
 rock
 rock matrix
 density
 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.)
 Active
Links
 239000011435 rock Substances 0.000 title claims abstract description 124
 210000002356 Skeleton Anatomy 0.000 title claims abstract description 56
 230000015572 biosynthetic process Effects 0.000 title claims abstract description 41
 238000005755 formation reaction Methods 0.000 title claims abstract description 41
 239000011159 matrix material Substances 0.000 claims abstract description 78
 238000004458 analytical method Methods 0.000 claims abstract description 76
 238000004364 calculation method Methods 0.000 claims description 10
 239000012530 fluid Substances 0.000 claims description 7
 230000000875 corresponding Effects 0.000 claims description 6
 238000005259 measurement Methods 0.000 claims description 4
 238000010835 comparative analysis Methods 0.000 claims description 2
 239000011148 porous material Substances 0.000 claims description 2
 230000037250 Clearance Effects 0.000 claims 1
 230000035512 clearance Effects 0.000 claims 1
 238000011161 development Methods 0.000 abstract description 3
 239000004575 stone Substances 0.000 description 10
 238000010586 diagram Methods 0.000 description 3
 238000011160 research Methods 0.000 description 3
 238000000034 method Methods 0.000 description 2
 210000000988 Bone and Bones Anatomy 0.000 description 1
 BVKZGUZCCUSVTDUHFFFAOYSAL Carbonate dianion Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSAyMjYuMzg4LDY0LjY4MzUgTCAxODUuODM1LDg4LjA5NjcnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTAnIGQ9J00gMTg1LjgzNSw4OC4wOTY3IEwgMTQ1LjI4MywxMTEuNTEnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiMzQjQxNDM7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTEnIGQ9J00gMTQ1LjI4MywxMTEuNTEgTCAxMDQuNzMsODguMDk2Nycgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAxMDQuNzMsODguMDk2NyBMIDY0LjE3Nyw2NC42ODM1JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0yJyBkPSdNIDEzMS4zNDEsMTExLjUxIEwgMTMxLjM0MSwxNTkuNTM0JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoyLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0yJyBkPSdNIDEzMS4zNDEsMTU5LjUzNCBMIDEzMS4zNDEsMjA3LjU1OCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSAxNTkuMjI0LDExMS41MSBMIDE1OS4yMjQsMTU5LjUzNCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6Mi4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSAxNTkuMjI0LDE1OS41MzQgTCAxNTkuMjI0LDIwNy41NTgnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjIuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+Cjx0ZXh0IHg9JzI1NC4wMjEnIHk9JzYxLjgwMTMnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8dGV4dCB4PScyODEuNjIxJyB5PSc0NS44MDEzJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+LTwvdGV4dD4KPHRleHQgeD0nMTIuNTQzOCcgeT0nNjEuODAxMycgY2xhc3M9J2F0b20tMicgc3R5bGU9J2ZvbnQtc2l6ZTo0MHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjx0ZXh0IHg9JzQwLjE0MzgnIHk9JzQ1LjgwMTMnIGNsYXNzPSdhdG9tLTInIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4tPC90ZXh0Pgo8dGV4dCB4PScxMzMuMjgzJyB5PScyNzAuOTI3JyBjbGFzcz0nYXRvbS0zJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHBhdGggY2xhc3M9J2JvbmQtMCcgZD0nTSA2Ni45MDgsMTUuNzA4OCBMIDUzLjQ4OCwyMy40NTY4JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0wJyBkPSdNIDUzLjQ4OCwyMy40NTY4IEwgNDAuMDY4MSwzMS4yMDQ5JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0xJyBkPSdNIDQwLjA2ODEsMzEuMjA0OSBMIDI2LjY0ODEsMjMuNDU2OCcgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6IzNCNDE0MztzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMScgZD0nTSAyNi42NDgxLDIzLjQ1NjggTCAxMy4yMjgyLDE1LjcwODgnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiNFODQyMzU7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTInIGQ9J00gMzYuMzQ1NiwzMS4yMDQ5IEwgMzYuMzQ1Niw0Ny4xMTg0JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojM0I0MTQzO3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8cGF0aCBjbGFzcz0nYm9uZC0yJyBkPSdNIDM2LjM0NTYsNDcuMTE4NCBMIDM2LjM0NTYsNjMuMDMxOScgc3R5bGU9J2ZpbGw6bm9uZTtmaWxsLXJ1bGU6ZXZlbm9kZDtzdHJva2U6I0U4NDIzNTtzdHJva2Utd2lkdGg6MS4wcHg7c3Ryb2tlLWxpbmVjYXA6YnV0dDtzdHJva2UtbGluZWpvaW46bWl0ZXI7c3Ryb2tlLW9wYWNpdHk6MScgLz4KPHBhdGggY2xhc3M9J2JvbmQtMicgZD0nTSA0My43OTA1LDMxLjIwNDkgTCA0My43OTA1LDQ3LjExODQnIHN0eWxlPSdmaWxsOm5vbmU7ZmlsbC1ydWxlOmV2ZW5vZGQ7c3Ryb2tlOiMzQjQxNDM7c3Ryb2tlLXdpZHRoOjEuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjEnIC8+CjxwYXRoIGNsYXNzPSdib25kLTInIGQ9J00gNDMuNzkwNSw0Ny4xMTg0IEwgNDMuNzkwNSw2My4wMzE5JyBzdHlsZT0nZmlsbDpub25lO2ZpbGwtcnVsZTpldmVub2RkO3N0cm9rZTojRTg0MjM1O3N0cm9rZS13aWR0aDoxLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxJyAvPgo8dGV4dCB4PSc2Ny44Mzg2JyB5PScyMC4wMzc1JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjE0cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPHRleHQgeD0nNzguMTEyNicgeT0nMTQuMDgxNScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZTo5cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+LTwvdGV4dD4KPHRleHQgeD0nMy4zNjM2NCcgeT0nMjAuMDM3NScgY2xhc3M9J2F0b20tMicgc3R5bGU9J2ZvbnQtc2l6ZToxNHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPk88L3RleHQ+Cjx0ZXh0IHg9JzEzLjYzNzYnIHk9JzE0LjA4MTUnIGNsYXNzPSdhdG9tLTInIHN0eWxlPSdmb250LXNpemU6OXB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPi08L3RleHQ+Cjx0ZXh0IHg9JzM1LjYwMTEnIHk9Jzc1Ljg3NDQnIGNsYXNzPSdhdG9tLTMnIHN0eWxlPSdmb250LXNpemU6MTRweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K [O]C([O])=O BVKZGUZCCUSVTDUHFFFAOYSAL 0.000 description 1
 235000013399 edible fruits Nutrition 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
 229910052500 inorganic mineral Inorganic materials 0.000 description 1
 239000011707 mineral Substances 0.000 description 1
 239000004576 sand Substances 0.000 description 1
Abstract
The invention provides a kind of Logging Identification Method of formation rock skeleton, belong to oil exploration, development field.The inventive method includes：The depth point that resistivity is found when tending to infinity, regard the log value of the depth point as rock matrix value；The rock matrix value includes interval transit time, compensation density and compensated neutron measured value；For the exploratory area for thering is core analysis to test, the check analysis of rock matrix value is carried out, the rock matrix value demonstrated by core experiment is obtained；The rock matrix value demonstrated by core experiment is applied to the unknown exploratory area of same lithology.The present invention is solved for target zone rock matrix is indefinite or mineralogical composition is complicated stratum, the problem of the problem of reservoir porosity can not accurately be solved is difficult to ask for rock matrix, the cost for obtaining complex lithology rock matrix is greatly reduced, the computational accuracy of reservoir porosity is improved.
Description
Technical field
The invention belongs to oil exploration, development field, it is related to the logging technique that identification reservoir porosity is calculated, and in particular to
A kind of Logging Identification Method of formation rock skeleton.
Background technology
Acoustic velocity logging and density log are the most frequently used lithologyporosity logging methods.Want to be determined with porosity logging
The lithology and porosity on stratum, just must be set up the response equation of acoustic velocity logging and density log.
At present popular acoustic velocity logging response equation, density log response equation form it is as follows successively：
It can thus be concluded that the porosity of clean rock is：
In formula：For rock porosity, decimal；Δ t is the acoustic wave of rock time difference of measurement, and unit is μ s/m；Δt_{ma}For rock
Stone skeleton interval transit time, unit is μ s/m；Δt_{f}For the interval transit time of blowhole fluid, unit is μ s/m；For density hole
Porosity, decimal；ρ_{b}To measure the bulk density of rock, unit is g/cm^{3}；ρ_{ma}For matrix density, unit is g/cm^{3}；ρ_{f}For
The bulk density of blowhole fluid, unit is g/cm^{3}。
Can be seen that from the calculation formula of porosity and accurately ask for acoustic wave of rock skeleton value or density skeleton value is accurately to ask
Take the necessary condition of formation porosity.
For the less Within Monominerals rock of shale content, the sound wave skeleton or density skeleton of rock are easier to determine, but
Because oilgas exploration and development object is complicated, the lithology of target zone tends to complicated, and rock matrix is really permanent often to turn into well log interpretation
New problem.Such as volcanic rock is different because of volcanic rock mineralogical composition, and the volcanic rocks skeletons of different regions exists inconsistent existing
As such as shown in table 1 (Karamay oilfield difference volcanic rock porosity relational expression and lithoelectric parameters (according to Xu Chunhua, 2007)).
A kind of simple and accurate log analysis method of the problems demand is solved.
Table 1
Solution is mainly geochemical well logging now and analysis of coring obtains rock matrix, and not only cost is high for the above method,
And be difficult to use, it is necessary to the replacement of cost effective method comprehensively in big work area.
The content of the invention
It is an object of the invention to solve problem present in abovementioned prior art, there is provided a kind of survey of formation rock skeleton
Well recognition methods, easy to operate, conclusion is reliable.
The present invention is achieved by the following technical solutions：
A kind of Logging Identification Method of formation rock skeleton, including：
The depth point that resistivity is found when tending to infinity, regard the log value of the depth point as rock matrix value；It is described
Rock matrix value includes interval transit time, compensation density and compensated neutron measured value；
For the exploratory area for thering is core analysis to test, the check analysis of rock matrix value is carried out, obtains testing by core experiment
The rock matrix value demonstrate,proved；
The rock matrix value demonstrated by core experiment is applied to the unknown exploratory area of same lithology.
The check analysis of the rock matrix value includes：
(1) the rock matrix value, is intersected into obtained rock bone with sound wave or density with using core analysis porosity
Frame value does comparative analysis, obtains both uniformity；
(2) the rock matrix value, is inputted into porosity calculation formula, formation porosity is asked for；In known core hole,
The computational accuracy of formation porosity according to the porosity of core analysis to being tried to achieve is examined again.
The utilization core analysis porosity intersects what obtained rock matrix value was realized in sound wave or density：
Core analysis porosity is classified by lithology：Under same lithology, according to the depth of core analysis porosity
Angle value, finds acoustic logging value or density log value corresponding under the depth in log data；
Core analysis porosity value and the acoustic logging value are done twodimentional cross plot, then with after fitting a straight line, y is obtained
=kx+b relation, wherein y refer to interval transit time, and x refers to porosity, and when porosity x is zero, b is exactly interval transit time skeleton
Value；
Equally, core analysis porosity value and the density log value are done twodimentional cross plot, then with after fitting a straight line,
Y=kx+b relation is obtained, wherein y refers to density, and x refers to porosity, and when porosity x is zero, b is exactly density skeleton value；
If the apparent density of core analysis, then the density value obtained in borehole log data is directly substituted with apparent density.
The porosity calculation formula is as follows：
Wherein,
In formula：For rock porosity, decimal；Δ t is the acoustic wave of rock time difference of measurement, and unit is μ s/m；Δt_{ma}For rock
Stone skeleton interval transit time, unit is μ s/m；Δt_{f}For the interval transit time of blowhole fluid, unit is μ s/m；For density hole
Porosity, decimal；ρ_{b}To measure the bulk density of rock, unit is g/cm^{3}；ρ_{ma}For matrix density, unit is g/cm^{3}；ρ_{f}For
The bulk density of blowhole fluid, unit is g/cm^{3}。
The computational accuracy refers to calculate obtained formation porosity and core analysis porosity one using rock matrix value
The height of cause property.
It is this that the computational accuracy of formation porosity of the porosity according to core analysis to being tried to achieve is examined again
What sample was realized：
Obtained formation porosity will be calculated and do relative error analysis with core analysis porosity, if relative error magnitudes are small
In minimum reference value, then the rock core skeleton value is the rock matrix value demonstrated by core experiment.
The minimum reference value takes 10%.
It is described that the rock matrix value demonstrated by core experiment is applied into the unknown exploratory area of same lithology is such reality
Existing：
The rock matrix value demonstrated by core experiment is substituted into the porosity calculation formula, same lithology is obtained
The porosity in unknown exploratory area；The unknown exploratory area of same lithology refers to the lithology identical with the exploratory area for having core analysis to test
Well, the interval transit time skeleton value or density skeleton value or neutron skeleton value in the exploratory area are unknown.
The resistivity tends to infinity and refers to that formation resistivity is more than 2000 ohm meters.
Compared with prior art, the beneficial effects of the invention are as follows：
The present invention is solved for target zone rock matrix is indefinite or mineralogical composition is complicated stratum, it is impossible to accurate to solve
The problem of the problem of reservoir porosity, is difficult to ask for rock matrix；
This method is workable, the porosity calculated using this method, the core analysis porosity through known core hole
Verified, research conclusion accurately and reliably, greatly reduces the cost for obtaining complex lithology rock matrix, improves reservoir pore space
The computational accuracy of degree.
Brief description of the drawings
Fig. 11 is the rock matrix figure of tufa stone in log identification stratum.
Fig. 12 is the rock matrix figure of rhyolite in log identification stratum.
Fig. 21 is the core analysis density and porosity graph of a relation of rhyolite.
Fig. 22 is the core analysis density and porosity graph of a relation of tufa stone.
Fig. 3 is that X borehole logging tools explain porosity and the porosity comparison diagram of core analysis.
Fig. 4 is the step block diagram of the inventive method.
Fig. 5 is the porosity and core analysis porosity relative error schematic diagram that X wells are calculated.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings：
Formation rock skeleton ask for for it is accurate ask for formation porosity for, one to be emphasis in well log interpretation and
Difficult point.Because target zone rock matrix is indefinite or mineralogical composition is complicated, and use geochemical well logging or analysis acquisition rock of coring
Skeleton, not only cost is high, and is difficult to use comprehensively in big work area, so research speculates the rock for studying stratum with log
Stone skeleton, is a new approach.The foundation of this method is that tight formation porosity goes to zero, and resistivity tends to be infinitely great,
The porosity for now thinking the measuring point is negligible, therefore the measuring point can be considered formation rock skeletal point.This side
Method is workable, and in tight sand, carbonate rock and igneous reservoirs can all have preferable application value there is provided can
The formation rock skeleton acquiring method leaned on.
What the present invention was realized in：For target zone rock matrix is indefinite or mineralogical composition bad ground, it can attempt
The rock matrix for studying stratum is speculated using log.It speculates that foundation is that tight formation porosity goes to zero, and resistance
Rate tends to be infinitely great, and the porosity for now thinking the measuring point is negligible, therefore the measuring point can be considered formation rock
Skeletal point.The formation rock skeleton value obtained using this method, can be intersected with using core analysis porosity with sound wave or density
Obtained skeleton value is contrasted；The porosity precision calculated using this method, can use the core analysis hole of known core hole
Porosity is verified, so as to realize the purpose accurately calculated for the porosity on unknown rock matrix stratum.
The well logging recognition of rock matrix：Lowest porosity log value corresponding to pure lithology point can be considered close to rock matrix
Log value (the rock matrix log value is the theoretical value speculated, and the value is confirmed by lot of experimental data.Such as Fig. 21 and figure
Shown in 22).When its corresponding resistivity tends to infinity, it is believed that porosity is ignored, the log value of this depth point
(i.e. interval transit time, compensation density and compensated neutron measured value) can be considered rock matrix value.
As shown in figure 4, the inventive method includes：
(1) well logging of rock matrix value is determined：The depth that formation resistivity is more than 2000 ohm meters is read in practical application
Point corresponding interval transit time, compensation density and compensated neutron measured value (these three values are all obtained from log), this
Three values are rock matrix value；According to log quality or other aspects in general.If this three curve qualities are all
It is very good, three skeleton values can be read；If there is log quality problems, then uptostandard log is chosen,
Read skeleton value.
(2) check analysis of rock matrix value.The checking is divided into two steps：
First, the rock matrix value for determining the above method is obtained with being intersected using core analysis porosity with sound wave or density
Rock matrix value (core analysis porosity intersects with sound wave and is obtained by rock matrix：First have to a core analysis hole
Porosity presses lithology different classifications.Under same lithology, according to the depth value of core analysis porosity, found in log data
Corresponding acoustic logging value under the depth.Then core analysis porosity value and sound wave log value are done twodimentional cross plot, one
As can be appreciated that lithologic analysis porosity value increases with the increase of interval transit time, both approximately linear relations.Intended with straight line
After conjunction, y=kx+b relation is obtained.General y refers to interval transit time, and x refers to porosity.So when porosity x is zero, b is just
It is interval transit time skeleton value.Equally, core analysis porosity is intersected with density is also obtained by density skeleton value.Certainly such as
Fruit has the apparent density of core analysis, and the density value obtained in borehole log data can be directly substituted with apparent density.) do to score
Analysis (the relatively direct rock matrix value rock matrix value that namely this method introduction is obtained that is obtained by log and logical
Cross the uniformity that core analysis porosity intersects obtained rock matrix value with sound wave, density.Present invention side is confirmed by contrasting
The feasibility of method.)；
2nd, for the exploratory area for having core analysis to test, the rock matrix value demonstrated by core experiment is inputted into hole
Calculation formula is spent, formation porosity is asked for.In known core hole, according to the porosity of core analysis to abovementioned porosity calculation
Precision (rock matrix calculates obtained formation porosity and the height of core analysis porosity uniformity) is examined (to use again
This kind of method calculates obtained formation porosity and does relative error analysis with core analysis porosity, and relative error average is less than
10% is believed that the formation porosity degree of accuracy that the rock matrix obtained with this kind of method calculating is obtained is higher.Such as accompanying drawing 5
Shown, 9.11%) its relative error average is.
Porosity calculation formula is as follows：
Wherein,
In formula：For rock porosity, decimal；Δ t is the acoustic wave of rock time difference of measurement, and unit is μ s/m；Δt_{ma}For rock
Stone skeleton interval transit time, unit is μ s/m；Δt_{f}For the interval transit time of blowhole fluid, unit is μ s/m；For density hole
Porosity, decimal；ρ_{b}To measure the bulk density of rock, unit is g/cm^{3}；ρ_{ma}For matrix density, unit is g/cm^{3}；ρ_{f}For
The bulk density of blowhole fluid, unit is g/cm^{3}。
(3) with the application in the unknown exploratory area of lithology：When abovementioned well logging calculates porosity and core analysis porosity goodness of fit height
(formation porosity does relative error analysis with core analysis porosity, and relative error average is less than 10% i.e. it is believed that the goodness of fit
It is high), meeting well log interpretation precision, (formation porosity does relative error with core analysis porosity, and its average thinks less than 10%
It is to meet well log interpretation precision) when, this method can be used for the unknown exploratory area of lithology (lithology identical well.The skeleton value of reading
It is same lithology with uniformity.It is unknown to refer to that interval transit time skeleton value or density skeleton value or neutron skeleton value are unknown.)
Porosity logging calculate, otherwise can not use this method.Other methods can only be explored again.
The present invention obtains the rock matrix of complex lithology by using log, solves and is difficult to ask for complex lithology
The present situation of rock matrix is taken, this method is workable, the porosity calculated using this method, the rock core through known core hole point
Analysis porosity is verified that research conclusion accurately and reliably, greatly reduces the cost for obtaining complex lithology rock matrix, improved
The computational accuracy of reservoir porosity.
1st, log determines stratum rock matrix
Fig. 11 and Fig. 12 are the schematic diagrames that log recognizes formation rock skeleton.When the resistivity of fine and close lithology is more than
During 2000 ohm meter, it is believed that its porosity can be neglected, porosity measured value now is lithology skeleton value.According to the party
The different lithology skeleton value of 2 classes can be read in method, practical study on log.Density is respectively 2.62g/cm^{3}With
2.67g/cm^{3}.It is combined afterwards with geology, confirmation density skeleton is rhyolite for 2.62g/cm3's, and density skeleton is 2.67g/
Cm3's is that (geology determines that the area has rhyolite and two kinds of tufa stone is main to tufa stone by core observation combination rockmineral analysis
Lithology.The depth point of the abovementioned two skeleton of well logging recognition is respectively in both lithology coring depth sections.).
2nd, the formation rock skeleton and core analysis porosity that well logging is determined intersect obtained rock matrix with density and done pair
Than analysis
As shown in Figure 2, for rhyolite, intersected, linearly closed with core analysis porosity with core analysis density
It is to be(linear relationship is the porosity and core analysis with core analysis to coefficient correlation
These scatterplots fit what is come after density intersection.Can all there is individual coefficient correlation after fitting, to show the goodness of fit of fitting.) R=
0.879.WhenWhen, rhyolite density skeleton value is 2.626g/cm^{3}, the density skeleton 2.62g/cm determined with well logging^{3}Kissing
Close.For tufa stone, intersected with core analysis porosity with core analysis density, obtaining linear relationship isCoefficient R^{2}=0.9893.WhenWhen, tufa stone density skeleton value is 2.67g/
cm^{3}, the density skeleton 2.67g/cm determined with well logging^{3}It coincide very much.Therefore deduce that, for rock matrix is indefinite or ore deposit
The target zone of thing complicated component, this method that rock matrix value is determined with log is accurate believable.
3rd, the porosity precision calculated according to the porosity of core analysis the rock matrix determined with well logging is examined again
Test.
As seen from Figure 3, porosity and core analysis porosity (right number that stratum rock matrix is calculated are determined with well logging
4th column, shown in stain) it is basically identical.Therefore deduce that, the formation rock skeleton determined with logging well accurately and reliably, drops significantly
The low cost for obtaining complex lithology rock matrix, improves the computational accuracy of reservoir porosity.
Abovementioned technical proposal is one embodiment of the present invention, for those skilled in the art, at this
On the basis of disclosure of the invention application process and principle, it is easy to make various types of improvement or deformation, this is not limited solely to
Invent the method described by abovementioned embodiment, therefore previously described mode is preferred, and and without limitation
The meaning of property.
Claims (4)
1. a kind of Logging Identification Method of formation rock skeleton, it is characterised in that：Methods described includes：
The depth point that resistivity is found when tending to infinity, regard the log value of the depth point as rock matrix value；The rock
Skeleton value includes interval transit time, compensation density and compensated neutron measured value；
For the exploratory area for thering is core analysis to test, the check analysis of rock matrix value is carried out, obtains demonstrating by core experiment
Rock matrix value；
The rock matrix value demonstrated by core experiment is applied to the unknown exploratory area of same lithology,
Wherein, the check analysis of the rock matrix value includes：
(1) the rock matrix value, is intersected into obtained rock matrix value with sound wave or density with using core analysis porosity
Comparative analysis is done, both uniformity are obtained；
(2) the rock matrix value, is inputted into porosity calculation formula, formation porosity is asked for；In known core hole, according to
The computational accuracy of formation porosity of the porosity of core analysis to being tried to achieve is examined again,
Wherein, it is this that the computational accuracy of formation porosity of the porosity according to core analysis to being tried to achieve is examined again
What sample was realized：Obtained formation porosity will be calculated and do relative error analysis with core analysis porosity, if relative error magnitudes
Less than minimum reference value, then the rock matrix value is the rock matrix value demonstrated by core experiment；
Wherein, the minimum reference value takes 10%；
Wherein, the rock matrix value demonstrated by core experiment is applied into the unknown exploratory area of same lithology to be realized in
's：The rock matrix value demonstrated by core experiment is substituted into the porosity calculation formula, same lithology is obtained unknown
The porosity in exploratory area；The unknown exploratory area of same lithology refers to the lithology identical well with the exploratory area for having core analysis to test,
The interval transit time skeleton value or density skeleton value or neutron skeleton value in the exploratory area are unknown；
Wherein, the resistivity tend to infinity refer to formation resistivity be more than 2000 ohm meters；
Wherein, uptostandard log is chosen, the log value for reading the depth point is used as the rock matrix value.
2. the Logging Identification Method of formation rock skeleton according to claim 1, it is characterised in that：The utilization rock core point
Analysis porosity intersects what obtained rock matrix value was realized in sound wave or density：
Core analysis porosity is classified by lithology：Under same lithology, according to the depth value of core analysis porosity,
Acoustic logging value or density log value corresponding under the depth is found in log data；
Core analysis porosity value and the acoustic logging value are done twodimentional cross plot, then with after fitting a straight line, y=kx is obtained
+ b relation, wherein y refer to interval transit time, and x refers to porosity, and when porosity x is zero, b is exactly interval transit time skeleton value；
Equally, core analysis porosity value and density log value are done twodimentional cross plot, then with after fitting a straight line, obtains y=
Kx+b relation, wherein y refer to density, and x refers to porosity, and when porosity x is zero, b is exactly density skeleton value；
If the apparent density of core analysis, then the density value obtained in borehole log data is directly substituted with apparent density.
3. the Logging Identification Method of formation rock skeleton according to claim 2, it is characterised in that：The porosity calculation
Formula is as follows：
Wherein,
In formula：For rock porosity, decimal；Δ t is the acoustic wave of rock time difference of measurement, and unit is μ s/m；Δt_{ma}For rock matrix
Interval transit time, unit is μ s/m；Δt_{f}For the interval transit time of blowhole fluid, unit is μ s/m；It is small for density porosity
Number；ρ_{b}To measure the bulk density of rock, unit is g/cm^{3}；ρ_{ma}For matrix density, unit is g/cm^{3}；ρ_{f}For rock pore
The bulk density of clearance flow body, unit is g/cm^{3}。
4. the Logging Identification Method of formation rock skeleton according to claim 1, it is characterised in that：The computational accuracy is
Refer to the height that obtained formation porosity and core analysis porosity uniformity is calculated using rock matrix value.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201310218843.4A CN104213899B (en)  20130604  20130604  A kind of Logging Identification Method of formation rock skeleton 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201310218843.4A CN104213899B (en)  20130604  20130604  A kind of Logging Identification Method of formation rock skeleton 
Publications (2)
Publication Number  Publication Date 

CN104213899A CN104213899A (en)  20141217 
CN104213899B true CN104213899B (en)  20170804 
Family
ID=52095768
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201310218843.4A Active CN104213899B (en)  20130604  20130604  A kind of Logging Identification Method of formation rock skeleton 
Country Status (1)
Country  Link 

CN (1)  CN104213899B (en) 
Families Citing this family (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN104948164B (en) *  20150505  20170616  中国海洋石油总公司  The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter 
CN105158802B (en) *  20150821  20181123  西南石油大学  Lacustrine Basins Gravity Flow Sediments well logging quantitative identification method 
CN107301255B (en) *  20160415  20210105  中国石油化工股份有限公司  Method for calculating reservoir fracture porosity 
CN107795320B (en) *  20160905  20210105  中国石油化工股份有限公司  Calculation method for horizontal well carbonate reservoir parameters 
CN107605471B (en) *  20170721  20200609  中国石油天然气股份有限公司  Lakefacies carbonate rock lithology well logging quantitative identification method 
CN107762483B (en) *  20171102  20201013  中国石油天然气集团公司  Fluid identification method based on correlation coefficient and envelope area of logging curve 
CN110146409B (en) *  20190509  20200707  中国地质大学（武汉）  Method for determining shale gas saturation based on density 
Family Cites Families (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN102096107B (en) *  20091209  20121017  中国石油天然气股份有限公司  Method for evaluating permeability of reservoir layer according to interval transit time and density inversed pore flat degree 

2013
 20130604 CN CN201310218843.4A patent/CN104213899B/en active Active
Also Published As
Publication number  Publication date 

CN104213899A (en)  20141217 
Similar Documents
Publication  Publication Date  Title 

CN104636588B (en)  Calculate the method and device of organic carbon content in hydrocarbon source rock  
Tingay et al.  Evidence for overpressure generation by kerogentogas maturation in the northern Malay Basin  
CN102678106B (en)  Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument  
CN103792338B (en)  A kind of determination method and apparatus of content of organic carbon of hydrocarbon source rock  
CN103046868B (en)  Integrated horizontal well geosteering method  
Passey et al.  A practical model for organic richness from porosity and resistivity logs  
Chaudhry  Oil well testing handbook  
CN103256046B (en)  Unconventionaloil pool hides method and the device that horizontal well stitches the simulation of long fracturing parameter entirely  
US8217337B2 (en)  Evaluating a reservoir formation  
CN103149589B (en)  Pyrogenic rock Petroleum Exploration Methods and device  
US7363164B2 (en)  Method of evaluating fluid saturation characteristics in a geological formation  
CN102426390B (en)  Method for determining reserve volume of nonhomogeneous sandstone reservoir  
US6721661B2 (en)  Method of distinguishing types of geologic sedimentation  
EP2965126B1 (en)  Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data  
CN103592690B (en)  Method based on electric imaging logging porosity spectrum information automatic Reservoir Recognition crack  
Close et al.  Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada  
Zang et al.  World stress map database as a resource for rock mechanics and rock engineering  
Cassiani et al.  Modeling unsaturated flow in a layered formation under quasisteady state conditions using geophysical data constraints  
US8359184B2 (en)  Method, program and computer system for scaling hydrocarbon reservoir model data  
US9341557B2 (en)  Method and system for permeability calculation using production logs for horizontal wells, using a downhole tool  
EP2798376B1 (en)  Insitu characterization of formation constituents  
CN103760081A (en)  Gas pool prediction method and system of carbonate reservoir based on pore structure characteristics  
Worthington  Net Pay: What is it? What does it do? How do we quantify it? How do we use it?  
Lai et al.  Brittleness index estimation in a tight shaly sandstone reservoir using well logs  
EP1672392B1 (en)  Method for determining the water saturation of an underground formation 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
GR01  Patent grant 