TOP
|
特許
|
意匠
|
商標
特許ウォッチ
Twitter
他の特許を見る
公開番号
2024122933
公報種別
公開特許公報(A)
公開日
2024-09-09
出願番号
2024028503
出願日
2024-02-28
発明の名称
複雑な微細構造におけるマルチスケール反応流
出願人
ダッソー システムズ アメリカス コーポレイション
代理人
弁理士法人谷・阿部特許事務所
主分類
G16C
60/00 20190101AFI20240902BHJP(特定の用途分野に特に適合した情報通信技術)
要約
【課題】複雑な微細構造におけるマルチスケール反応流を提供する。
【解決手段】実施形態は反応流システムの挙動を決定する。一つのこうした実施形態は、反応流システムの複数のモデルを定義し、各定義されたモデルは、それぞれのスケールで反応流システムを表す。反応流システムの速度場は、定義された複数のモデルのうちの(第一のそれぞれのスケールでの)第一のモデルを使用して決定され、反応流システムの拡散率は、定義された複数のモデルのうちの(第二のそれぞれのスケールでの)第二のモデルを使用して決定される。次いで、反応流システムのための複数の反応パラメータが定義される。次に、反応流システムの挙動は、決定された速度場、決定された拡散率、および定義された複数の反応パラメータを反応輸送解決器への入力として使用することによって自動的に決定される。
【選択図】図1
特許請求の範囲
【請求項1】
反応流システムの挙動を決定するコンピュータ実装方法であって、
反応流システムの複数のモデルを定義することであって、各定義されたモデルが、それぞれのスケールで前記反応流システムを表す、定義することと、
前記定義された複数のモデルのうちの、第一のそれぞれのスケールでの、第一のモデルを使用して、前記反応流システムの速度場を決定することと、
前記定義された複数のモデルのうちの、第二のそれぞれのスケールでの、第二のモデルを使用して、前記反応流システムの拡散率を決定することであって、速度場を前記決定すること、および拡散率を決定することが、一つ以上のデジタルプロセッサによって自動的に実施される、決定することと、
前記反応流システムのための複数の反応パラメータを定義することと、
前記決定された速度場、前記決定された拡散率、および前記定義された複数の反応パラメータを反応輸送解決器への入力として使用することによって、前記反応流システムの前記挙動を自動的に決定することと、を含む方法。
続きを表示(約 960 文字)
【請求項2】
所与のスケールがマイクロスケール、分子スケール、または地下スケールである、請求項1に記載の方法。
【請求項3】
前記定義された複数のモデルのうちの少なくとも一つのモデルが、前記反応流システムの特性を示す幾何学的モデルである、請求項1に記載の方法。
【請求項4】
前記反応流システムの前記複数のモデルの所与のモデルを定義することが、
一つ以上の不均一な表面反応のモデルを定義することと、
前記一つ以上の不均一な表面反応の前記定義されたモデルについて、速度則を鉱物溶解および沈殿の関数としてモデル化することと、
前記モデル化された速度則および一つ以上の均一系バルク反応のモデルに基づいて、所与のモデルを定義することと、を含む、請求項1に記載の方法。
【請求項5】
前記第一のモデルを使用して前記反応流システムの前記速度場を前記決定することが、
前記反応流システム内の材料の画像を受信することと、
前記画像を複数の相にセグメント化することであって、各相が材料、固体、または流体を表す、セグメント化することと、
前記第一のモデルを使用して、前記複数の相に基づいて、前記反応流システムの前記速度場を決定することであって、前記第一のモデルが単相流体流モデルである、決定することと、を含む、請求項1に記載の方法。
【請求項6】
前記材料が多孔性である、請求項5に記載の方法。
【請求項7】
前記材料が一つ以上の破砕をさらに含む、請求項5に記載の方法。
【請求項8】
前記材料がナノ多孔性粘土材料である、請求項5のいずれか一項に記載の方法。
【請求項9】
前記材料が、Cu
0.66
[Al
3.33
Mg
0.66
][Si
8
]O
20
[OH]
4
の式を有するワイオミング産モンモリロナイトである、請求項5に記載の方法。
【請求項10】
前記決定された速度場が多相速度場である、請求項1に記載の方法。
(【請求項11】以降は省略されています)
発明の詳細な説明
【技術分野】
【0001】
本出願は、複雑な微細構造におけるマルチスケール反応流に関する。
続きを表示(約 14,000 文字)
【背景技術】
【0002】
(関連出願)
本出願は、2023年2月28日出願の米国仮特許出願第63/487,509号の利益を主張する。上記出願の全教示は、参照により本明細書に組み込まれる。
【0003】
反応流モデリングは、流体流および化学反応を伴うプロセスを分析するためのツールである。他の実施例の中でも特に、こうしたプロセスには、地中炭素回収、地下水浄化、CO
2
注入による原油増進回収(EOR)、構造材料の浸出・腐食・汚損による原位置採掘、リサイクル、持続可能なエネルギー用途のための電気化学システム(電池、燃料電池、電解槽など)が含まれる。
【先行技術文献】
【非特許文献】
【0004】
Abbott, J.W., Hanke, F., 2022. Kinetically Corrected Monte Carlo-Molecular Dynamics Simulations of Solid Electrolyte Interphase Growth. J. Chem. Theory Comput. 18, 925-934. https://doi.org/10.1021/acs.jctc.1c00921
Abd, A.S., Abushaikha, A.S., 2021 Reactive transport in porous media: a review of recent mathematical efforts in modeling geochemical reactions in petroleum subsurface reservoirs. SN Appl. Sci. 3, 1-28. https://doi.org/10.1007/s42452-021-04396-9
Agmon, N., 1995 CHEMICAL PHYSICS The Grotthuss mechanism. Chem. Phys. Lett. 50, 456-462.
Akkermans, R. L. C., Spenley, N. A., Robertson, S. H. 2021. COMPASS III: automated fitting workflows and extension to ionic liquids. Molecular Simulation, 47, 540-551
Al-Khulaifi, Y., Lin, Q., Blunt, M.J., Bijeljic, B., 2017. Reaction Rates in Chemically Heterogeneous Rock: Coupled Impact of Structure and Flow Properties Studied by X-ray Microtomography. Environ. Sci. Technol. 51, 4108-4116. https://doi.org/10.1021/acs.est.6b06224
Allen, M. P., Tildesley, D. J. 1987. Computer simulate of Liquids;Clarendon Press,Oxford Science Publications
Beckstein, O., Sansom, M.S.P., 2003. Liquid-vapor oscillations of water in hydrophobic nanopores. Proc. Natl. Acad. Sci. U. S. A. 100, 7063-7068. https://doi.org/10.1073/pnas.1136844100
Boek, E.S., 2014. Molecular dynamics simulations of interlayer structure and mobility in hydrated Li-, Na- and K-montmorillonite clays. Mol. Phys. 112, 1472-1483. https://doi.org/10.1080/00268976.2014.907630
Casewit, C.J., Colwell, K.S., Rappe, A.K., 1992. Application of a Universal Force Field to Main Group Compounds. J. Am. Chem. Soc. 114, 10046-10053. https://doi.org/10.1021/ja00051a042
Chatterjee A.,Ebina T.,Mizukami F.,2005. Effects of Water on the Structure and Bonding of Resorcinol in the Interlayer of Montmorillonite Nanocomposite: A Periodic First Principle Study. Phys. Chem. B, 109, 15, 7306-7313 https://doi.org/10.1021/jp045775z
Chen, L., He, Y.L., Tao, W.Q., Zelenay, P., Mukundan, R., Kang, Q., 2017. Pore-scale study of multiphase reactive transport in fibrous electrodes of vanadium redox flow batteries. Electrochim. Acta 248, 425-439. https://doi.org/10.1016/j.electacta.2017.07.086
Cygan, R.T., Liang, J.J., Kalinichev, A.G., 2004. Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field. J. Phys. Chem. B 108, 1255-1266. https://doi.org/10.1021/jp0363287
Damiani, L.H., Kosakowski, G., Glaus, M.A., Churakov, S. V., 2020. A framework for reactive transport modeling using FEniCS-Reaktoro: governing equations and benchmarking results. Comput. Geosci. 24, 1071-1085. https://doi.org/10.1007/s10596-019-09919-3
Chemley, J.J., Dahowski,R.T., Davidson,C.L., Wise, M.A., Gupta, N., Kim, S.H., Malone, E.L., 2006. Carbon dioxide capture and geologic storage-A core element of a global energy technology strategy to address climate change. Jt. Glob. Chang. Res. Inst. 1-37.
Essmann, U., Perera, L., Berkowitz, M.L., Darden, T., Lee, H., Pedersen, L.G., 1995. A smooth particle mesh Ewald method. J. Chem. Phys. 103, 8577-8593. https://doi.org/10.1063/1.470117
Fager A., Crouse B., Sun G., Xu R., Freed D., 2019 “ Evaluation of Directly Simulated WAG Hysteresis at Pore Scale and its Effect on Injectivity Index” SPE 195734, SPE Offshore Europe Conf. & Exhib., Aberdeen, UK
Fager A., Otomo H., Salazar-Tio R., Balasubramanian G., Crouse B., Zhang R., Chen H., Schembre-McCabe J., 2021. “Multi-scale Digital Rock: Application of a multi-scale multi-phase workflow to a Carbonate reservoir rock”, Int. Symp. Soc. Core Analysts
Freed, D. 1998 “Lattice-Boltzmann Method for Macroscopic Porous Media Modeling.” International Journal of Modern Physics C 09 (08): 1491-1503. https://doi.org/10.1142/S0129183198001357
Holmboe, M., Bourg, I.C., 2014. Molecular dynamics simulations of water and sodium diffusion in smectite interlayer nanopores as a function of pore size and temperature. J. Phys. Chem. C 118, 1001-1013. https://doi.org/10.1021/jp408884g
Jerauld, G., Fredrich J., Lane N., Sheng Q., Crouse B., Freed D., Fager A., and Xu R. 2017. “Validation of a Workflow for Digitally Measuring Relative Permeability.” In. Society of Petroleum Engineers. https://doi.org/10.2118/188688-MS
Kang, Q., Lichtner, P.C., Zhang, D., 2007. An improved lattice Boltzmann model for multicomponent reactive transport in porous media at the pore scale. Water Resour. Res. 43, 1-12. https://doi.org/10.1029/2006WR005551
Kosakowski, G., Churakov, S. V., Thoenen, T., 2008. Diffusion of Na and Cs in montmorillonite. Clays Clay Miner. 56, 190-206. https://doi.org/10.1346/CCMN.2008.0560205
Krevor, S., Blunt, M.J., Benson, S.M., Pentland, C.H., Reynolds, C., Al-Menhali, A., Niu, B., 2015. Capillary trapping for geologic carbon dioxide storage - From pore scale physics to field scale implications. Int. J.Greenh. Gas Control 40, 221-237. https://doi.org/10.1016/j.ijggc.2015.04.006
Krynicki, K., Green, C.D., Sawyer, D.W., 1978. Pressure and temperature dependence of self-diffusion in water. Faraday Discuss. Chem. Soc. 66, 199-208. https://doi.org/10.1039/DC9786600199
Li, L., Peters,C.A., Celia, M.A., 2006. Upscaling geochemical reaction rates using pore-scale network modeling. Adv. Water Resour. 29, 1351-1370. https://doi.org/10.1016/j.advwatres.2005.10.011
Mackay, E.J., 2013. Modelling the injectivity, migration and trapping of CO2 in carbon capture and storage (CCS), Geological Storage of Carbon Dioxide (co2): Geoscience, Technologies, Environmental Aspects and Legal Frameworks. Woodhead Publishing Limited. https://doi.org/10.1533/9780857097279.1.45
MacQuarrie, K.T.B., Mayer, K.U., 2005. Reactive transport modeling in fractured rock: A state-of-the-science review. earth-Science Rev. 72, 189-227. https://doi.org/10.1016/j.earscirev.2005.07.003
Maes, J., Menke, H.P., 2021. GeoChemFoam: Direct Modelling of Multiphase Reactive Transport in Real Pore Geometries with Equilibrium Reactions. Transp. Porous Media 139, 271-299. https://doi.org/10.1007/s11242-021-01661-8
Molins, S., Soulaine, C., Prasianakis, N.I., Abbasi, A., Poncet, P., Ladd, A.J.C., Starchenko, V., Roman, S., Trebotich, D., Tchelepi, H.A., Steefel, C.I., 2021. Simulation of mineral dissolution at the pore scale with evolving fluid-solid interfaces: review of approaches and benchmark problem set. Comput. Geosci. 25, 1285-1318. https://doi.org/10.1007/s10596-019-09903-x
Noiriel, C., 2015. Resolving Time-dependent Evolution of Pore-Scale Structure, Permeability and Reactivity using X-ray Microtomography. Rev. Mineral. Geochemistry 80, 247-285.
Noiriel, C., Daval, D., 2017. Pore-Scale Geochemical Reactivity Associated with CO2 Storage: New Frontiers at the Fluid-Solid Interface. Acc. Chem. Res. 50, 759-768. https://doi.org/10.1021/acs.accounts.7b00019
Norgate, T.E., Jahanshahi, S., Rankin, W.J., 2007. Assessing the environmental impact of metal production processes. J. Clean. Prod. 15, 838-848. https://doi.org/10.1016/j.jclepro.2006.06.018
Nose, S. 1984. A unified formulation of the constant temperature molecular-dynamics methods. Journal of Chemical Physics. 81 (1): 511-519. Bibcode:1984JChPh..81..511N. doi:10.1063/1.447334.)
Oliveira, T.D.S., Blunt, M.J., Bijeljic, B., 2019. Modelling of multispecies reactive transport on pore-space images. Adv. Water Resour. 127, 192-208. https://doi.org/10.1016/j.advwatres.2019.03.012
Parkhurst, D.L., Wissmeier, L., 2015. PhreeqcRM:A reaction module for transport simulators based on the geochemical model PHREEQC. Adv. Water Resour. 83, 176-189. https://doi.org/10.1016/j.advwatres.2015.06.001
Rahromostaqim, M., Sahimi, M., 2020. Molecular Dynamics Study of the Effect of Layer Charge and Interlayer Cations on Swelling of Mixed-Layer Chlorite-Montmorillonite Clays. J. Phys. Chem. C 124, 2553-2561. https://doi.org/10.1021/acs.jpcc.9b10919
Rigby, D., 2004. Fluid density predictions using the COMPASS force field. Fluid Phase Equilib. 217, 77-87. https://doi.org/10.1016/j.fluid.2003.08.019
Sinclair, L., Thompson, J., 2015. In situ leaching of copper: Challenges and future prospects. Hydrometallurgy 306-324.
Smith, D.E., 1998. Molecular computer simulations of the swelling properties and interlayer structure of cesium montmorillonite. Langmuir 14, 5959-5967. https://doi.org/10.1021/la980015z
Sposito, G., Skipper, N.T., Sutton, R., Park, S.H., Soper, A.K., Greathouse, J.A., 1999. Surface geochemistry of the clay minerals. Proc. Natl. Acad. Sci. U. S. A. 96, 3358-3364. https://doi.org/10.1073/pnas.96.7.3358
Sridhar, N., 2017. Local corrosion chemistry - A review. Corrosion 73, 18-30. https://doi.org/10.5006/2246
Steefel, C.I., Appelo, C.A.J., Arora, B., Jacques, D., Kalbacher, T., Kolditz, O., Lagneau, V., Lichtner, P.C., Mayer, K.U., Meeussen, J.C.L., Molins, S., Moulton, D., Shao, H., Simunek, J., Spycher, N., Yabusaki, S.B., Yeh, G.T., 2015. Reactive transport codes for subsurface environmental simulation, Computational Geosciences. https://doi.org/10.1007/s10596-014-9443-x
Steefel, C.I., DePaolo, D.J., Lichtner, P.C., 2005. Reactive transport modeling: An essential tool and a new research approach for the Earth sciences. Earth Planet. Sci. Lett. 240, 539-558. https://doi.org/10.1016/j.epsl.2005.09.017
Sun G., Sun Z., Fager A., Crouse B., 2022. “Pore-scale Analysis of CO2-brine displacement in Bereatone and its Implications to CO2 Injectivity”, Int. Symp. Soc. Core Analysts
Teich-McGoldrick, S.L., Greathouse, J.A., Jove-Colon, C.F., Cygan, R.T., 2015. Swelling Properties of Montmorillonite and Beidellite Clay Minerals from Molecular Simulation: Comparison of Temperature, Interlayer Cation, and Charge Location Effects. J. Phys. Chem. C 119, 20880-20891. https://doi.org/10.1021/acs.jpcc.5b03253
Wu, Y., Tepper, H.L., Voth, G.A., 2006. Flexible simple point-charge water model with improved liquid-state properties. J. Chem. Phys. 124 https://doi.org/10.1063/1.2136877
Xu, T., Sonnenthal, E., Spycher, N., Pruess, K., 2006. TOUGHREACT - A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration. Comput. Geosci. 32, 145-165. https://doi.org/10.1016/j.cageo.2005.06.014
Xu, R., Crouse B., Freed D., Fager A., Jerauld G., Lane N., and Sheng Q.. 2018. “Continuous vs Discontinuous Capillary Desaturation and Implications for IOR/EOR”, SCA 2018-066, Int. Symp. Soc. Core Analysts, Trondheim, Norway
Yoon, H., Kang, Q., Valocchi, A.J., 2015. Lattice Boltzmann-based approaches for pore-scale reactive transport. Pore Scale Geochemical Process. 80, 393-431.https://doi.org/10.2138/rmg.2015.80.12
【発明の概要】
【0005】
反応流モデリングツールは存在するが、これらのツールは改善から利益を享受することができる。実施形態は、こうした改善された機能を提供する。
【0006】
反応流モデリングには多くの業界において用途があるが、この分野での進歩の大部分は、地下システムにおける地球化学モデリングの開発を通して行われてきた。それ故に、本明細書に開示の例示的な実施形態は、このクラスで選択された代表的な用途、すなわち銅の原位置浸出の反応流モデリングを例示する。しかしながら、実施形態は、銅(金属全般)の原位置浸出のシミュレーション/モデリングに限定されず、それよりむしろ任意の反応流システムを分析するために実施形態を使用することができることに留意されたい。原位置浸出は、標的鉱物を溶解して関心の金属(この実施例において銅)を抽出または生成するための酸の地下循環を伴う。反応流モデリングを通して流体および流体-岩石界面の組成物のリアルタイムでの変化をモデル化することは、原位置浸出による銅生成の最適化、および代替的な設計および動作条件の考えられる結果の予測を可能にし、プロセスの最適化を容易にする。
【0007】
反応流モデリングには、関心の物理学のスケールに応じて、(i)連続レベルモデル
42,47
、および(ii)細孔スケールモデル
28,29
という、モデルの二つの広範なカテゴリーが存在する。細孔スケールで微細構造は完全に分解され、体積当たりの表面反応と連動した有効種輸送は、微細構造内で直接シミュレーションされ、入力は材料の特性である。連続スケールで微細構造は分解されず、モデルは有効な輸送、微細構造によって制限される反応速度、材料空隙率、または材料の表面対体積比などの特定の特性の入力を必要とする。実際の細孔スケール構造への入力としてのX線マイクロトモグラフィ(microCT)などの三次元(3D)顕微鏡撮像法の使用が可能であるが、反応流モデリングにおいて主流ではない。現実世界の用途では、実験測定を必要とせずに、または単純化されたモデルの仮説を制限することなく、正確な結果を生成するために、両方の側面、スケール、および3D撮像をモデルに含める必要がある。反応流に関する現在の解決策におけるこのマルチスケールの側面の欠如はまた、材料、反応、拡散パラメータに対するより小さい分子レベルのモデリングを含む。現在の反応流解決器における別の制限は、多相流(油、水、ガスなどの複数の流体流)を同時に記述することに関する。
【0008】
要約すると、既存の細孔スケールモデルは、未分解の細孔を通る流れ(例えば、マルチスケール流シミュレーション)を考慮せず、多相流を正確にモデル化することはできない。既存の連続スケールモデルは、透水性および空隙率に基づいて、簡略化された流れと輸送の方程式を考慮する。従って、フィールドスケールの問題をシミュレートするためのアップスケーリングと併せて細孔スケールで多相およびマルチスケール反応流を正確に解決するための解決策はない。文献の細孔スケールモデルと連続スケールモデルの両方は、材料、反応、拡散パラメータについてオープンソースデータベースに依存する。最後に、一次反応後の流れ溶液の二次種分化を正確に考慮することは、連続モデルにおいて大抵の場合で一般的であるが、細孔スケールモデルにおいては、ある程度制限される。「二次種分化を考慮する」とは、水化学を正確にモデル化する、または反応流モデルにおいて均一系反応を含むことも指すことに留意されたい。
【0009】
本明細書に記載の通りの反応流の用途を正確にモデル化するために、実施形態は、堅牢な多相(液体および気体の連成流)およびマルチスケール(分解細孔および未分解細孔を通る流れ)流れシミュレータ、およびマルチスケール(細孔スケールおよび連続スケール)反応シミュレータを有する反応流モデルを実施する。3D撮像ベースの微細構造モデルも、特性相関に対する単純近似を回避するために、実施形態で使用される。さらに、反応流モデリングは、材料、反応、拡散パラメータからの多くの入力を利用することに留意されたいが、これらをすべての用途のために文献およびデータリポジトリから見つけるのは容易ではない。この問題を解決するために、実施形態は、化学的特性を計算する分子モデリング構成要素の使用を通して化学的特性を得る。なおもさらに、本明細書に記載の反応流モデリングの実施形態は、反応流モデルにおける不均一系反応に加えて、均一系反応を含む。
【0010】
実施形態は、反応流システムの挙動を決定するためのコンピュータ実装方法およびシステムを対象とする。一つのこうした実施形態は、反応流システムの複数のモデルを定義し、各定義されたモデルは、それぞれのスケールで反応流システムを表す。次いで、反応流システムの速度場は、定義された複数のモデルのうちの(第一のそれぞれのスケールでの)第一のモデルを使用して決定され、反応流システムの拡散率は、定義された複数のモデルのうちの(第二のそれぞれのスケールでの)第二のモデルを使用して決定される。一実施形態において、速度場を決定すること、および拡散率を決定することは、一つ以上のデジタルプロセッサによって自動的に実施される。次に、反応流システムのための複数の反応パラメータが定義される。次に、反応流システムの挙動は、決定された速度場、決定された拡散率、および定義された複数の反応パラメータを反応輸送解決器への入力として使用することによって自動的に決定される。
(【0011】以降は省略されています)
この特許をJ-PlatPatで参照する
関連特許
個人
抗体パスポート
1か月前
株式会社リンク
心理分析システム
2か月前
花王株式会社
美容システム
6か月前
個人
支援システム、支援プログラム
5か月前
個人
自動車乗員の健康管理システム
2か月前
TOTO株式会社
排泄情報管理システム
7か月前
TOTO株式会社
生体情報管理システム
7か月前
富士通株式会社
リスクと診断
12日前
株式会社Y4.com
情報処理装置
1か月前
株式会社I.W.G
医療支援システム
1か月前
住友化学株式会社
組成提案装置
5か月前
株式会社タカゾノ
作業時間計測システム
2日前
株式会社アルメックス
遠隔医療システム
6か月前
住友化学株式会社
組成提案システム
1か月前
学校法人近畿大学
腫瘍の分類推定システム
7か月前
富士通株式会社
介入方法及び介入プログラム
2か月前
株式会社Smart119
情報処理システム
10日前
トヨタ自動車株式会社
方法
6か月前
住友化学株式会社
学習装置及び予測装置
5か月前
株式会社レイヤード
情報処理システム
7か月前
株式会社BiPSEE
情報処理装置及びプログラム
2か月前
株式会社医療情報技術研究所
医療機関受付システム
5か月前
富士通株式会社
表示制御方法及び表示制御プログラム
4か月前
富士通株式会社
機器特定方法及び機器特定プログラム
2か月前
株式会社医用工学研究所
手術時間比較システム
6か月前
株式会社レイヤード
患者管理システム及びプログラム
5か月前
日本電気株式会社
児童支援システム及び児童支援方法
4か月前
個人
プログラム、情報処理方法及び情報処理装置
6か月前
三幸電子株式会社
医薬品販売管理システム
2か月前
個人
アミデーション反応のシミュレーション方法
4か月前
ロゴスサイエンス株式会社
ヘルスケアシステムおよびその方法
1か月前
株式会社 サン・システム
手術器材セット管理システム
4か月前
株式会社ミラボ
スケジュール案内装置、及びプログラム
2か月前
アサイクル株式会社
管理方法、管理システム、管理装置
4か月前
TIS株式会社
処方箋管理システム
1か月前
株式会社ウェルナス
健康管理システム
4か月前
続きを見る
他の特許を見る