1.一种基于T2谱分析的砂岩注水采油剩余油分布测定方法,其特征在于包括如下步骤:步骤一、对取芯岩心进行洗油、烘干、抽真空并饱和地层水,对饱和水状态下的岩心进行NMR测试,确定饱和水岩心T2谱Pw及T2谱的孔隙度累积曲线Sw;
步骤二、对饱和水岩心进行离心建立束缚水,对饱和水离心后的岩心进行NMR测试,确定饱和水离心后岩心T2谱Pwc及T2谱的孔隙度累积曲线Swc;
步骤三、烘干饱和水离心后的岩心、抽真空并饱和原油,对饱和油状态下的岩心进行NMR测试,确定饱和油岩心T2谱Po;
步骤四、对饱和油岩心进行洗油、烘干、抽真空并依次饱和地层水、原油,然后对初始饱和油状态下的岩心进行NMR测试,确定初始饱和油状态下的岩心T2谱Pio及T2谱的孔隙度累积曲线Sio;
步骤五、开展注水驱替/渗吸置换采油实验,采油结束后对岩心进行NMR测试,确定驱替/渗吸置换采油后岩心T2谱Pro及T2谱的孔隙度累积曲线Sro;
步骤六、对驱替/渗吸置换采油后岩心进行洗油、烘干,对岩心进行高压压汞测试,确定岩心的孔喉分布;
步骤七、基于实验数据计算得到岩心剩余油分布,具体计算方法为:
①依据Pw、Po及高压压汞测试结果完成水的横向弛豫时间T2w与岩心孔喉半径r的转换及油的横向弛豫时间T2o与r的转换,确定T2w‑r转换关系式及T2o‑r转换关系式;
②从Pio中读取信号量大于0的最大T2值to1,从Pw中读取信号量大于0的最大T2值twm,从Pw中读取大于twm且信号量为0的最小T2值to2,从Sio中读取to1、to2对应的累积孔隙度值sio1、sio2,从Sro中读取to1、to2对应的累积孔隙度值sro1、sro2,根据①确定的T2o‑r转换关系式将to1、to2转换为孔喉半径r1、r2,r1>r2,则半径为r2~r1孔喉中的初始含油量sio(1~2)=sio1‑sio2,剩余油含量sro(1~2)=sro1‑sro2,原油采收率k1~2=1‑sro(1~2)/sio(1~2);
③由①确定的T2w‑r转换关系式确定r1、r2分别对应的T2w值tw1、tw2,由①确定的T2o‑r转换关系式确定tw1、tw2分别对应的孔喉半径r3、r4,r3>r4,从Sio中读取tw1、tw2对应的累积孔隙度值sio3、sio4,从Swc中读取tw1、tw2对应的累积孔隙度值swc1、swc2,则半径为r4~r3孔喉中初始含油量sio(3~4)=sio3‑sio4‑(swc1‑swc2);从Sro中读取tw1、tw2对应的累积孔隙度值sro3、sro4,从Sw中读取tw1、tw2对应的累积孔隙度值sw1、sw2,则半径为r4~r3孔喉中剩余油含量sro(3~4)=sro3‑sro4‑[sw1‑sw2‑sro(1~2)],原油采收率k3~4=1‑sro(3~4)/sio(3~4);
④从Swc中读取twm对应的累积孔隙度值swcm,从Sw中读取twm对应的累积孔隙度值swm,则半径为r3~r2孔喉中初始含油量sio(2~3)=sio2‑sio3‑(swcm‑swc1),剩余油含量sro(2~3)=sro2‑sro3‑(swm‑sw1),原油采收率k2~3=1‑sro(2~3)/sio(2~3);
⑤由①确定的T2w‑r转换关系式确定r3、r4分别对应的T2w值tw3、tw4,由①确定的T2o‑r转换关系式确定tw3、tw4分别对应的孔喉半径r5、r6,r5>r6,从Sio中读取tw3、tw4对应的累积孔隙度值sio5、sio6,从Swc中读取tw3、tw4对应的累积孔隙度值swc3、swc4,则半径为r6~r5孔喉中初始含油量为sio(5~6)=sio5‑sio6‑(swc3‑swc4);从Sro中读取tw3、tw4对应的累积孔隙度值sro5、sro6,从Sw中读取tw3、tw4对应的累积孔隙度值sw3、sw4,则半径为r6~r5孔喉中剩余油含量sro(5~6)=sro5‑sro6‑[sw3‑sw4‑sro(3~4)],原油采收率k5~6=1‑sro(5~6)/sio(5~6);半径为r5~r4孔喉中初始含油量sio(4~5)=sio4‑sio5‑(swc2‑swc3),剩余油含量sro(4~5)=sro4‑sro5‑[sw2‑sw3‑sro(2~3)],原油采收率k4~5=1‑sro(4~5)/sio(4~5);
⑥采用步骤⑤的计算方法,依次计算半径为r2n+2~r2n+1、r2n+1~r2n孔喉中初始含油量、剩余油含量及原油采收率,n≥3,具体为:由T2w‑r转换关系式确定r2n‑1、r2n分别对应的T2w值tw2n‑1、tw2n,由T2o‑r转换关系式确定tw2n‑1、tw2n分别对应的孔喉半径r2n+1、r2n+2,r2n+1>r2n+2,从Sio中读取tw2n‑1、tw2n对应的累积孔隙度值sio2n+1、sio2n+2,从Swc中读取tw2n‑1、tw2n对应的累积孔隙度值swc2n‑1、swc2n,则半径为r2n+2~r2n+1孔喉中初始含油量为sio(2n+1~2n+2)=sio2n+1‑sio2n+2‑(swc2n‑1‑swc2n);从Sro中读取tw2n‑1、tw2n对应的累积孔隙度值sro2n+1、sro2n+2,从Sw中读取tw2n‑1、tw2n对应的累积孔隙度值sw2n‑1、sw2n,则半径为r2n+2~r2n+1孔喉中剩余油含量sro(2n+1~2n+2)=sro2n+1‑sro2n+2‑[sw2n‑1‑sw2n‑sro(2n‑1~2n)],原油采收率k2n+1~2n+2=1‑sro(2n+1~2n+2)/sio(2n+1~2n+2);
半径为r2n+1~r2n孔喉中初始含油量sio(2n~2n+1)=sio2n‑sio2n+1‑(swc2n‑2‑swc2n‑1),剩余油含量sro(2n~2n+1)=sro2n‑sro2n+1‑[sw2n‑2‑sw2n‑1‑sro(2n‑2~2n‑1)],原油采收率k2n~2n+1=1‑sro(2n~2n+1)/sio(2n~2n+1),确定未计算的其他尺寸孔喉内的剩余油含量原油采收率
式中szro为
根据岩心实验结果测定的剩余油量;szio为根据岩心实验结果测定的初始含油量;
sro(2n‑1~2n)为半径为r2n~r2n‑1孔喉中剩余油含量;sro2n为从Sro中读取的tw2n‑2对应的累积孔隙度值,Sro是驱替/渗吸置换采油后岩心T2谱的孔隙度累积曲线;swc2n‑2为从Swc中读取的tw2n‑2对应的累积孔隙度值,Swc是饱和水并离心后岩心T2谱的孔隙度累积曲线;sro(2n‑2~2n‑1)是半径为r2n‑1~r2n‑2孔喉中剩余油含量;sro(2n~2n+1)是半径为r2n+1~r2n孔喉中剩余油含量;
sro(2n+1~2n+2)是半径为r2n+2~r2n+1孔喉中剩余油含量;当kra<0.01时计算结束,未计算的其他尺寸孔喉内的原油未动用。