CN119415811A 融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法 （中国石油大学(华东)）

融合光纤应变监测和裂缝扩展模拟的压裂本发明公开了一种融合光纤应变监测和裂2建立基于状态变量的状态控制方程，并根据建立的状态控制方程根据裂缝扩展模型的系统误差的协方差矩阵，计算所述状态根据监测到的光纤应变数据、计算出的下一时间步状态变量对应的对裂缝单元状态进行更新，并在更新完后再次进入到预测2.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，;Dnn为n时刻的缝宽估计值；3.如权利要求2所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，4.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，3;;5.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，6.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，;n+1为n+1时刻的观测方程的系数矩阵；7.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，;n+1为n+1时刻监测到的光纤应变数据；Dnrn为步骤S300中预测出的n+1时刻的缝宽估计值；为更新后的n+1时刻的缝宽估计值，即为通过状态更新方程计算出来的n+14;n+1为n+1时刻的观测方程的系数矩阵；Ie为应变的测量误差的协方差矩阵；根据计算出的状态控制方程的下一时间步预测误差的协方差;8.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，9.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，10.如权利要求1所述的融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法，5[0004]上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技[0005]本发明的主要目的在于提供了一种融合光纤应变监测和裂缝扩展模拟的压裂裂对裂缝单元状态进行更新，并在更新完后再次进入到预测下一时间步的缝宽分[0007]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压;6Dnn为n时刻的缝宽估计值；[0008]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压[0009]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压;;[0010]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压所述建立分布式光纤应变的观测方程，计算下一时间步状态变量对应的光纤应变的步骤;7[0011]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压;n+1为n+1时刻的观测方程的系数矩阵；[0012]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压;n+1为n+1时刻监测到的光纤应变数据；Dnrn为步骤S300中预测出的n+1时刻的缝宽估计值；;n+1为n+1时刻的观测方程的系数矩阵；Ie为应变的测量误差的协方差矩阵；;[0013]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压[0014]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压8[0015]优选地，在所述融合光纤应变监测和裂缝扩展模拟的压[0019]图1为本发明融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法的流程91示意出了本发明融合光纤应变监测和裂缝扩展模拟的压裂裂缝动态反演方法的流程图。[0026]步骤S100处初始化状态变量参数。其中状态变量参数包括地层参数和流体参;；;;;3x2323。;Dnn为n时刻的缝宽估计值；;;;[0039]根据上述分布式光纤应变的观测方程，计算下一时间步状态变量对应的光纤应;n+1为n+1时刻的观测方程的系数矩阵；n+1为n+1时刻监测到的光纤应变数据；Dnrn为步骤S300中预测出的n+1时刻的缝宽估计值；;n+1为n+1时刻的观测方程的系数矩阵；Ie为应变的测量误差的协方差矩阵。;[0045]每个时间步均通过模型误差和测量误差的最小化得到状态参数(裂缝宽度)的最[0048]主要地层参数，弹性模量E为35GPa，泊松比v为0.2，临界应力