一种测井数据的溯源方法与流程

1.本发明涉及一种测井数据的溯源方法，属于地球物理勘探数据处理技术领域。


背景技术：

2.测井数据在石油地质研究、油气田勘探和开发领域有着重要的作用，其记录了随深度不断变化岩层的电、磁、声、核等参数的数值。利用测井数据可以解释评价储层性质，计算油气储层的孔隙度、渗透率和饱和度等数据。测井数据是油田生产重要的涉密数据，如何对这些数据资产进行溯源，提高测井数据的信息安全，是当前亟需解决的课题。为了管理和利用测井数据，部分油田建立了具有测井数据检索和下载的信息系统。虽然系统可以记录下载用户的用户名、ip地址等信息，但是由于不同用户下载的测井数据相同，导致一旦测井数据发生泄密，无法准确溯源泄密测井数据的用户。
3.目前，数据的溯源，主要使用标注法对数据进行溯源。通过标注法对数据进行溯源，其技术实现在于将指示不同用户信息的标注与数据一起传播，当数据泄密需要溯源时，可以通过标注溯源泄密数据的用户。例如标注a指示用户甲，当甲下载数据时，其下载的数据内将参杂标注a的信息。一旦用户甲下载的数据发生泄密，通过数据内参杂的标注a，很容易溯源泄密用户为甲。而测井数据为可编辑的文本格式，如果使用标注法对测井数据进行溯源，授权用户可通过文本编辑很容易擦除标注信息，导致无法准确溯源泄密测井数据的用户。


技术实现要素：

4.本发明的目的是提供一种测井数据的溯源方法，以解决目前的标注溯源法无法有效对测井数据进行溯源的问题。
5.本发明为解决上述技术问题而提供一种测井数据的溯源方法，该方法包括以下步骤：
6.1)获取原始测井数据并进行存储；
7.2)当有授权用户下载测井数据时，根据该授权用户的下载信息生成若干位整数信息，所述的授权用户的下载信息包括用户名、登录ip、所下载测井数据信息和下载时间；
8.3)将存储的原始测井数据中的至少一条测井曲线进行分段处理，在每段选取部分深度进行叠加处理，将叠加处理后的测井数据提供给该授权用户；
9.叠加处理是在选取的部分深度的各深度点对应测井曲线上增加相应的叠加量，所述叠加量是根据生成的整数信息按照构建的叠加函数确定的，；所述的叠加函数是基于原始测井曲线对应深度的数值、整数信息和叠加系数确定的，叠加系数与原始测井曲线对应深度的数值有关；
10.4)当需要对测井数据进行溯源时，将待溯源的测井数据上传给后台服务器，后台服务器将待溯源的测井数据中的各测井曲线与原始的测井数据中的各测井曲线进行对照比较，找到经过叠加处理的测井曲线，从中确定出相应的叠加量和，并根据所述叠加函数确
定出对应的表征用户下载信息的整数信息，基于该整数信息确定对应的用户下载信息，从而实现了对测井数据的溯源。
11.本发明通过利用授权用户下载时的下载信息生成整数信息，根据整数信息确定和原始测井数据确定出相应的叠加量，将测井曲线进行分段处理，在每段选取部分深度按照得到的叠加量进行叠加处理；在需要对测井数据进行溯源时，将待溯源的测井数据与原始测井数据进行比较，确定出叠加量并反推出相应的整数信息，由此该整数信息确定出相应用户下载信息，从而实现对该测井数据的溯源。通过上述过程，本发明能够准确地实现对测井数据的溯源。此外，本发明叠加量的计算是基于测井曲线原始数据，因而叠加后的数据不会对测井数据的应用造成影响；叠加处理时不会更改测井数据体的结构，因而不会影响到专业地质软件对叠加后的测井数据体进行解析及应用。
12.进一步地，所述步骤3)中若进行叠加处理的测井曲线为至少两条时，不同测井曲线选取的叠加处理对应深度段不同。
13.进一步地，所述步骤3)中若进行叠加处理的测井曲线为至少两条时，不同测井曲线选取的叠加处理的深度段相同。
14.进一步地，所述的叠加函数为：
15.y
result
＝y
source
+k*n
[0016][0017][0018]
其中y
source
为原始测井曲线某个深度点的数值，y
result
为叠加后的测井曲线对应深度点的数值，n为整数信息分解后的某位数值，k为叠加系数。
[0019]
进一步地，所述的叠加函数为：
[0020]yresult
＝y
source-k*n
[0021][0022][0023]
其中y
source
为原始测井曲线某个深度点的数值，y
result
为叠加后的测井曲线对应深度点的数值，n为整数分解后的某位数值，k为叠加系数。
[0024]
进一步地，为了有效保证溯源的若干位整数不会被滤波程序和曲线分割程序破坏，所述步骤3)叠加处理是按照设定叠加重复量和叠加间隔量进行的，所述的叠加重复量指的是按深度顺序连续对测井数据进行相同的叠加处理的个数，所述的叠加间隔量指的是对测井数据进行叠加处理的间隔数。
[0025]
进一步地，所述的叠加重复量为3，叠加间隔量为2。
附图说明
[0026]
图1是本发明测井数据溯源方法的流程图。
具体实施方式
[0027]
下面结合附图对本发明的具体实施方式作进一步地说明。
[0028]
本发明通过表征下载用户信息的整数信息确定原始测井数据的叠加量，其中确定的叠加小于待叠加的测井曲线段的原始测井数据的5％，按照确定叠加对原始测井数据进行叠加，将叠加后的测井数据提供给该下载用户进行下载；当需要对测井数据溯源时，将待溯源的测井数据上传给服务器，由服务器将待溯源的测井数据与原始测井数据进行比较，计算出叠加的偏移值，根据该偏移值对应的整数信息确定出下载的用户信息，从而实现对测井数据的溯源，该方法的实现流程如图1所示，下面结合具体的测井数据对本发明实施过程进行详细说明。
[0029]
1.对原始测井数据进行存储。
[0030]
本实施例将原始测井数据a存储到后台服务器中，由后台服务器对测井数据a进行管理。本实施例选用某井4214.500米至4217.375米的测井数据作为原始数据，如表1所示，该图中的第一列为深度数据，第二列为ac曲线，第三列为azim曲线，第四列为devi数据。
[0031]
表1
[0032]
[0033][0034]
2.生成表征用户下载的整数信息。
[0035]
当某个授权用户下载该井的测井数据时，后台服务器根据该授权用户的信息生成若干位整数a，假设本实施例生成5位整数39384，该数据记录了登录用户及其下载测井数据的信息，包括但不仅限于用户名、登录ip、所下载测井数据信息和下载时间。
[0036]
3.构建叠加函数。
[0037]
本发明利用整数信息x对用户要下载的测井原始数据中的测井曲线进行叠加处理。将x位整数按位数进行分解，确定与原始测井数据a叠加的偏移值n1、n2、n3
……
nx，且其偏移值小于待叠加的测井曲线段的原始测井数据a的5％。为此，本发明构建了相应的叠加函数。
[0038]
本实施例中利用的叠加函数有两个参数，分别为叠加重复量和叠加间隔量，本实施例中默认叠加重复量a为3，叠加间隔量b为2，其中叠加重复量a为3指的是按深度顺序连续对3个深度点的测井数据进行相同的叠加处理，叠加间隔量b为2指的是每隔两个深度点对测井数据进行叠加处理。
[0039]
叠加处理采用的叠加函数为：
[0040]yresult
＝y
source
+k*n
[0041]
本实例中的整数为39384，则分解完后n有5个取值，分别为3、9、3、8和4，k为叠加系数，k的取值与原始测井数据的大小有关，具体采用的确定公式为：
[0042][0043][0044]
其中y
source
为原始测井曲线某个深度点的数值，数值大小按照测井曲线的标准单位来取值，例如假设测井曲线为ac曲线，其单位为。
[0045]
作为其他实施方式，叠加处理采用的叠加函数也可以为：
[0046]yresult
＝y
source-k*n
[0047]
4.将存储的原始测井数据中的至少一条测井曲线进行分段处理。
[0048]
本实施例中的测井曲线包括有ac曲线、azim曲线和devi曲线，本发明可以对这三条曲线都进行叠加处理，也可以只对其中的某一条或者两条进行叠加处理。在进行叠加处理时，需要将测井曲线进行分段，在每一段选取其中部分深度点进行叠加处理。
[0049]
若需要进行叠加处理的测井曲线不止一条，则在进行分段时，不同测井曲线的分段可以不同，也可以相同。
[0050]
本实施例中将ac曲线分成了三段，4178.5-4181.75为第一段、4193.375-4196.625为第二段、4220.75-4224.25为第三段；其他曲线可以按照ac曲线的分段方式，也可以采用其他的分段方式，例如每一段对应的深度范围与ac曲线不同；段数可以是四段、五段等。
[0051]
5根据生成的整数信息按照构建的叠加函数对选取的测井曲线各段中的部分深度点进行叠加处理。
[0052]
对本实施例而言，假设生成的整数信息为39384，一共有五位，将整数逐位分解，本实例是将5位整数39384分解为3(n1)、9(n2)、3(n3)、8(n4)和4(n5)，按照设定的叠加函数、叠加重复量和叠加间隔量对其中任意一条测井数据的各段进行叠加处理。
[0053]
以表1中的ac曲线为例，本实施例中的ac曲线取值是在40-60之间，可以确定出对应的r值在0.1-1之间，按照叠加系数k的公式，此时的叠加系数k取值为0.1，即在对ac曲线的各段进行叠加处理时，叠加系数k＝0.1。同理，假设测井曲线为azim曲线，则叠加处理时的叠加系数k＝1；假设测井曲线为devi曲线，则叠加处理时的叠加系数k＝0.01。
[0054]
以ac曲线第一段为例，由于默认叠加重复量a为3，叠加间隔量b为2，按照深度顺序第一个深度叠加点选取4178.625，第二个深度叠加点选取4178.75，第三个深度叠加点选取4178.875，然后间隔两个深度点(4179和4179.125)后再选取三个深度叠加点，以此类推，每隔两个深度点选取三个深度叠加点。
[0055]
在具体叠加时，由于本实施例中的整数被分解成3、9、3、8和4，即n可以等于3、9、3、8和4，根据叠加函数，在原始测井曲线的上的叠加量可以是0.3、0.9、0.3、0.8和0.4，本发明在前三个深度叠加点的ac测井曲线的上均叠加0.3，按照上述顺序，第四、五和六深度叠加点的ac测井曲线的上均叠加0.9，第七、八和九深度叠加点的ac测井曲线的上均叠加0.3，第十、十一和十二深度叠加点的ac测井曲线的上均叠加0.8，第十三、十四和十五深度叠加点
的ac测井曲线的上均叠加0.4，最后三个深度点的ac测井曲线均不进行叠加。通过上述过程可完成对ac测井曲线的第一段的叠加处理，结果如表2所示。
[0056]
表2
[0057]
#depthac原始数据ac叠加处理后ac差值4178.556.07156.0710.0004178.62542.35842.6580.3004178.7553.73154.0310.3004178.87557.08557.3850.300417942.84842.8480.0004179.12542.07742.0770.0004179.2543.08143.9810.9004179.37554.77855.6780.9004179.553.70554.6050.9004179.62543.23743.2370.0004179.7542.69342.6930.0004179.87540.41340.7130.300418055.22955.5290.3004180.12546.62146.9210.3004180.2547.99147.9910.0004180.37541.36541.3650.0004180.548.82349.6230.8004180.62545.40546.2050.8004180.7542.67943.4790.8004180.87543.41043.4100.000418151.84451.8440.0004181.12552.75553.1550.4004181.2558.74059.1400.4004181.37547.39647.7960.4004181.546.28946.2890.0004181.62557.36957.3690.0004181.7549.63849.6380.000
[0058]
通过同样的方式，可以对其他两端进行叠加处理，结果如表3和表4所示。
[0059]
表3
[0060]
[0061][0062]
表4
[0063][0064][0065]
按照同样的方式，可以实现对其他测井曲线中每段的叠加处理，需要注意的是，需根据测井曲线数值的大小确定叠加的数量级，即确定k值。
[0066]
本实施例中上述过程均在后台服务器中实现，后台服务器将服务器中叠加后的所有测井曲线a重新打包为含有上述原始测井数据和叠加偏移值的测井数据c，并提供给当前
操作的用户下载。
[0067]
5.进行数据溯源。
[0068]
当需要对测井数据进行溯源时，将待溯源的测井数据c上传给后台服务器，后台服务器将待溯源的测井数据c与原始的测井数据a进行对照比较，确定出经过叠加的数据段，从中检测出相应的叠加量，将叠加量进行组合，并按照设定的叠加重复量和叠加间隔量以及叠加函数确定出对应的表征用户下载信息的整数x，通过该整数在服务器中检索到下载该测井数据的用户名、登录ip、所下载测井数据信息和下载时间，实现对测井数据的溯源。
[0069]
对本实施例而言，若待溯源的测井数据为表2中经过叠加处理的ac测井数据，将表2中的测井数据和表1中的同样深度的原始测井数据进行比较，可以确定出相应的叠加量，如表2中的ac差值所在的列。根据叠加重复量3和叠加间隔量2可以确定出其包含有多个整数位，其中前三个叠加量为0.3，对应的叠加系数为0.1，根据叠加函数确定出前三个叠加量对应的整数位为3，同理，可以确定出后续的整数位分别为9、3、8和4，将各整数位组合就可以得到整数39384，基于该整数信息确定对应的用户下载信息，从而实现了对测井数据的溯源。
[0070]
通过上述过程，本发明通过对测井曲线进行叠加，并将叠加量与n位整数进行关联，从而实现测井数据的溯源的方法。本发明在对测井曲线进行叠加时，不会更改测井数据体的结构，因而不会影响到专业地质软件对叠加后的测井数据体进行解析及应用；其叠加量的计算是基于测井曲线原始数据，因而叠加后的数据不会对测井数据的应用造成影响；同时设置了叠加重复量和叠加间隔量，能有效保证溯源的n位整数不会被滤波程序和曲线分割程序破坏，提高了溯源的成功率。