专利摘要

一种适用于欠平衡钻井的地层特性随钻解释方法，其特征在于，所述方法利用井底 压力和液相出口流量作为测量参数进行地层特性随钻解释，解释参数的确定过程包括： 步骤1，获取模拟井的基本参数，构建油气水三相变质量流动模型； 步骤2，构建井筒水力学系统的状态方程与测量方程，所述状态方程根据所述解释参数 建立，所述测量方程根据油气水三相变质量流动模型建立； 步骤3，根据所述状态方程与测量方程，基于无迹卡尔曼滤波的解释算法进行所述解释 参数的预测、更新和校正，得到所述解释参数的值； 所述步骤2中所述状态方程和测量方程中，以向量形式的状态向量表示所述解释参数， 以向量形式的测量向量表示所述测量参数； 所述状态方程表示相邻时刻点的所述状态向量的关系，所述测量方程表示的同一时刻 点所述状态向量和测量向量的关系； 所述步骤1中获取模拟井的基本参数后，确定所述油气水三相变质量流动模型的过程 为： 步骤101，建立油气水三相变质量流动模型的控制方程，所述控制方程包括注入气质量 守恒方程、钻井液质量守恒方程、产出气质量守恒方程、产出油质量守恒方程和油气水三相 动量守恒方程； 步骤102，将所述控制方程进行离散化处理； 步骤103，建立油气水三相变质量流动的辅助方程，所述辅助方程包括：漂移流方程、PR 状态方程以及地层渗流方程； 步骤104，针对所述油气水三相变质量流动模型进行数值计算，所述数值计算的过程为 通过所述状态向量得到所述测量向量的过程； 所述步骤102中，采用四点差分显示格式对所述控制方程进行离散化处理： 轴向空间网格划分为定步长；时间网格划分为两个阶段：在气体运移到井口之前，基于 气液界面追踪得到实时变化的时间步长；在气体运移到井口之后的拟稳态过程中，时间步 长设为定值； 所述控制方程中的一阶空间导数采用一阶迎风格式，一阶时间导数采用四点中心差 分； 所述步骤104包括： k+1 步骤10401，估算k+1时刻井底压力P0 ，带入前一时刻模拟计算中给定的地层压力pr和 地层渗透率K； 步骤10402，结合所述地层渗流方程，计算给定储层流体组分下的侵入量； k+1(0) 步骤10403，假设节点j处k+1时刻的压力Pj ，用所述PR状态方程求解气相和液相组 分的密度和黏度的物性参数； k+1(0) 步骤10404，假设节点j处k+1时刻的气相体积分数αj ，由质量守恒方程计算气液相 各组分速度； k+1 k+1 k+1(0) 步骤10405，结合所述漂移流方程计算新的气相体积分数αj ，判断|αj ‑αj | 否成立，是，执行步骤10406，否则返回所述步骤10404重新迭代计算，ε为根据精度要求设置 的阈值； k+1 步骤10406，将已知参数代入所述动量守恒方程，计算新的节点j处k+1时刻的压力Pj ， k+1 k+1(0) k+1(0) 若|Pj ‑Pj | 量，否则返回所述步骤10403重新迭代计算； * * * 步骤10407，节点一直循环至井口，得到新的井口回压Pc ，若|Pc‑Pc| 确，继续下一时刻的计算，否则返回所述步骤10401。 2.根据权利要求1所述的方法，其特征在于，所述状态方程和测量方程为： xk＝f(xk‑1)+wk‑1； zk＝h(xk)+vk； xk和xk‑1分别为k和k‑1时刻的状态向量，f(·)为非线性系统状态函数；zk为测量向量，h (·)为所述油气水三相变质量流动模型；wk‑1、vk分别为系统对应时刻的过程噪声和测量噪 声。 3.根据权利要求1所述的方法，其特征在于，所述步骤3包括： 步骤301，设定初始时刻的所述状态向量和误差协方差矩阵； 步骤302，根据选择的比例修正对称采样策略，由前一时刻的状态向量的估计值和误差 协方差矩阵预测当前时刻的状态向量的均值和误差协方差矩阵； 步骤303，利用测量方程传递Sigma点，预测当前时刻的测量向量的值、自协方差矩阵和 互协方差矩阵； 步骤304，利用步骤303中预测的当前时刻的测量向量的值、自协方差矩阵和互协方差 矩阵，对当前时刻的状态向量的均值和误差协方差矩阵进行校正更新。 4.根据权利要求1所述的方法，其特征在于，所述步骤302包括： 步骤30201，由前一时刻状态向量的估计值 和误差协方差矩阵Pk‑1来计算Sigma点 (m) (c) xi,k‑1、均值加权所用的权值Wi 和协方差加权所用的权值Wi ： 2 其中，λ为比例因子，λ＝α(n+κ)‑n；n为状态向量的维数；κ为可调参数，对于高斯分布， 当n≤3,取κ＝3‑n，当n>3,取κ＝0；α为扩展因子，影响Sigma点围绕x的分布状态，取值范围 为[0,1]；β为描述状态向量先验分布信息的参数，对于高斯分布，最优值为2； 步骤30202，利用所述状态方程传递Sigma点为xi,k/k‑1： xi,k/k‑1＝f(xi,k‑1) i＝0,1,…,2n； 步骤30303，预测当前时刻的状态均值 和协方差矩阵Pk/k‑1： 其中，Qk‑1为前一时刻的过程噪声协方差矩阵。 5.根据权利要求1所述的方法，其特征在于，所述步骤303中，利用测量方程传递Sigma 点为zi,k/k‑1：zi,k/k‑1＝h(xi,k‑1)i＝0,1,…,2n； 预测当前时刻的测量向量的值 自协方差矩阵 和互协方差矩阵 的公式为： 其中，Rk为当前时刻的测量噪声协方差矩阵。 6.根据权利要求1所述的方法，其特征在于，所述步骤304包括： 步骤30401，根据自协方差矩阵和互协方差矩阵计算卡尔曼增益矩阵Kk： 步骤30402，基于步骤303中获取的当前时刻的所述测量向量的值，对当前时刻的状态 向量的均值 和误差协方差矩阵Pk进行校正更新：

一种适用于欠平衡钻井的地层特性随钻解释方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。