获得电荷参数的方法、分子力学模拟结果的方法及装置与流程

1.本技术属于分子模拟技术领域，具体涉及一种获得电荷参数的方法、分子力学模拟结果的方法及装置。


背景技术：

2.建立于经典力学基础上的分子力学以其速度优势在很多领域有着广泛的应用，例如，药物设计和材料设计领域。分子力场是分子力学中不可或缺的组成部分，分子力场中参数的好坏直接决定分子力学模拟结果的精度。传统的分子力场中包含成键项(例如，键长、键角、二面角等)和非键项(例如，电荷作用、范德华作用等)等。在非共价药物分子和蛋白质分子结合的过程中，非键项尤其是电荷作用起至关重要的作用。
3.而分子力场计算电荷作用需要力场中的电荷参数，且目前普遍采用点电荷模型来计算获得电荷参数。现在一般点电荷获得的方式是基于气相，也就是单个分子在孤立的环境中拟合量子化学方法产生的空间电势获得。但是在真实的生物应用场景中，所需要模拟的过程常常在溶液(一般是水)或蛋白环境中。这些环境会影响小分子本身的电势分布，从而造成和气相中不同的空间电势。即在气相中拟合获得的电荷参数可能会在实际的复杂环境中带来不精确的结果。


技术实现要素：

4.本技术提供一种获得电荷参数的方法、分子力学模拟结果的方法及装置，通过考虑溶质分子所处的环境来拟合获得电荷参数，以提升溶质分子在所处环境中分子力学模拟结果的精确度。
5.为解决上述技术问题，本技术采用的一个技术方案是：提供一种获得电荷参数的方法，包括：s1、对溶质分子进行结构优化以获得至少一个目标分子，并获得每个所述目标分子的初始静电势电荷；s2、针对每个所述目标分子，基于所述目标分子和多个环境分子构建模拟系统；其中，所述目标分子被多个所述环境分子环绕；s3、利用所述初始静电势电荷对所述模拟系统进行分子动力学模拟，以至少获得模拟过程中多个时刻下的所述环境分子的坐标信息；s4、基于多个时刻下的所述环境分子的坐标信息以及所述环境分子中各个原子的电荷量获得所述目标分子的第一静电势电荷；s5、判断所述第一静电势电荷与所述初始静电势电荷是否符合第一预设条件；s6、若是，将所述第一静电势电荷作为所述目标分子被所述环境分子环绕时的电荷参数；s7、否则，将所述第一静电势电荷作为所述初始静电势电荷并返回至步骤s3。
6.为解决上述技术问题，本技术采用的另一个技术方案是：提供一种获得分子力学模拟结果的方法，包括：利用上述任一实施例中所述的方法获得溶质分子的电荷参数；基于所述电荷参数获得所述溶质分子和所述环境分子的分子力学模拟结果。
7.为解决上述技术问题，本技术采用的另一个技术方案是：提供一种获得电荷参数的装置，包括：第一获得模块，用于对溶质分子进行结构优化以获得至少一个所述目标分
子，并获得每个所述目标分子的初始静电势电荷；构建模块，与所述第一获得模块连接，用于针对每个所述目标分子，基于所述目标分子和多个环境分子构建模拟系统；其中，所述目标分子被多个所述环境分子环绕；第二获得模块，与所述构建模块连接，用于利用所述初始静电势电荷对所述模拟系统进行分子动力学模拟，以至少获得模拟过程中多个时刻下的所述环境分子的坐标信息；第三获得模块，与所述第二获得模块连接，用于基于多个时刻下的所述环境分子的坐标信息以及所述环境分子中各个原子的电荷量获得所述目标分子的第一静电势电荷；判断模块，与所述第三获得模块连接，用于判断所述第一静电势电荷与所述初始静电势电荷是否符合第一预设条件；执行模块，与所述判断模块连接，用于在所述判断模块判断为是时，将所述第一静电势电荷作为所述目标分子被所述环境分子环绕时的电荷参数；以及用于在所述判断模块判断为否时，将所述第一静电势电荷作为所述初始静电势电荷并返回至所述利用所述初始静电势电荷对所述模拟系统进行分子动力学模拟的步骤。
8.为解决上述技术问题，本技术采用的另一个技术方案是：提供一种电子设备，包括相互耦接的存储器和处理器，所述存储器中存储有程序指令，所述处理器用于执行所述程序指令以实现上述任一实施例中所述的获得电荷参数的方法，或所述的获得分子力学模拟结果的方法。
9.为解决上述技术问题，本技术采用的另一个技术方案是：提供一种存储装置，存储有能够被处理器运行的程序指令，所述处理器用于执行所述程序指令以实现上述任一实施例中所述的获得电荷参数的方法，或所述的获得分子力学模拟结果的方法。
10.区别于现有技术情况，本技术的有益效果是：本技术获得电荷参数的方法的过程包括：对溶质分子进行结构优化以获得至少一个目标分子，并获得每个所述目标分子的初始静电势电荷；针对每个所述目标分子，基于所述目标分子和多个环境分子构建模拟系统；其中，所述目标分子被多个所述环境分子环绕；利用所述初始静电势电荷对所述模拟系统进行分子动力学模拟，以至少获得模拟过程中多个时刻下的所述环境分子的坐标信息；基于多个时刻下的所述环境分子的坐标信息以及所述环境分子中各个原子的电荷量获得所述目标分子的第一静电势电荷；判断所述第一静电势电荷与所述初始静电势电荷是否符合第一预设条件；若是，将所述第一静电势电荷作为所述目标分子被所述环境分子环绕时的电荷参数；否则，将所述第一静电势电荷作为所述初始静电势电荷并返回至所述利用所述初始静电势电荷对所述模拟系统进行分子动力学模拟的步骤。通过上述方式，本技术通过考虑环境分子对溶质分子的电荷参数的影响，以提升溶质分子在所处环境中分子力学模拟结果的精确度。
附图说明
11.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图，其中：
12.图1为本技术获得电荷参数的方法一实施方式的流程示意图；
13.图2为图1中步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的一实施方式的流程示意图；
14.图3为图1中步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的一实施方式的流程示意图；
15.图4为图1中步骤s3对应的一实施方式的流程示意图；
16.图5为图1中步骤s4对应的一实施方式的流程示意图；
17.图6为本技术获得分子力学模拟结果的方法一实施方式的流程示意图；
18.图7为本技术获得电荷参数的装置一实施方式的结构示意图；
19.图8为本技术电子设备一实施方式的结构示意图；
20.图9为本技术存储装置一实施方式的结构示意图。
具体实施方式
21.下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本技术的一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本技术保护的范围。
22.请参阅图1，图1为本技术获得电荷参数的方法一实施方式的流程示意图，该方法包括：
23.s1：对溶质分子进行结构优化以获得至少一个目标分子，并获得每个目标分子的初始静电势电荷。
24.具体地，在一个实施方式中，请参阅图2，图2为图1中步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的一实施方式的流程示意图。上述步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的过程具体包括：
25.s201：基于多个环境分子构建隐式溶剂模型。
26.可选地，环境分子可以包括有机溶剂分子或无机溶剂分子或蛋白质分子等；其中，有机溶剂分子包括甲醇、乙醇、苯、乙醚、二氯甲烷、丙酮、四氢呋喃、乙酸乙酯、乙腈或甲苯等；无机溶剂分子包括水、氨、氟化氢、或二氧化硫等。上述仅仅只是列举了几种比较常见的溶剂分子，在其他实施例中，可以根据实际应用场景选择对应的环境分子；且在实际应用场景中环境分子一般择一引用，即环境分子可以为上述所列举的分子中一个。另一可选地，上述隐式溶剂模型可以为onsager(昂萨格模型)、pcm(polarizable continuum model，极性连续介质模型)、cpcm(conductor-like polarizable continuum model，导体极性连续介质模型)、ipcm(isodensity polarizable continuum model，等密度极性连续介质模型)、scipcm(self-consistent isodensity polarizable continuum model，自洽等密度极性连续介质模型)、cosmo(conductor-like screening model，类导体遮蔽模型)、smd(solvation model based on density，基于密度的溶剂模型)等。
27.简单来说，隐式溶剂模型是将溶剂分子(即环境分子)当成一种连续介质的模型，即不具体描述溶质分子附近的环境分子的具体结构和分布，而是把周围环境简单地当成可极化的连续介质来考虑的。这种考虑溶剂效应的好处是可以表现溶剂的平均效应而不需要像显式溶剂模型那样需要考虑各种可能的溶剂层分子的排布方式，而且不至于令计算耗时增加很高。
28.s202：将溶质分子放置于隐式溶剂模型中进行结构优化，获得局部能量最低的目
标分子。
29.可选地，可以利用开源的量子化学计算软件(如，nwchem、gamess、gaussian等)中的hartree-fock(hf，哈特里-福克)方法进行结构优化，以获得局部能量最低的目标分子。以nwchem软件为例，在操作nwchem软件时可以直接选取关键词“task hf energy”即可获得目标分子。
30.进一步，在获得目标分子之后，还可以利用当前的量子化学计算软件获得目标分子的初始静电势(esp)电荷q0、以及目标分子中各个原子的坐标。以nwchem软件为例，在操作nwchem软件时可以直接选取关键词“task esp”即可获得目标分子的初始静电势电荷q0。
31.在又一个实施方式中，当溶质分子的分子结构较大时，其可能包括多种构象，此时可以针对不同的构象获得对应的目标分子，以提高后续获得的溶质分子的电荷参数的精确度。具体而言，请参阅图3，图3为图1中步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的一实施方式的流程示意图。上述步骤s1中对溶质分子进行结构优化以获得至少一个目标分子的过程具体包括：
32.s301：对溶质分子进行二面角旋转以获得多个不同的构象分子。
33.具体地，由公知常识可知，构象的改变不涉及共价键的改变，不同的构象分子的差异仅仅在于立体结构上的不同。
34.s302：对每个构象分子进行结构优化以获得对应的优化结构分子。
35.具体地，上述步骤s302的实现过程可以为：首先，基于多个环境分子构建隐式溶剂模型；该步骤与上述步骤s201类似，在此不再详述。接着，将构象分子放置于隐式溶剂模型中进行结构优化，获得局部能量最低的目标分子；该步骤与上述步骤s202类似，在此不再详述。
36.s303：将所有优化结构分子中的相同立体结构分子进行合并，以获得至少一个目标分子。
37.可选地，所有构象分子对应的优化结构分子的立体结构可能都相同，此时经过合并最终获得的目标分子的个数为一。例如，溶质分子a进行二面角旋转后形成的构象分子包括a、a1、a2、a3，构象分子a、a1、a2、a3对应的优化结构分子均为b，则此时目标分子为b。
38.或者，所有构象分子对应的优化结构分子的立体结构可能包括至少两种，此时经过合并最终获得的目标分子的个数为至少两种。例如，溶质分子a进行二面角旋转后形成的构象分子包括a、a1、a2、a3，构象分子a、a1对应的优化结构分子均为b，构象分子a2、a3对应的优化结构分子均为b1，则此时目标分子为b和b1。
39.s2：针对每个目标分子，基于目标分子和多个环境分子构建模拟系统；其中，目标分子被多个环境分子环绕。
40.具体地，在一个实施方式中，上述步骤s2中基于目标分子和多个环境分子构建模拟系统的步骤包括：基于多个环境分子构建显式溶剂模型；将目标分子放置于显式溶剂模型中，获得模拟系统。可选地，上述显式溶剂模型可以为tip3p(transferable intermolecular potential 3point，可迁移的分子间势的3点模型)、tip4p(transferable intermolecular potential 4point，可迁移的分子间势的4点模型)、tip5p(transferable intermolecular potential 5point，可迁移的分子间势的5点模型)等。显式溶剂模型是将每个溶剂分子(即环境分子)作为单独的分子进行模拟的模型；即显式溶剂模型中每个环境
分子都是一个单独的分子，利用显式溶剂模型可以考虑环境分子中每个原子对溶质分子的电荷参数的影响，以提高后续获得的电荷参数的准确性。
41.s3：利用初始静电势电荷对模拟系统进行分子动力学模拟，以至少获得模拟过程中多个时刻下的环境分子的坐标信息。
42.具体地，请参阅图4，图4为图1中步骤s3对应的一实施方式的流程示意图，上述步骤s3具体包括：
43.s401：利用初始静电势电荷初始化分子动力学模型。
44.具体地，可以利用分子动力学软件gromacs中内置的动力学模型进行动力学模拟；且此时动力学模型中的电荷参数为初始静电势电荷q0。
45.s402：基于初始化后的分子动力学模型对模拟系统进行动力学模拟，且在模拟过程中保持目标分子中各个原子的坐标不变，仅允许环境分子的各个原子的坐标改变。
46.可选地，在本实施例中，目标分子中各个原子的坐标可以与经过步骤s1中获得的目标分子的坐标相同；即可以保持步骤s1和步骤s3中目标分子的位置不变。
47.s403：在预定模拟时间范围内，每间隔预定时间间隔保存当前时刻下所有环境分子中各个原子的坐标。
48.可选地，预定模拟时间范围可以为2纳秒，预定时间间隔可以为10皮秒，此时可以获得200个时刻下所有环境分子中各个原子的坐标。
49.上述保持目标分子中各个原子坐标不变的方式可以使得后续步骤s4中获得第一静电势电荷的方式较为简单，且计算量较小。
50.s4：基于多个时刻下的环境分子的坐标信息以及环境分子中各个原子的电荷量获得目标分子的第一静电势电荷。
51.可选地，当步骤s3中保持目标分子的各个原子的坐标不变时，请参阅图5，图5为图1中步骤s4对应的一实施方式的流程示意图。上述步骤s4具体包括：
52.s501：将各个时刻下的所有环境分子中各个原子的坐标作为目标分子外部的多个点电荷的坐标，且将各个时刻下的所有环境分子中各个原子的电荷量与多个时刻的个数的第二比值作为对应位置处的点电荷的电荷量。
53.例如，当前溶质分子外围有100个环境分子，环境分子中每个原子的电荷量可以相同或不同。例如，当环境分子为水分子时，每个水分子有3个原子，水分子中氧原子的电荷量为-0.8，水分子中氢原子的电荷量为0.4。那么同一时刻下当前环境中有3*100个点电荷。当步骤s3中获得200个时刻下的环境分子坐标时，则当前溶剂分子外围共用300*200＝60000个点电荷，每个点电荷位置处的电荷量为该点电荷坐标处对应的原子的电荷量除以200。
54.s502：基于目标分子外部的多个点电荷的坐标及其电荷量拟合获得第一静电势电荷。
55.可选地，可以通过量子化学hf方法的单点能(single point energy)计算并拟合得到第一静电势电荷q1。
56.另一可选地，当步骤s3中动力学模拟过程中目标分子的位置也发生移动时，则此时可以先将多个时刻下的模拟系统的位置进行移动，此时模拟系统中的目标分子和环境分子会同步移动。当各个时刻下的目标分子的位置相同时，进入步骤s501。
57.s5：判断第一静电势电荷与初始静电势电荷是否符合第一预设条件。
58.具体地，上述步骤s5的具体实现过程可以为：判断第一静电势电荷q1与初始静电势电荷q0之间的平均方差是否小于阈值；可选地，该阈值可以为0.01等。该设计方式计算量较小。
59.当然，在其他实施例中，上述步骤s5的具体实现过程也可为：判断第一静电势电荷q1与初始静电势电荷q0之间的标准差是否小于阈值，本技术对此不作限定。
60.s6：若是，将第一静电势电荷作为目标分子被环境分子环绕时的电荷参数。
61.s7：否则，将第一静电势电荷作为初始静电势电荷，并返回至步骤s3。
62.在上述设计方式中，通过考虑环境分子对溶质分子的电荷参数的影响，以提升溶质分子在所处环境中分子力学模拟结果的精确度。此外，一般而言，隐式溶剂模型计算速率较快，但准确率相比与显式溶剂模型差一些；显式溶剂模型比隐式溶剂模型在考虑溶质分子和溶剂分子之间的相互作用力时具有更高的精度，但是会消耗更多的计算资源，计算效率较慢。本技术中在步骤s1中先利用计算速率快但准确率没有这么高的隐式溶剂模型进行计算，然后再利用计算速率慢但准确率高的显式溶剂模型进行计算，能够提高整个方法的效果，包括效率和准确度。此外，当步骤s1中获得的目标分子的个数为至少两个时，上述步骤s7之后，还可以包括：获得所有目标分子对应的电荷参数的和值；将和值与目标分子的个数的第一比值作为溶质分子被环境分子环绕时的电荷参数。该设计方式可以进一步提高溶质分子的电荷参数的精确度。
63.请参阅图6，图6为本技术获得分子力学模拟结果的方法一实施方式的流程示意图，上述方法具体包括：
64.s601：获得溶质分子的电荷参数。
65.可选地，在本实施例中，上述步骤s601的实现过程可参见上述任一实施例中所提及的获得电荷参数的方法。
66.s602：基于电荷参数获得溶质分子和环境分子的分子力学模拟结果。
67.在一个应用场景中，当利用步骤s601中所获得电荷参数计算一些与溶剂化相关的性质如水合自由能(即分子从气相溶解于水中所需的自由能)时，以所获得的第一结果；相比利用气相方法获得的电荷参数计算获得的第二结果而言，第一结果相比第二结果更接近实验值。
68.例如，将本技术所提供的获得电荷参数的方法应用于freesolv数据集中，将力场中的电荷参数(am1-bcc模型产生)替换为经过图1中的流程产生的考虑水环境的电荷参数，而保留力场中其他参数，最后水合自由能计算值和实验值的平均方差为1.12kcal/mol。当电荷参数采用气相方法获得时，所计算获得的水合自由能与实验值的平均方差为1.51kcal/mol。较小的方差说明计算值和实验值更加符合，本技术所提供的新的电荷参数获得流程带来更精确的计算值。
69.请参阅图7，图7为本技术获得电荷参数的装置一实施方式的结构示意图，上述获得电荷参数的装置包括第一获得模块10、构建模块12、第二获得模块14、第三获得模块16、判断模块18和执行模块11。
70.其中，第一获得模块10用于对溶质分子进行结构优化以获得至少一个目标分子，并获得每个目标分子的初始静电势电荷。构建模块12与第一获得模块10连接，用于针对每个目标分子，基于目标分子和多个环境分子构建模拟系统；其中，目标分子被多个环境分子
环绕。第二获得模块14与构建模块12连接，用于利用初始静电势电荷对模拟系统进行分子动力学模拟，以至少获得模拟过程中多个时刻下的环境分子的坐标信息。第三获得模块16与第二获得模块14连接，用于基于多个时刻下的环境分子的坐标信息以及环境分子中各个原子的电荷量获得目标分子的第一静电势电荷。判断模块18与第三获得模块16连接，用于判断第一静电势电荷与初始静电势电荷是否符合第一预设条件。执行模块11与判断模块18连接，用于在判断模块判断为是时，将第一静电势电荷作为目标分子被环境分子环绕时的电荷参数；以及用于在判断模块判断为否时，将第一静电势电荷作为初始静电势电荷并返回至利用初始静电势电荷对模拟系统进行分子动力学模拟的步骤。
71.可选地，上述第一获得模块10具体用于基于多个环境分子构建隐式溶剂模型；将溶质分子放置于隐式溶剂模型中进行结构优化，获得局部能量最低的目标分子。或者，第一获得模块10具体用于对溶质分子进行二面角旋转以获得多个不同的构象分子；对每个构象分子进行结构优化以获得对应的优化结构分子；将所有优化结构分子中的相同立体结构分子进行合并，以获得至少一个目标分子。
72.另一可选地，上述构建模块12具体用于基于多个环境分子构建显式溶剂模型；将目标分子放置于显式溶剂模型中，获得模拟系统。
73.另一可选地，上述第二获得模块14具体用于利用初始静电势电荷初始化分子动力学模型；基于初始化后的分子动力学模型对模拟系统进行动力学模拟，且在模拟过程中保持目标分子中各个原子的坐标不变，仅允许环境分子的各个原子的坐标改变；在预定模拟时间范围内，每间隔预定时间间隔保存当前时刻下所有环境分子中各个原子的坐标。
74.另一可选地，上述第三获得模块16具体用于将各个时刻下的所有环境分子中各个原子的坐标作为目标分子外部的多个点电荷的坐标，且将各个时刻下的所有环境分子中各个原子的电荷量与多个时刻的个数的第二比值作为对应位置处的点电荷的电荷量；基于目标分子外部的多个点电荷的坐标及其电荷量拟合获得第一静电势电荷。
75.另一可选地，上述判断模块18具体用于判断第一静电势电荷与初始静电势电荷之间的平均方差是否小于阈值。
76.此外，当目标分子的个数为至少两个时，上述装置还可以包括第四获得模块，与执行模块18连接，用于获得所有目标分子对应的电荷参数的和值；将和值与目标分子的个数的第一比值作为溶质分子被环境分子环绕时的电荷参数。
77.请参阅图8，图8为本技术电子设备一实施方式的结构示意图，该电子设备包括：相互耦接的存储器22和处理器20，存储器22中存储有程序指令，处理器20用于执行程序指令以实现上述任一获得电荷参数的方法，或获得分子力学模拟结果的方法。具体地，电子设备包括但不限于：台式计算机、笔记本电脑、平板电脑、服务器等，在此不做限定。此外，处理器20还可以称为cpu(center processing unit,中央处理单元)。处理器20可能是一种集成电路芯片，具有信号处理能力。处理器20还可以是、通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit,asic)、现场可编程门阵列(field-programmable gate array,fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。另外，处理器20可以由集成电路芯片共同实现。
78.请参阅图9，图9为本技术存储装置一实施方式的结构示意图，该存储装置30存储
有能够被处理器运行的程序指令300，程序指令300用于实现上述任一获得电荷参数的方法，或获得分子力学模拟结果的方法。
79.在本技术所提供的几个实施例中，应该理解到，所揭露的方法和装置，可以通过其它的方式实现。例如，以上所描述的装置实施方式仅仅是示意性的，例如，模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性、机械或其它的形式。
80.作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施方式方案的目的。
81.另外，在本技术各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
82.集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(processor)执行本技术各个实施方式方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
83.以上所述仅为本技术的实施例，并非因此限制本技术的专利范围，凡是利用本技术说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本技术的专利保护范围内。