WO2008041361A1 - Dispositif de mesure et procédé de mesure utilisant un procédé de résonance magnétique nucléaire - Google Patents

Dispositif de mesure et procédé de mesure utilisant un procédé de résonance magnétique nucléaire Download PDF

Info

Publication number
WO2008041361A1
WO2008041361A1 PCT/JP2007/001055 JP2007001055W WO2008041361A1 WO 2008041361 A1 WO2008041361 A1 WO 2008041361A1 JP 2007001055 W JP2007001055 W JP 2007001055W WO 2008041361 A1 WO2008041361 A1 WO 2008041361A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic field
current
sample
coil
pulse
Prior art date
Application number
PCT/JP2007/001055
Other languages
English (en)
Japanese (ja)
Inventor
Kuniyasu Ogawa
Tomoyuki Haishi
Kohei Ito
Original Assignee
Keio University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Keio University filed Critical Keio University
Priority to JP2008537406A priority Critical patent/JP5170686B2/ja
Publication of WO2008041361A1 publication Critical patent/WO2008041361A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341—Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415—Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/3808—Magnet assemblies for single-sided MR wherein the magnet assembly is located on one side of a subject only; Magnet assemblies for inside-out MR, e.g. for MR in a borehole or in a blood vessel, or magnet assemblies for fringe-field MR
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/383—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using permanent magnets
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE

Definitions

  • the present invention relates to a measuring apparatus and a measuring method using a nuclear magnetic resonance method, and more particularly to a technique for measuring a current at a specific portion of a sample using the nuclear magnetic resonance method.
  • Non-patent Documents 1 and 2 include conventional methods for measuring the surface current distribution of a sample.
  • Non-Patent Document 1 when measuring the current distribution in the surface direction of a fuel cell, the electrodes are divided into “divided electrodes”, and the individual currents flowing through each divided electrode are insulated. A method of measuring is described.
  • Non-Patent Document 2 describes a method of measuring the strength of a magnetic field using a Hall element.
  • the Hall element is an element having a characteristic that the electric resistance of the element changes according to the magnetic field strength applied to the element.
  • this Hall element is moved closer to the fuel cell and spatially scanned to measure a spatial map of the magnetic field strength and analyze it as an inverse problem to obtain the current distribution.
  • a method has been proposed.
  • Non-patent document 1 Kazuo Onda et al., “Mechanical properties measurement of solid polymer fuel cell and analysis / measurement of current distribution”, Proceedings of the 1st 3rd Fuel Cell Symposium, 2006, p. 2 3 4-2 3 7
  • Non-Patent Document 2 Masaaki Izumi, Yuji Goto, “Research and Development on Measurement Technology and Modeling of Polymer Electrolyte Fuel Cell”, NEDO Fuel Cell ⁇ Summary of Hydrogen Technology Development Interim Report Meeting, February 1999 2 Announced on 7th, p. 3 9-4 0
  • a split electrode is incorporated. It is necessary to manufacture a fuel cell, and since it is actually measured with a device for measurement, there is a possibility that the measurement result may be different from the actual machine that does not use the divided electrode, and the reliability of the experimental data is improved There was room for. In addition, each time a new cell is designed and manufactured, the split electrodes must be designed and manufactured again, which is not practical in terms of increasing development costs.
  • the magnetic field generated by the current flowing in the electrode is measured, but the magnetic field strength is almost equal to the strength of the geomagnetism. It is a weak value.
  • the Hall element is required to have high resolution and high reproducibility.
  • the Hall element when a Hall element is used for measurement of a fuel cell, the Hall element is sensitive to temperature changes, and is installed in and around a fuel cell that generates heat, and a magnetic field is generated by the Hall element.
  • the relationship between the current or resistance value flowing through the Hall element measured at each temperature and the applied magnetic field strength can be used as a calibration curve in advance and applied to the fuel cell.
  • the temperature of the hole element itself must be measured with very high accuracy, and a very time-consuming method of calculating the magnetic field from the calibration curve must be taken.
  • a static magnetic field application unit that applies a static magnetic field to the sample
  • a small RF coil smaller than the sample which applies an oscillating magnetic field for excitation to the sample and obtains a nuclear magnetic resonance signal generated at a specific location of the sample; Calculating a difference between the frequency of the nuclear magnetic resonance signal acquired by the small RF coil and the frequency of the excitation oscillating magnetic field, and calculating a current of the specific portion of the sample from the difference;
  • a measuring device is provided.
  • a static magnetic field applying unit that applies a static magnetic field to the solid polymer electrolyte membrane; and an oscillating magnetic field for excitation applied to the solid polymer electrolyte membrane, and generated at a specific location of the solid polymer electrolyte membrane.
  • a plurality of small RF coils that are smaller than the solid polymer electrolyte membrane,
  • a difference between the frequency of the nuclear magnetic resonance signal acquired by the small RF coil and the frequency of the oscillating magnetic field for excitation is calculated, and from the difference, the solid polymer electrolyte membrane A current distribution acquisition unit for acquiring an in-plane current distribution;
  • a measuring device is provided.
  • a method of locally measuring the current at a specific part of a sample using a nuclear magnetic resonance method is a method of locally measuring the current at a specific part of a sample using a nuclear magnetic resonance method.
  • an excitation oscillating magnetic field is applied to a specific location of the sample placed in a static magnetic field using a small RF coil smaller than the sample, and a nuclear magnetic resonance signal generated at the specific location is acquired.
  • a measurement method is provided.
  • an excitation oscillating magnetic field is applied locally using a small RF coil smaller than the sample, and emitted from the location where the excitation oscillating magnetic field is applied.
  • the obtained nuclear magnetic resonance signal is obtained, and the current at a specific portion of the sample is obtained from the obtained nuclear magnetic resonance signal.
  • the measurement accuracy can be improved by using a nuclear magnetic resonance signal with high frequency resolution when obtaining the current.
  • the measurement using the Hall element described above with reference to Non-Patent Document 2 is performed. Compared with methods that use absolute values, the effects of changes in the environment around the element, such as the temperature environment, and the need for a calibration curve can be reduced, so measurement accuracy can be further improved.
  • the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the oscillating magnetic field for excitation can be specifically obtained as (i i) below with respect to (i) below.
  • the RF oscillator as an oscillating magnetic field for excitation has (remembered) “frequency that is the reference when no current is flowing”
  • the above (i i) is measured, for example, as a phase change amount.
  • the frequency of the excitation oscillating magnetic field is equal to the frequency of the nuclear magnetic resonance signal under the static magnetic field created only by the magnet, and the current flows.
  • the "static magnetic field” is a completely stable magnetic field as long as the magnetic field is stable in time so that the nuclear magnetic resonance signal and current can be stably acquired. It may not be a magnetic field, and there may be some variation within that range.
  • the obtained current is expressed as a spatial distribution. Therefore, it can be expressed as the current density divided by the area where the current flows.
  • the measurement apparatus of the present invention may further include a detection unit that detects a real part and an imaginary part of the nuclear magnetic resonance signal, and the current calculation unit is detected by the detection unit.
  • the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the excitation oscillating magnetic field may be calculated using the imaginary part and the imaginary part.
  • the real part and the imaginary part of the nuclear magnetic resonance signal are detected, and the frequency of the nuclear magnetic resonance signal and the frequency using the real part and the imaginary part are detected.
  • the difference from the frequency of the oscillating magnetic field for excitation can also be calculated.
  • the frequency difference can be obtained more simply and reliably.
  • the excitation vibration magnetic field is used as a reference.
  • the amount of change in phase at a certain time interval of the nuclear magnetic resonance signal or the difference between the two frequencies may be calculated.
  • a storage unit that stores information indicating the correlation between the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the excitation vibration magnetic field and the current, for example, calibration curve data
  • the current calculation unit may acquire the information corresponding to the sample to be measured from the storage unit, and calculate the current based on the information.
  • the sample to be measured in the present invention can be, for example, a film.
  • the measurement apparatus of the present invention includes a plurality of the small RF coils, and the plurality of small RF coils apply the excitation oscillating magnetic field to a plurality of locations of the sample, and the nuclear magnetic resonance.
  • a signal may be acquired, and the current calculation unit may be configured to calculate currents at the plurality of locations of the sample.
  • the small RF coil applies, for example, the pulsed oscillating magnetic field, and also outputs a FID (Free Inductive Decay) signal corresponding to the oscillating magnetic field for excitation.
  • the current calculation unit can acquire the real part and the imaginary part of the FID signal.
  • the angle of tilting as a reference) is arbitrary.
  • the small RF coil can apply an excitation oscillating magnetic field in the following sequence, for example, and can also acquire an echo signal corresponding to the excitation oscillating magnetic field.
  • the excitation oscillating magnetic field is a pulse sequence including the above (a) and (b), and the current calculation unit uses the echo echo method to obtain the real part and imaginary part of the echo signal.
  • the phase of the signal can be converged.
  • measurement errors due to magnetic field inhomogeneities can be effectively reduced. For this reason, the measurement accuracy of the real part and the imaginary part of the nuclear magnetic resonance signal can be further improved.
  • the real part and the imaginary part of the echo signal are detected, and the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the excitation oscillating magnetic field is calculated using the real part and the imaginary part. You may calculate it.
  • FID signal and "echo signal” are nuclear magnetic resonance signals that correspond to the excitation oscillating magnetic field and can detect the real part and the imaginary part. Any signal that functions as a signal may be used.
  • the "pulse sequence" in the present invention is a sequence that defines a timing diagram for setting a time at which an excitation oscillating magnetic field is applied and its interval.
  • the timing diagram also includes a procedure table for performing necessary operations in time series.
  • another sequence in which a step of applying a 180 ° pulse is added at a time just before the 90 ° pulse (a) may be executed.
  • the intensity of the N MR signal acquired with the 90 ° pulse (a) with the intensity of the N MR signal acquired by selecting the time at the 1 80 ° pulse (b) as appropriate, the irradiation from the RF coil It can be determined whether the intensity of the oscillating magnetic field for excitation corresponds exactly to 90 ° and 1800 °.
  • the intensity of the two pulses is 1 to 2, or the irradiation energy is 1 to 4, or the pulse application time is 1 to 2, and the magnetization vector is excited to 90 ° and 1 80 ° respectively.
  • an RF signal generation unit that generates an RF signal for generating the excitation oscillating magnetic field in the small RF coil, and an echo signal obtained by the small RF coil are detected, and An echo signal detection unit that sends an echo signal to the current calculation unit; and a branch unit that connects the small RF coil, the RF signal generation unit, and the echo signal detection unit; and the small RF coil and the RF
  • a switch circuit that switches between a state in which the signal generation unit is connected and a state in which the small RF coil and the echo signal detection unit are connected can be provided.
  • the pulse sequence of the excitation oscillating magnetic field may include the following (a), (b), and (c).
  • N is a natural number.
  • the current at a specific part of the sample is measured using an echo signal corresponding to the pulse (b) or (c), and Using a plurality of echo signals corresponding to the pulses of b) and (c), the amount of the protic solvent in the sample at the specific location can be measured.
  • the measurement apparatus of the present invention includes a solvent amount calculation unit that calculates the amount of the protonic solvent in the sample based on the nuclear magnetic resonance signal acquired by the small RF coil, A switching unit that switches between a first measurement mode for measuring current and a second measurement mode for measuring the amount of the protonic solvent in the sample, and when in the first measurement mode, the current calculation The unit calculates the current at the specific location of the sample based on the difference between the frequency of the nuclear magnetic resonance signal acquired by the small RF coil and the frequency of the oscillating magnetic field for excitation, and the second measurement mode
  • the solvent amount calculation unit may calculate the amount of the protonic solvent at the specific location in the sample based on the nuclear magnetic resonance signal acquired by the small RF coil. so That.
  • the small RF coil acquires an echo signal corresponding to the excitation oscillating magnetic field, and the solvent amount calculation unit calculates from the intensity of the echo signal. , Calculate the T 2 relaxation time constant and calculate ⁇ 2 relaxation time constant From the number, the amount of the protic solvent at a specific location in the sample can also be calculated.
  • the small RF coil applies the excitation oscillating magnetic field in a pulse sequence including (a), (b) and (c), and An FID signal corresponding to the pulse of (a) or an echo signal corresponding to the pulse of (b) or (c) is acquired, and in the second step, the FID signal corresponding to the pulse of (a) or the The difference between the frequency of the nuclear magnetic resonance signal and the frequency of the oscillating magnetic field for excitation is calculated using the real part and the imaginary part of the echo signal corresponding to the pulse of (b) or (c). ) And above, in the small RF coil applies the excitation oscillating magnetic field in a pulse sequence including (a), (b) and (c), and An FID signal corresponding to the pulse of (a) or an echo signal corresponding to the pulse of (b) or (c) is acquired, and in the second step, the FID signal corresponding to the pulse of (a) or the The difference between the frequency of the nuclear magnetic resonance signal and the frequency of the oscillating magnetic field for excitation is calculated using
  • the T 2 relaxation time constant is calculated from the intensity of a plurality of echo signals corresponding to the pulse of (c), and the proton solvent at the specific location in the sample is calculated from the calculated ⁇ ⁇ 2 relaxation time constant.
  • the amount of can also be calculated. More specifically, the above
  • the current can be measured with the first pulse-corresponding signal in (b), and the moisture content can be measured using the subsequent n-pulse-corresponding signal group in (c).
  • the measurement of the two physical quantities may be performed at the same time, or may be performed at different timings such as alternating.
  • the acquisition of an echo signal for calculating the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the oscillating magnetic field for excitation, and the echo signal for calculating the T 2 relaxation time constant Acquisition may be performed alternately.
  • the nuclear magnetic and acquisition of echo signals for calculating the difference between the frequency of the resonance signal and the frequency of the excitation oscillating magnetic field, T 2 relaxation time E co one for calculating the constants Signal acquisition may be performed simultaneously.
  • a small RF coil is used to (i) apply an excitation vibration magnetic field locally, and (ii) obtain an echo signal emitted from a location where the excitation vibration magnetic field is applied.
  • the T 2 relaxation time constant (lateral relaxation) Sum the time constant) and measure the water content based on this. Since the pulse-echo method is applied with a small RF coil to limit the area to be measured, the local water content can be measured in a short time.
  • the 90 ° pulse is in the first phase, and the n number of the 180 ° pulses are in the second phase shifted by 90 ° from the first phase. Can do.
  • magnetic field inhomogeneities of the static magnetic field and excitation oscillating magnetic field occur, which may cause measurement errors in the T 2 relaxation time constant.
  • the pulse sequence with the above configuration is 90 ° with the first phase as a 180 ° pulse. Since the one in the shifted second phase is used, by applying a 180 ° pulse, the nuclear magnetization is reversed in the rotating coordinate system, and this causes a measurement error factor due to the above magnetic field inhomogeneity. It will be resolved.
  • Second phase 1 8 0 ° pulses since it is periodically applied, so that every time the measurement error factors are eliminated, it is possible to reliably obtain accurate T 2 relaxation time constant.
  • a plurality of the small RF coils are provided, and in the second measurement mode, an excitation oscillating magnetic field is applied to a plurality of locations of the sample and an echo signal corresponding to the excitation oscillating magnetic field is acquired.
  • the moisture content distribution of the sample may be presented based on the moisture content at the plurality of locations of the sample.
  • a gradient magnetic field application unit that applies a gradient magnetic field to a sample, and a proton in the sample based on a nuclear magnetic resonance signal acquired by the small RF coil.
  • a mobility calculation unit for calculating the mobility of the organic solvent, and a first measurement mode for measuring the current of the sample and a third measurement mode for measuring the mobility of the protonic solvent in the sample.
  • a switching unit; and in the third measurement mode, the small RF coil applies the excitation vibration magnetic field to the sample, and the nuclear magnetic resonance corresponding to the excitation vibration magnetic field and the gradient magnetic field.
  • the signal is acquired, and the mobility calculation unit supports different gradient magnetic fields.
  • the mobility of the specific portion of the sample may be calculated based on the information of the nuclear magnetic resonance signal obtained as described above.
  • a magnetic resonance signal is acquired, and in the first step, application of a gradient magnetic field to the sample is executed according to a predetermined pulse sequence; in the third step, the first step It is also possible to perform the application of the gradient magnetic field of step with different sizes according to a predetermined pulse sequence
  • the small RF coil is used to (i) apply the excitation oscillating magnetic field and gradient magnetic field locally, and (ii) apply the excitation oscillating magnetic field and gradient magnetic field.
  • the nuclear magnetic resonance signal emitted from the sample is acquired, and the mobility at a specific location of the sample is measured from the NMR signal obtained corresponding to different gradient magnetic fields. Since the spin echo method and gradient magnetic field NMR method are applied by limiting the measurement target area with a small RF coil, the local mobility of the protonic solvent in a predetermined region of the sample can be reduced in a short time. It can be measured.
  • the "mobility” measured by the present invention refers to a physical property value representing the movement of a protonic solvent in a sample.
  • Such physical property values include parameters such as self-diffusion coefficient and mobility (movement speed). Main departure According to Ming, one of these parameters can be obtained.
  • different gradient magnetic fields includes a case where one gradient magnetic field is zero, that is, measurement performed without applying a gradient magnetic field.
  • the first step may be performed without applying a gradient magnetic field.
  • the gradient magnetic field application unit may take various modes. For example, it can be a gradient magnetic field application coil arranged away from a small RF coil, or a planar coil provided in the same plane as the small RF coil. Alternatively, a pair of gradient magnetic field application coils arranged with a small RF coil in between may be used. Alternatively, these configurations may be arbitrarily combined.
  • the pair of gradient magnetic field applying coils may have a substantially half-moon planar shape, and half-moon strings may be arranged to face each other toward the small RF coil.
  • a substantially half-moon shape means that a pair of planar coils have a chord-like straight line region, and a gradient magnetic field that inclines in a direction perpendicular to the straight line region by arranging them in a sample. The configuration can be applied, and if the application of such a gradient magnetic field is possible, the plane shape of the moon shape of the coil may be larger or smaller than half a month.
  • the sample may include a matrix made of a solid or a gel.
  • the solvent amount calculation unit and the mobility calculation unit can be configured to calculate the amount and mobility of the protonic solvent contained in the matrix, respectively.
  • Examples of such a sample include a film containing moisture, for example, a solid electrolyte film used for a fuel cell or the like.
  • the protonic solvent refers to a solvent that dissociates itself to generate proton.
  • Protonic solvents include, for example, water; alcohols such as methanol and ethanol;
  • Carboxylic acids such as acetic acid
  • Liquid ammonia is mentioned.
  • water and alcohol are solvents that can more easily measure the mobility or the amount of solvent in the present invention.
  • FIG. 1 is a flowchart showing a procedure for measuring current in the embodiment.
  • FIG. 2 is a diagram for explaining a compensation function of the CPMG method.
  • FIG. 3 is a flowchart showing a procedure for measuring a water content in the embodiment.
  • FIG. 4 A diagram for explaining the principle of acquiring the N MR signal by the spin echo method.
  • FIG. 5 is a diagram showing an example of a pulse sequence for self-diffusion coefficient measurement.
  • FIG. 6 is a flowchart showing the procedure for measuring the self-diffusion coefficient in the embodiment.
  • FIG. 7 is a diagram showing a schematic configuration of a measuring apparatus according to an embodiment.
  • FIG. 8 is a diagram showing an example of an LC circuit that performs application of an oscillating magnetic field for excitation and detection of an NMR signal of the measurement apparatus in the embodiment.
  • FIG. 9 is a diagram showing a configuration of a switch unit of the measuring apparatus in the embodiment.
  • FIG. 10 is a diagram for explaining a phase difference shift of an N M R signal.
  • FIG. 11 is a perspective view showing an arrangement example of a plurality of small RF coils of the measuring apparatus in the embodiment.
  • FIG. 12 is a diagram showing a configuration of an output unit of the measurement apparatus in the embodiment.
  • FIG. 13 is a diagram showing a schematic configuration of a measuring apparatus according to an embodiment.
  • FIG. 14 is a diagram showing a configuration of a control unit of the measurement apparatus in the embodiment.
  • FIG. 15 is a flowchart showing a measurement procedure in the embodiment.
  • FIG. 16 is a diagram showing a schematic configuration of a measuring apparatus according to an embodiment.
  • FIG. 17 is a diagram showing a configuration of a G coil of the measuring apparatus in the embodiment.
  • FIG. 18 is a flowchart showing a measurement procedure in the embodiment.
  • FIG. 19 is a diagram showing a schematic configuration of a measuring apparatus in an embodiment.
  • FIG. 20 is a diagram for explaining a method of measuring current in the example.
  • FIG. 21 is a view showing a sample used in an example.
  • FIG. 22 is a diagram showing an echo signal in the example.
  • FIG. 23 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 24 is a diagram showing an echo signal in the example.
  • FIG. 25 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 26 is a diagram showing an echo signal in the example.
  • FIG. 27 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 28 is a diagram showing the relationship between current and frequency shift in the example.
  • FIG. 29 is a diagram showing an echo signal in the example.
  • FIG. 30 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 31 is a diagram showing the relationship between current and frequency shift in the examples.
  • FIG. 32 is a diagram for explaining a method of measuring current in the example.
  • FIG. 33 is a diagram showing a part of a configuration of a current measuring apparatus in an example.
  • FIG. 34 is a diagram for explaining a method of measuring current in the example.
  • FIG. 35 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 36 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 37 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 38 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 39 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 40 is a diagram showing a relationship between current and frequency shift in the example.
  • FIG. 41 is a diagram showing a configuration of a measurement apparatus in the embodiment.
  • FIG. 42 is a diagram showing a time change of the F I D signal in the example.
  • FIG. 43 is a diagram showing a time change of the phase difference in the example.
  • FIG. 44 is a diagram for explaining the time change of the phase difference in the example.
  • FIG. 45 is a diagram showing the time change of the F I D signal in the example.
  • FIG. 46 is a diagram showing a time change of the phase difference in the example.
  • FIG. 47 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 48 is a diagram showing a time change of the F I D signal in the example.
  • FIG. 49 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 50 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 51 is a diagram showing a time change of the F I D signal in the example.
  • FIG. 52 is a diagram showing a temporal change in phase difference in the example.
  • FIG. 53 is a diagram showing the time change of the phase difference in the example.
  • FIG. 54 is a diagram showing the relationship between current and frequency shift in the example.
  • FIG. 55 is a diagram showing directions of a static magnetic field H 0 and a magnetic field H i used for the analysis in the example.
  • FIG. 56 is a perspective view showing the positions of the copper plate, RF detection coil, and water sample in the magnetic field analysis in the example.
  • FIG. 57 is a diagram showing a coordinate system used for magnetic field analysis in an example.
  • FIG. 58 is a diagram showing an analysis result in the example.
  • FIG. 59 is a diagram showing the relationship between the frequency shift amount ⁇ and the current I in the example.
  • FIG. 60 is a diagram showing the configuration of the small surface coil used in the example.
  • FIG. 61 is a diagram showing the relationship between the current I and the frequency shift amount ⁇ in the example.
  • FIG. 62 is a cross-sectional view showing a schematic configuration of ⁇ ⁇ in the example.
  • FIG. 63 is a diagram showing current flowing through the collector electrode and the PEM in the example.
  • FIG. 64 is a diagram showing an analysis result of a frequency shift amount ⁇ in the example.
  • FIG. 65 is a diagram showing a configuration of a small surface coil used in Examples.
  • FIG. 66 is a diagram showing a coil placed on a force-pong mesh in an example.
  • FIG. 67 is a diagram showing a measurement result of a frequency shift amount ⁇ in the example.
  • FIG. 68 is a cross-sectional view showing a schematic configuration of ⁇ ⁇ ⁇ used in the analysis in this example.
  • FIG. 69 is a diagram showing an analysis result of a frequency shift amount in the example.
  • FIG. 70 is a cross-sectional view showing a schematic configuration of M EA in an example.
  • FIG. 71 is a diagram showing a frequency shift obtained by measurement and analysis in the example.
  • FIG. 72 is a cross-sectional view showing the arrangement of M EA and small coils in the example.
  • FIG. 73 is a diagram showing the measurement timing of P G S E and C P M G in an example.
  • FIG. 74 is a diagram showing a time change of a current flowing through M E A in the P G S E measurement of the example.
  • FIG. 75 is a diagram showing the time change of the voltage applied to M E A in the P G S E measurement of the example.
  • FIG. 76 is a diagram showing the measurement result of the frequency shift amount and the analysis result of the frequency shift amount in the PGSE measurement of the example.
  • FIG. 77 is a diagram showing the time change of the frequency shift amount in the PGSE measurement of the example.
  • FIG. 78 is a diagram showing temporal changes in echo signal intensity measured by the anode coil in the CPMG measurement of the example.
  • FIG. 79 is a diagram showing temporal changes in echo signal intensity measured by the cathode side coil in the CPMG measurement of the example.
  • FIG. 80 is a diagram showing the relationship between the water content in P EM and the signal intensity in the example.
  • the current measurement mode is also referred to as a first measurement mode.
  • Figure 1 is a flowchart showing an overview of the current measurement procedure.
  • the following steps are performed sequentially, and the current at a specific part of the sample is measured locally using the nuclear magnetic resonance (N M R) method.
  • N M R nuclear magnetic resonance
  • the motion of nuclear magnetization can be detected as an NMR signal by the spin resonance phenomenon of a nucleus placed in a magnetic field. If the NMR signal is measured using a small surface coil (small RF coil), local NMR measurement around the coil becomes possible.
  • Step 3 0 1 Place the sample in the space where the magnet is placed and apply a static magnetic field to the sample.
  • Step 3 0 3 Apply a vibrating magnetic field for excitation to a specific part of a sample placed in a static magnetic field using a small RF coil that is smaller than the sample, and generate a nuclear magnetic resonance (NMR) signal generated at the specific part.
  • NMR nuclear magnetic resonance
  • Step 3 0 5 Calculate the difference between the frequency of the nuclear magnetic resonance signal acquired in Step 3 0 3 and the frequency of the oscillating magnetic field for excitation.
  • Step 3 07 Obtain the current at a specific part of the sample from the difference obtained in Step 3 0 5, and
  • Step 3 0 9 Then, output the result.
  • Step 3 0 3 Application of high frequency pulse for excitation and acquisition of NMR signal
  • a high frequency pulse applied to the measurement target nucleus in the sample is applied as an oscillating magnetic field for excitation.
  • the NMR signal emitted from the measurement target nucleus in the sample is acquired by the nuclear magnetic resonance phenomenon caused by the oscillating magnetic field for excitation.
  • the NMR signal is specifically an echo signal corresponding to the excitation high-frequency pulse. It is preferable that the phase of the echo signal is converged so that the frequency difference in Step 3 0 5 can be obtained reliably. Further, it is preferable to apply the high-frequency pulse in a pulse sequence in which the phases of the echo signals are matched.
  • the NMR signal is detected by separating the real part and the imaginary part by a phase sensitive detection method. As a result, the calculation of the frequency difference in step 300 is easily performed.
  • the difference (frequency shift) between the frequency of the NMR signal acquired in step 303 and the frequency of the oscillating magnetic field for excitation is obtained.
  • the phase difference ⁇ 0 is obtained by calculating the real part and the imaginary part of the echo signal acquired by the phase sensitive detection method. Then, the frequency shift ⁇ is converted as a phase difference ⁇ per unit time.
  • the current is calculated from the frequency difference ⁇ obtained in Step 3 0 5.
  • the calculation principle of the current will be explained.
  • the resonance frequency ⁇ of nuclear magnetization is directly proportional to the magnetic field strength H. If the magnetic resonance signal is acquired by the small detection coil, the small detection coil In other words, the magnetic field strength H in the region measured by the magnetic field is indirectly measured as the magnetic resonance frequency ⁇ .
  • the frequency difference ⁇ obtained in step 305 is obtained.
  • the current j flowing through the sample can be obtained.
  • the NMR detection method has a frequency resolution of the order of p pm, and this makes it possible to capture changes in magnetic field strength with high resolution and high sensitivity.
  • the frequency of the oscillating magnetic field for excitation is 43 MHz, a resolution of about 10 Hz can be obtained sufficiently.
  • the spin echo method is used, and the high frequency pulse for excitation is, for example, a pulse sequence including a plurality of pulses including the following (a) and (b).
  • the resonance excited magnetization vector M- y relaxes with time.
  • the time change of the magnetic resonance signal actually observed at this time is relaxed by the spin-lattice relaxation time constant T and another time constant ⁇ 2 * that cannot be expressed by the spin-spin relaxation time constant T 2 alone.
  • T spin-lattice relaxation time constant
  • ⁇ 2 * time constant that cannot be expressed by the spin-spin relaxation time constant T 2 alone.
  • spin echo is a method for correcting the phase shift due to magnetic field inhomogeneity as a sample or device characteristic. This is because the phase of the magnetization vector M is disturbed on the xy plane by applying a 180 ° excitation pulse with twice the excitation pulse intensity after the 90 ° excitation pulse. middle reverses the disturbance of the phase, after 2 Te time is a method of obtaining an echo signal to get on and converges the phase T 2 decay curves on.
  • the 180 ° excitation pulse applied in (b) above may be either the X direction or the Y direction, either the 180 ° excitation pulse. Can be used.
  • the phase of the magnetization vector can be converged and an echo signal as strong as possible can be obtained.
  • an echo signal the real and imaginary parts of the NMR signal can be detected with higher accuracy, and the amount of phase change from the reference frequency can be determined reliably.
  • the moisture content measurement mode is hereinafter also referred to as a second measurement mode.
  • FIG. 3 is a flowchart showing an outline of moisture content measurement.
  • the sample is placed in the space where the magnet is arranged, and a static magnetic field is applied to the sample (S 10 02).
  • an excitation oscillating magnetic field (high frequency pulse) is applied to the sample via a small RF coil, and the corresponding N MR signal (echo signal) is acquired (S 104).
  • a T 2 relaxation time constant is calculated from this echo signal (S 10 06). And from the obtained T 2 relaxation time constant, to measure the local water content in the sample (S 1 0 8). Specifically, to get the de-one data showing the correlation between the water content and the T 2 relaxation time constant of the sample, from this data and the T 2 relaxation time constant, contact to a particular location in the sample Find the local water content. After that, the result is output (S 1 1 0). By performing the above procedure (Step 104 to Step 110) via each small RF coil, the distribution of moisture content can be grasped.
  • step 104 to step 108 will be specifically described.
  • Step 104 Applying excitation RF pulse and acquiring N MR signal
  • the excitation high-frequency pulse in step 104 is preferably a pulse sequence composed of a plurality of pulses, and an echo signal group corresponding to the pulse sequence is acquired. In this way, the T 2 relaxation time constant can be obtained accurately.
  • the pulse sequence preferably includes the following (a), (b) and (c).
  • n 1 80 ° pulses (n is a natural number) that starts after the lapse of time 2 of the pulse in (b) and is applied at intervals of time 2
  • Hydrogen nuclei placed in a static magnetic field have a net magnetization vector in the direction along the static magnetic field (for convenience, the Z direction), and RF of a specific frequency (this is called the resonance frequency).
  • the resonance frequency By irradiating a wave from the outside in the X-axis direction perpendicular to the Z-axis, the magnetization vector tilts in the positive direction of the Y-axis, and a nuclear magnetic resonance signal (referred to as the NMR signal) can be observed.
  • the excitation pulse in the X-axis direction irradiated to obtain the maximum intensity N MR signal is called a 90 ° pulse.
  • the magnetization vector After tilting in the positive direction of the Y axis with a 90 ° pulse, irradiate an external 1 80 ° excitation pulse in the ⁇ axis direction '' after a while, and reverse the magnetization vector ⁇ with the axis symmetric '' Let As a result, after two hours, the magnetization vector converges on the “positive direction” of the ⁇ axis, and an N MR signal with a large amplitude is observed.
  • Figures 2 (a) to 2 (d) illustrate the compensation function of the spin echo method.
  • the coordinates shown in the figure are a rotating coordinate system.
  • the nuclear magnetization P which had a more advanced phase, has a phase that is later than Q, so that both nuclear magnetizations reach the y 'axis at the same time. (Fig. 2 (d)).
  • the magnetization vector converges on the “positive direction” of the Y axis, and an echo signal having a large amplitude is observed. Furthermore, in the above (c), after 3 hours, the magnetization vector is irradiated with an external 180 ° excitation pulse in the “Y-axis direction” and converged again on the “positive direction” of the vertical axis. Then, an echo signal with a large amplitude is observed after 4 hours. In addition, continue to irradiate 1 80 ° pulse at the same interval of two.
  • Step 1 06 (Calculation of T 2 relaxation time constant)
  • p is the density distribution of the target nuclide as a function of the position (X, y, z)
  • TR is the 90 ° excitation pulse repetition time (from about 10 Oms to about 10 s)
  • TE is the echo time (2ms, about 1 ms to 10 Oms)
  • A is a constant that represents the RF coil detection sensitivity and device characteristics such as amplifier.
  • Step 1 06 as described above, Step 1 04 plurality of echo signal groups riding on the obtained T 2 decay curves on the (on 2, 4 Te, Te 6, ⁇ ⁇ ⁇ ) an exponential By fitting with ⁇ 2, the relaxation time constant ⁇ 2 can be obtained from the above equation ( ⁇ ).
  • Step 1 08 (Calculation of water content)
  • the water content is calculated from the T 2 relaxation time constant.
  • the amount of water in the sample and the T 2 relaxation time constant have a positive correlation, and the ⁇ 2 relaxation time constant increases as the amount of water increases. Since this correlation varies depending on the type and form of the sample, it is desirable to prepare a calibration curve for a sample of the same type as the sample to be measured whose moisture concentration is known in advance. In other words, it is desirable to measure the relationship between the moisture content and the ⁇ 2 relaxation time constant for a plurality of standard samples with known moisture content, and obtain a calibration curve representing this relationship in advance. By referring to the calibration curve created in this way, the amount of water in the sample can be calculated from the measured value of the ⁇ 2 relaxation time constant.
  • the mobility measurement mode is also called the third measurement mode below.
  • a pair of gradient magnetic field pulses pulsed gradient magnetic field
  • N nuclei move due to Brownian motion or diffusion, and the phase of the nuclear spin does not converge, so the intensity of the N MR signal decreases.
  • the self-diffusion coefficient of a specific molecular species can be measured by associating the gradient magnetic field pulse changed stepwise with the decrease in the intensity of the NMR signal. This is the principle of measurement of the self-diffusion coefficient by the PGS E method.
  • FIG. 5 is a diagram showing an example of a PGS E sequence used for measuring the self-diffusion coefficient.
  • a pair of gradient magnetic field pulses Gz with the same application time and intensity are applied to the spin echo sequence described above with reference to FIG. N
  • a spin echo signal is acquired as the MR signal.
  • the peak intensity S of the obtained N MR signal is the applied pulse gradient magnetic field intensity G z
  • D, ⁇ , and G z indicate the pulse width of the gradient magnetic field pulse, the time interval between the pair of gradient magnetic field pulses, and the magnetic field gradient (z direction) of the gradient magnetic field pulse, respectively.
  • ⁇ represents the gyromagnetic ratio, which is a value intrinsic to the nucleus. For example, in the case of hydrogen nucleus 1 H, the gyromagnetic ratio 42. 577 X 1 0 2 [1 / ga
  • FIG. 6 is a flowchart for measuring the mobility of a specific portion of a sample using the PGS E method as described above, and includes the following steps.
  • an oscillating magnetic field for excitation and a gradient magnetic field are applied to the same region in the sample, and a corresponding N MR signal is acquired via a small RF coil (S 204).
  • Fig. 6 shows a case where no gradient magnetic field is applied in step 202, but in step 202, a gradient magnetic field having a magnitude different from that in step 204 is applied in accordance with a predetermined pulse sequence. May be executed. At this time, for example, it is preferable to set the magnitude of the gradient magnetic field in step 202 to a value close to zero.
  • the self-diffusion coefficient D is calculated from a plurality of N MR signals obtained by changing the gradient of the pulse gradient magnetic field stepwise (S206).
  • step 206 based on the self-diffusion coefficient D calculated in step 206, other parameters indicating the mobility of water in the sample may be calculated. After that, the result is output (S208).
  • Step 202 to Step 208 By performing such operations (Step 202 to Step 208) via each small RF coil, the distribution of the self-diffusion coefficient can be grasped.
  • Step 202 and Step 204 Application of excitation oscillating magnetic field, application of gradient magnetic field and acquisition of N M R signal
  • an excitation oscillating magnetic field and a gradient magnetic field are applied to the sample according to a predetermined pulse sequence.
  • the excitation oscillating magnetic field is a pulse sequence composed of a plurality of pulses
  • the gradient magnetic field is a pair of pulse sequences corresponding to the excitation oscillating magnetic field.
  • step 202 the gradient magnetic field is set to a value close to zero or close to zero, and in step 204, a predetermined gradient magnetic field is applied.
  • the pulse sequence is preferably composed of the following (a) to (d).
  • (b) the time to finish applying the gradient magnetic field pulse and (d) the time to start applying the gradient magnetic field pulse are as follows.
  • (C) 1 80 ° pulse (even though it has a width of 120 microseconds. Considering 60 microseconds at the center as the axis of symmetry), the time is equal ((34.5 m s- (1) 5 ms) / ⁇ 6.5 ms), and
  • the N MR signal corresponding to the pulse sequence is measured.
  • the peak intensity S of the NMR signal is measured by the spin echo method. Specifically, as shown in Fig. 5, the peak intensity S of the echo signal that appears in time 2 is measured.
  • the peak intensity S may be an average value of NMR signal intensities measured not only in the 2 hour NMR signal intensity but also in the surrounding time. This method can reduce variations in measured values caused by noise included in the NMR signal.
  • step 204 the step of executing the excitation oscillating magnetic field and the gradient magnetic field according to a predetermined pulse sequence, and the step of acquiring the NMR signal corresponding to this pulse sequence are performed once or a plurality of times. Run once.
  • step 206 the self-diffusion coefficient D of water at a specific location of the sample is obtained from the peak intensity of the NMR signal obtained in steps 202 and 204.
  • the self-diffusion coefficient D of the proton is expressed by the above formula (I I) using the peak intensity S of the NMR signal obtained by the PGS E method.
  • the sample The self-diffusion coefficient D of the inner proton can be obtained.
  • the same location in the sample is measured by changing the magnitude of the gradient magnetic field G, and I n ( By plotting the relationship between S / S 0 ) and 1 2 DA 2 d G 2 , the self-diffusion coefficient D can be obtained from the slope of the plot.
  • (B) measurement of moisture content and (C) measurement of mobility described above may be performed by switching measurement modes. It is also possible to calculate the distribution of water molecule movement based on the water content calculated in each measurement mode and the water mobility.
  • FIG. 7 is a diagram showing a schematic configuration of the measuring apparatus 300 according to the present embodiment.
  • Each component of the measuring apparatus 300 is realized by an arbitrary combination of hardware and software, mainly a CPU, a memory, and a program that implements the components shown in FIG. It will be understood by those skilled in the art that there are various variations in the implementation method and apparatus.
  • Measuring device 300 is a device that locally measures the current at a specific location of sample 1 15 using the NMR method.
  • a static magnetic field application unit (magnet 1 1 3) that applies a static magnetic field to sample 1 1 5, an oscillating magnetic field for excitation to sample 1 1 5, and N generated at a specific location on sample 1 1 5 Small RF coil 1 1 4 smaller than sample 1 1 5 to acquire MR signal, and
  • the current calculation unit 303 calculates the difference ⁇ between the frequency of the N MR signal acquired by the small RF coil 1 1 4 and the frequency of the oscillating magnetic field for excitation, and calculates the current at a specific location of the sample 1 1 5 from the difference. Is provided.
  • the measurement apparatus 300 includes an RF oscillator 10 02, a modulator 104, an RF amplifier 10 06, a preamplifier 1 1 2, a detector 301, an A / D converter 1 1 8, and a switch unit 1.
  • the measuring device 300 may have a structure described later with reference to FIG.
  • Sample 1 1 5 is a sample to be measured.
  • the sample 1 15 can be in various forms such as a membrane, a solid such as a lump, a liquid, agar, a gel such as a jelly.
  • the measurement result of local water content can be obtained stably.
  • a membrane having a property of retaining moisture in the membrane such as a solid electrolyte membrane, is used as a sample, the measurement result can be obtained more stably.
  • the sample mounting table 1 1 6 is a table on which the sample 1 1 5 is mounted, and a sample having a predetermined shape and material can be used.
  • Magnet 1 1 3 applies a static magnetic field to sample 1 1 5 (S 3 0 1 in FIG. 1). With this static magnetic field applied, an oscillating magnetic field for excitation is applied to the sample, and the current is measured.
  • the small RF coil 1 1 4 applies an excitation oscillating magnetic field to a specific part of the sample 1 1 5 and acquires an NMR signal corresponding to the excitation oscillating magnetic field (S 3 0 3 in Fig. 1). )
  • the NMR signal is a high-frequency pulse for the excitation oscillating magnetic field to generate nuclear magnetic resonance.
  • the small RF coil 1 1 4 is preferably 1/2 or less, more preferably 10 or less, of the size of the entire sample. By using such a size, it becomes possible to accurately measure the local mobility of the proton solvent in the sample in a short time.
  • the sample size can be, for example, the projected area when the sample is placed, and the exclusive area of the small RF coil 1 14 is preferably 1/2 or less of the projected area, more preferably By setting it to 1/10 or less, accurate measurement can be performed in a short time.
  • the size of the small RF coil 1 14 is preferably, for example, a diameter of 1 O mm or less. As the small RF coil 1 1 4, for example, the one shown in FIG. 33 (a) described later in the embodiment can be used.
  • the planar coil as shown in the figure, it is possible to limit the measurement region and perform local measurement.
  • the measurement area of such a spiral coil has a width of about the coil diameter and a depth of the coil radius, for example.
  • this coil is different from an ordinary solenoid coil, and is flat, so that NMR signals can be acquired simply by pasting on a flat sample.
  • the small RF coil 1 1 4 is not limited to a flat spiral coil, and various types of RF coils can be used.
  • a planar 8-shaped coil (sometimes called a butterfly coil, Doub I e _D type coil, etc.) is also available.
  • the figure 8 coil includes two spiral coils, and the NMR from the sample, even if the coil's spiral axis is parallel to the direction of the main magnetic field of the magnet, or when both are at an angle. The signal can be detected.
  • the spiral coil has sensitivity in the axial direction of the wound coil, whereas the figure 8 coil has sensitivity in the same plane as the wound coil.
  • the oscillating magnetic field (exciting oscillating magnetic field) applied by the small RF coil 1 1 4 is the RF oscillator 1 02, the modulator 1 04, the RF amplifier 1 06, the pulse controller 1 08, It is generated by the linkage of the switch unit 1 61 and the small RF coil 1 1 4.
  • the RF pulse generator that generates an RF pulse for generating an oscillating magnetic field for excitation in the small RF coil 1 14 includes an RF oscillator 10 2, a modulator 1 04, and an RF amplifier 1 06. Composed.
  • the excitation oscillating magnetic field oscillated from the RF oscillator 102 is modulated by the modulator 104 based on the control by the pulse control unit 108 and becomes a pulse shape.
  • the generated RF pulse is amplified by the RF amplifier 106 and then sent to the small RF coil 1 14.
  • the reference frequency is matched to the resonance frequency of the NMR signal when no current is flowing. This resonance frequency is stored in the RF oscillator 102.
  • the pulse control unit 10 8 controls the above cooperation so that the excitation oscillating magnetic field applied to the sample 1 1 5 by the small RF coil 1 1 4 is executed according to the pulse sequence described above. To do.
  • the pulse control unit 1 0 8 is connected to the sequence table 1 2 7 and the time measuring unit 1 2 8, and the sequence data acquired from the sequence table 1 2 7 and the measurement time in the time measuring unit 1 2 8 are used. Based on this, a high frequency pulse is generated.
  • the sequenceable 1 2 7 stores the sequence data of the high frequency pulse when measuring the current.
  • the sequence table 1 27 stores a timing diagram in which the generation time and interval of the high-frequency pulse are set, and the intensity of the high-frequency pulse to be applied based on the timing diagram.
  • the small R F coil 1 1 4 applies this R F pulse to a specific location of the sample 1 1 5 placed on the sample mounting table 1 1 6.
  • the small RF coil 1 1 4 acquires the NM R signal of the applied RF pulse.
  • the NMR signal is, for example, an echo signal corresponding to an oscillating magnetic field for excitation.
  • the frequency of the echo signal changes from the reference frequency described above due to the magnetic field formed by the flow of current. For this reason, the current flowing through the sample 115 is obtained from the difference in the frequency of the measured echo signal by acquiring the relationship between the frequency change (difference) and the current value in advance.
  • the frequency difference can be obtained by converting the amount of phase change at a certain time interval per unit time.
  • the exciting oscillating magnetic field applied to the sample 1 1 5 by the small RF coil 1 1 4 is, for example,
  • the small RF coil 1 14 When the small RF coil 1 14 is used, it may be difficult to adjust the excitation pulse intensity of the above (a) and (b). For example, in the region to be measured, that is, the region surrounded by the small RF coil 1 1 4 It is difficult to excite so that the whole is at a uniform excitation angle, that is, the intensity ratio of the excitation magnetic field in (a) and (b) is constant. There is a case. If the excitation angle ratio in (a) and (b) varies, an appropriate spin echo signal cannot be obtained, and accurate measurement of the current becomes difficult.
  • (a) a 90 ° pulse is in the first phase
  • (b) 1 80 ° pulse is in a second phase that is 90 ° away from the first phase.
  • the N MR signal detection unit detects the N MR signal acquired by the small RF coil 1 14 and sends this NMR signal to the calculation unit 130.
  • the NMR signal detector includes a preamplifier 1 1 2, a detector 301 and an A / D converter 1 1 8. The detected NMR signal is amplified by the preamplifier 1 1 2 and then sent to the detector 301.
  • the detector 301 is configured to detect the real part and the imaginary part of the NMR signal by the phase sensitive detection method.
  • the phase difference between the sine wave and the cos wave of the fundamental wave that is the source of demodulation is precisely 90 degrees. It is preferable to adjust to. By adjusting the two fundamental waves so that the phase difference is exactly 90 degrees, the phase difference can be calculated more accurately using the real part and imaginary part tan- 1 described later. Note that the reference wave that is the source of demodulation is Made by the 90 ° hybrid described.
  • the detector 30 1 sends the detected real part and imaginary part to the A / D converter 1 1 8.
  • the A / D converter 1 1 8 converts the N MR signal from A / D and sends it to the data reception unit 1 3 1.
  • the calculation unit 130 including the current calculation unit 303 acquires the data sent to the data reception unit 1 3 1.
  • FIG. 8 is a diagram showing an example of such an LC circuit.
  • the coil part (inductance part) of the resonant circuit is a small RF coil with a diameter of 1.4 mm.
  • NMR nuclear magnetic resonance
  • the number density and spin relaxation time constant can be measured by detecting the movement of nuclear magnetization as an N MR signal by the spin resonance phenomenon of the nucleus placed in a magnetic field.
  • the spin resonance frequency in the magnetic field of 1 Tes Ia is about 43 MHz, and an LC resonance circuit as shown in Fig. 8 is used to selectively detect the frequency band with high sensitivity.
  • the switch section 1 6 1 is provided at the branch section connecting the small RF coil 1 1 4, the RF amplifier 1 06 and the preamplifier 1 1 2, and the small RF coil 1 1 4 and RF The first state where the signal generator (RF amplifier 106) is connected, and
  • the switch section 1 6 1 serves as a “transmission / reception switching switch”. The role of this is to transmit the excitation pulse amplified by RF p ower-amp to the small RF coil 1 1 4, disconnecting the receiving preamplifier 1 1 2 and protecting it from a large voltage, and N after excitation
  • the switch section 1 6 1 is required for the following reasons in order to handle weak signals.
  • the cross diode is a diode that is turned on when a voltage equal to or higher than a predetermined value is applied, and is turned off when the voltage is lower than the predetermined value.
  • the sample volume that can be detected by the small coil of this measurement system is smaller than that of the large coil. This detectable sample volume is approximately (coil inner area x coil radius depth).
  • This detectable sample volume is approximately (coil inner area x coil radius depth).
  • preamplifier 1 1 2 When using the high sensitivity preamplifier 1 1 2, preamplifier 1 1 2 must be disconnected so that preamplifier 1 1 2 can be protected from the high voltage excitation pulses sent to the small coil during transmission.
  • the excitation pulse intensity is very large and the loss at the cross diode can be considered to be negligibly small, but in the case of a small coil, the excitation pulse intensity is smaller than that of the large coil.
  • the loss at the cross diode cannot be ignored. For this reason, in order to obtain an appropriate excitation pulse intensity, a “transmission / reception switching switch” with minimal loss is required.
  • the RF detection coil for local measurement is miniaturized, and low noise during NMR reception is an important factor for ensuring the accuracy of measurement.
  • the noise that enters the preamplifier 1 1 2 mainly includes an RF wave transmission system. RF leakage from the RF amplifier 1 06 that amplifies the excitation pulse and high power There is noise generated by the amplifier.
  • the switch unit 1 61 can employ various configurations.
  • FIG. 9 is a circuit diagram showing an example of the configuration of the switch unit 1601.
  • the real part and the imaginary part of the N MR signal (echo signal) detected by the detector 301 are acquired by the data reception unit 1 31 and sent to the calculation unit 1 30.
  • the calculation unit 130 has a current calculation unit 303.
  • the current calculation unit 303 acquires the real part and the imaginary part of the echo signal detected by the detector 301, calculates the phase difference between the echo signal and the excitation oscillating magnetic field using these, and from this phase difference, the echo is calculated.
  • the difference between the frequency of the signal and the frequency of the oscillating magnetic field for excitation (frequency shift amount) ⁇ is calculated (S 305 in FIG. 1).
  • tan- 1 (Re / I mg) is calculated from the detected real part and imaginary part. This value corresponds to the phase difference ⁇ [rad] of the NMR signal.
  • ⁇ is the phase difference between the measured NMR signal and the reference wave (phase ⁇ 0 ) that travels at a frequency that does not change with time, as shown in Fig. 10.
  • the reference frequency is The resonance frequency of the NMR signal in a state where no flow is flowing is set in advance.
  • the current calculation unit 303 obtains ⁇ from the amount of change per unit time of the obtained phase difference ⁇ . Then, by referring to the relationship between ⁇ and the current, the value of the current of the sample 115 at the measurement location is calculated (S 307 in FIG. 1). In addition, the current calculation unit 303 may calculate the current density by dividing the obtained current value by the area where the current flows.
  • the measuring apparatus 300 includes a storage unit 305 that holds information indicating the correlation between the current and the frequency difference for each type of the sample 1 15.
  • the storage unit 305 stores, for example, data on correlation between the frequency difference ⁇ and current obtained experimentally. More specifically, this is calibration curve data of the frequency difference ⁇ and the current.
  • the current calculation unit 303 in the calculation unit 130 acquires calibration curve data corresponding to the sample to be measured from the storage unit 305, and calculates a current corresponding to the frequency difference ⁇ based on the calibration curve data.
  • the current calculated by the current calculation unit 303 is presented to the user by the output unit 135 (S309 in FIG. 1).
  • the presentation format can take various forms, and there are no particular restrictions on display, printer output, file output, etc.
  • FIG. 41 shows the RF oscillator 10 02, modulator 10 04, RF amplifier 106, pulse control unit 10 08, switch unit 161, small RF coil 1 1 4, in the measuring apparatus 300 shown in Fig. 7.
  • FIG. 3 is a diagram showing an example of a more detailed configuration regarding the cooperation of a preamplifier 1 1 2, a detector 301 and an A / D converter 1 1 8. This configuration can also be applied to the measuring apparatus shown in FIGS. 13 and 16 described later.
  • the modulator 104 includes a mixer 1 77, a mixer 1 79, and a synthesizer 1 81.
  • the detector 301 includes a mixer 1 83, a mixer — 1 85, and a distributor 1 87.
  • the A / D converter 1 1 8 includes a first _A / D converter 1 189 and a second A / D converter 191.
  • a 90 ° hybrid 1 7 1 and a distributor 1 73 are further arranged in this order between the RF oscillator 1 02 and the modulator 1 04, and the 9 0 ° hybrid 1 7 1 Distributor 1 75 is placed between the detector and detector 301 Has been.
  • the waveform output from RF oscillator 102 is changed to two waveforms of the same frequency but different in phase by 90 ° by 90 ° hybrid 1 71. Based on these two reference waveforms, the NMR signal is detected and becomes R e a I and I m i g n a r y components.
  • the two waveforms output from the 90 ° hybrid 1 7 1 are specifically a sine wave and a cos wave, and the two waveforms are accurately orthogonal to each other. It is an important point in seeking.
  • the signal name of A / D converter 1 1 8 is named Rea I and I ma ginary, but this is a convenient expression, I ma ginary And it may be the opposite of Real.
  • the phase obtained by a r c t a n is only shifted by ⁇ 90 °, and there is no problem in obtaining the “phase change” that increases or decreases with time.
  • the "N MR signal strength" required when determining the water content and mobility described later in the fourth and fifth embodiments is the component of the acquired Rea I and I ma ginary.
  • the static magnetic field uniformity in the measurement area is high, and the echo signal can be observed over a very long time.
  • the frequency shift amount can be measured with high frequency resolution.
  • the phase of the echo signal can be converged and acquired by the small RF coil 114.
  • the real part and the imaginary part of the echo signal can be detected and the phase difference ⁇ using them can be calculated more accurately.
  • the spin-echo method need not be used for the measurement of the frequency shift amount in other embodiments of the present specification, and the frequency shift amount can also be calculated from a simple FID (Free Induction Decay).
  • the spin echo method has a smaller measurement area than the FID, so that the uniformity of the static magnetic field can be further improved compared to the FID.
  • a static magnetic field may be created using a rod-shaped magnet, and the sensor unit is small and can be used as a current measurement probe that can be easily installed in the apparatus.
  • water content and “mobility of water molecules” in a sample such as a polymer membrane can be measured locally at the same place almost at the same time. it can.
  • This embodiment and the following embodiments can be applied to, for example, local current measurement of a solid electrolyte membrane of a fuel cell.
  • the current j measured by the method of the above embodiment is generated by the fuel cell.
  • the principle of forming a magnetic field is the same regardless of whether the current flows in a state where the electric field is flowing or when the direct current is applied to perform water electrolysis. Therefore, by measuring the increase / decrease ⁇ of the resonance frequency, the spatial current during power generation and water electrolysis operation of the fuel cell can be grasped.
  • the power generation state changes depending on the gas supply state, catalyst deterioration, and ionic conductivity of the polymer electrolyte membrane.
  • the hydrogen utilization rate is increased, the hydrogen concentration is high near the gas supply port and the generated current at that location is large.
  • the hydrogen concentration is low near the gas outlet and the generated current is also small. This is because “material transport loss” increases.
  • the “activation loss” increases and the generated current decreases.
  • the Pt catalyst deteriorates due to transient fluctuations such as starting and stopping of the fuel cell, resulting in spatial non-uniformity.
  • the ionic conductivity increases and decreases, the “ohm loss” changes, and the power generation current increases and decreases. Since this loss depends on the spatial distribution of water content, the current also has a spatial distribution even within a single polymer electrolyte membrane.
  • the state of power generation at each location can be determined, and what is the state of “material transport loss”, “activation loss”, and “ohm loss” spatially? Specifically, when the gas supply concentration, hydrogen utilization rate, gas supply pressure, humidification amount, and moisture content of the membrane are changed, the power generation state at which location changes. It is possible to examine in detail whether the resulting battery performance has increased or decreased. This can provide technical guidelines for improving battery performance.
  • the spatial current of the fuel cell can be grasped by measuring the increase / decrease ⁇ of the resonance frequency.
  • the CPMG method in addition to the current measurement, can be used to obtain the 2 relaxation time constant, whereby the amount of water in the sample can be calculated.
  • the measuring device 300 can be further provided with a gradient magnetic field coil, and by appropriately applying a gradient magnetic field to the sample, the mobility of the target molecule can be measured by using the PGS method with the gradient magnetic field applied. it can.
  • the measuring apparatus 300 described in the first embodiment may include a plurality of small RF coils 1 1 4.
  • the plurality of small RF coils 1 1 and 4 apply an excitation oscillating magnetic field to a plurality of locations of the sample 1 1 5 and acquire a nuclear magnetic resonance signal.
  • the current calculation unit 303 is configured to calculate currents at a plurality of locations of the sample 1 15.
  • FIG. 11 is a perspective view showing an arrangement example of a plurality of small RF coils 1 14.
  • the current distribution in the sample 1 1 5 can be measured.
  • the two-dimensional arrangement is performed along the surface of the sample 1 15, the two-dimensional current distribution on the sample surface can be obtained.
  • the three-dimensional current distribution in the sample can be obtained.
  • a current distribution calculation unit (not shown) that calculates the current distribution in the sample 1 15 based on the current calculation result in the calculation unit 1 3 0 force current calculation unit 3 0 3 is provided. It may be. Thereby, it is possible to apply the excitation oscillating magnetic field to a plurality of locations of the sample and acquire the corresponding NMR signals.
  • the current distribution calculation unit (not shown) calculates the current distribution in the sample based on the current at multiple points in the sample.
  • the output unit 1 3 5 outputs this distribution.
  • the output unit 1 3 5 may have the configuration shown in FIG.
  • the output unit 1 3 5 is a measurement data acquisition unit that acquires the current for each measurement region of the plurality of small RF coils 1 1 4 (FIG. 1 1) calculated by the current calculation unit 3 0 3. 1 3 5 A and a display section 1 3 5 B that displays the acquired current in a partitioned area of the same screen.
  • the screen is divided into a plurality of regions according to the arrangement position of small RF coil 1 1 4. Each area is displayed in a predetermined color according to the current in the measurement area of each small RF coil 1 14.
  • the plurality of areas of the display unit 1 3 5 B are divided into upper and lower parts, current is shown on one side (for example, the upper half), and in the sample 1 1 5 in the corresponding place on the lower half The amount of moisture may be indicated.
  • the measurement of the water content in the sample 1 15 will be described later in the fourth embodiment.
  • the measuring apparatus 300 described in the first embodiment may be an apparatus that acquires a current distribution in a plane of a solid polymer electrolyte membrane of a fuel cell by using a nuclear magnetic resonance method. Good.
  • Sample 1 15 is a solid polymer electrolyte membrane of a fuel cell.
  • the measuring device 300 is a measuring device for a fuel cell
  • the operating state of the fuel cell can be diagnosed. For example, when ⁇ ⁇ is measured in each region and compared with the theoretical analysis value, if the measured value of ⁇ ⁇ shows a different behavior from the analysis value only at a specific measurement location, There may be a defect in ME IV at the location.
  • ⁇ ⁇ is measured at a predetermined time interval for a plurality of measurement points, and the difference between the measured value of ⁇ and the analysis value becomes large at all the measurement points, ⁇ ⁇ ⁇ There is a possibility that the output is reduced overall.
  • the measurement accuracy of the frequency shift amount ⁇ can be further increased by the following procedure.
  • the reference frequency when the reference frequency does not change with time (constant) (for example, when a superconducting magnet or electromagnet is used as the magnet), the reference frequency is the first to be measured. It is only necessary to acquire it, and after that, it is only necessary to perform measurement when current is passed. In this case, there is an advantage that current fluctuation over time and transient response are easy to measure.
  • the magnetic field strength may fluctuate with time due to fluctuations in the temperature of the magnet.
  • the reference frequency cannot be said to be constant over time.
  • the fuel cell when the above embodiment is used for measurement of a fuel cell, the fuel cell does not immediately enter a steady state even if it starts power generation and starts flowing current, and it takes several seconds to several minutes to reach a steady state. May be required. Also, it may be difficult to apply to current measurement in “slowly occurring phenomena” where the load varies, the gas supply fluctuates, or moisture condenses in the gas diffusion layer. The reason is that the above-mentioned “NMR measurement with and without current at a short time interval where the reference frequency can be regarded as constant” cannot be performed.
  • the magnetic field distribution in the permanent magnet only rises and falls uniformly in the entire space. If the fluctuation of the nuclear magnetic resonance frequency is measured at a certain point, The nuclear magnetic resonance frequency of a place can be estimated by raising and lowering the fluctuation amount.
  • a plurality of small RF coils 1 14 are arranged as follows.
  • the case of measurement of power generation of a fuel cell is shown as an example, but the following method can be applied to the types of Sample 1 15 without any particular limitation.
  • Measurement is performed using a system that can receive NMR signals using multiple coils.
  • the fluctuation amount ⁇ (t1) at time t 1 is
  • ⁇ ⁇ (t1, x) ⁇ ourrent (t1, X) _OJ no — current (t1, X)
  • the current distribution should be analyzed based on ⁇ (t1, x) obtained as described above.
  • the following operational effects can be obtained. That is, by measuring the NMR signals for coil 1 and coil 2 at the same time, the reference frequency can be more accurately determined. Can be estimated. Thereby, the measurement accuracy of the frequency shift amount ⁇ ⁇ (t1, x) can be improved. For this reason, for example, when used for measurement of a fuel cell, it is not necessary to measure the fuel cell by switching between “no current in a short time”. This makes it possible to measure in a more realistic power generation situation.
  • Coil 1 and coil 2 need not be measured simultaneously. For example, measurement may be performed by alternately switching between coil 1 and coil 2.
  • the reference frequency w n for “no current” depending on the position of the coil. It is sufficient if _current (t1, x) can be estimated with a certain degree of accuracy.
  • current and water content are measured using a local magnetic resonance signal acquired by a small detection coil.
  • FIG. 13 is a diagram showing a schematic configuration of the measurement apparatus of the present embodiment.
  • the basic configuration of the apparatus shown in FIG. 13 is the same as that of the measuring apparatus 300 shown in FIG. 7 except that a solvent information calculation unit 3 09 is provided in the calculation unit 1 3 0. . Further, the difference is that a control unit 3 0 7 is provided instead of the pulse control unit 1 0 8 in FIG.
  • the solvent information calculation unit 3009 calculates information related to the solvent contained in the sample 1 15 and includes a water content calculation unit 1 32 in this embodiment. Based on the NMR signal acquired by the small RF coil 1 1 4 Calculate the amount of protonic solvent (water).
  • FIG. 14 is a diagram showing the configuration of the control unit 307 of the apparatus shown in FIG.
  • the control unit 307 includes a first measurement mode for measuring the current of the sample 1 15 and a second measurement for measuring the water content in the sample 1 15. Includes a switching unit (mode switching control unit 169) for switching between modes.
  • the operation signal receiving unit 129 connected to the mode switching control unit 169 receives the operator's request for the measurement mode. Then, the operation signal receiving unit 1 29 sends this request to the mode switching control unit 169.
  • the current measurement of the sample 1 15 is performed according to the procedure described in the above embodiment. That is, the current calculation unit 303 forces the calculation of the current at a specific location of the sample 1 15 based on the difference between the frequency of the N MR signal acquired by the small RF coil 1 1 4 and the frequency of the excitation oscillating magnetic field.
  • a small RF coil 1 1 4 Forced NMR signal (echo signal) corresponding to the vibration magnetic field for excitation is acquired, and the solvent amount calculation unit 1 32 (moisture amount calculation unit 1 32) However, it calculates the amount of the protonic solvent (water) in the sample 1 1 5 based on the echo signal acquired by the small RF coil 1 1 4. Specifically, the water content calculation unit 1 32 calculates the T 2 relaxation time constant from the intensity of the echo signal, and from the calculated ⁇ 2 relaxation time constant, the protonic solvent at a specific location in the sample 1 15 Calculate the amount of
  • a common pulse sequence is used for the first and second measurement modes.
  • the small RF coil 1 14 applies an oscillating magnetic field for excitation in a pulse sequence including the following (a) to (c).
  • n 1 80 ° pulses (n is a natural number) that starts after the lapse of time 2 of the pulse in (b) and is applied at intervals of time 2
  • the current calculation unit 303 obtains the real part and the imaginary part of the force echo signal and calculates the current.
  • small RF coil 1 1 4 force Acquires a plurality of echo signals corresponding to the pulses (b) and (c) above.
  • the water content calculation unit 1 32 calculates the T 2 relaxation time constant from the intensity of the plurality of echo signals.
  • FIG. 15 is a flowchart showing an example of a procedure for measuring current and water content. This measurement method includes the following steps.
  • Step 301 (Step 10 02): Applying a static magnetic field to sample 1 15;
  • Step 303 (Step 10 04): Pulse sequence including (a) to (c) above via small RF coil 1 1 4 Apply the excitation oscillating magnetic field and obtain the corresponding echo signal.
  • Step 305 Using the real part and imaginary part of the echo signal corresponding to the pulse (b) or (c) obtained in step 303, the frequency of the echo signal and the frequency of the oscillating magnetic field for excitation The difference between
  • Step 307 From the difference obtained in step 305, obtain the current at a specific part of the sample.
  • Step 1 06 Calculate the T 2 relaxation time constant from the intensities of multiple echo signals corresponding to the pulses (b) and (c) obtained in Step 303.
  • Step 1 08 ⁇ calculated in Step 1 06 2 Measure the local moisture content in the sample from the relaxation time constant.
  • Step 309 (Step 1 1 0): After that, the result is output.
  • Step 303 an echo signal is obtained for calculating the difference between the frequency of the nuclear magnetic resonance signal and the frequency of the oscillating magnetic field for excitation, and ⁇ 2 for calculating the relaxation time constant.
  • the acquisition of the echo signal may be performed simultaneously.
  • to calculate the difference calculation and T 2 relaxation time constant of the frequency both above (b)
  • step 3 0 3 step 1 0 4
  • an echo signal corresponding to the above pulses (b) and (c) is acquired.
  • the frequency difference is calculated by using the real part and the imaginary part of the echo signal corresponding to the pulse (b).
  • step 3 07 the T 2 relaxation time constant is calculated from the intensities of multiple echo signals corresponding to the pulses (b) and (c) above. Note that the pulse in (c) above can be used for both frequency difference calculation and ( 2) relaxation time constant calculation.
  • the present embodiment it is possible to measure not only the local current of the sample 1 15 such as a membrane, but also the water content with one apparatus by a series of measurements using a common pulse sequence. . For this reason, the state of the sample 1 15 during power generation or water electrolysis operation can be grasped in more detail.
  • the small RF coil 1 14 is a pulse that alternately repeats a pulse for current measurement (first measurement mode) and a pulse for moisture measurement (second measurement mode) multiple times.
  • An oscillating magnetic field for excitation can also be applied in a sequence. In other words, it performs an acquisition of echo signals for calculation unloading of the difference between the frequency and the frequency of the excitation oscillating magnetic field of the nuclear magnetic resonance signal, the acquisition of echo signals for calculating the T 2 relaxation time constant, alternating You can also In this way, the local current and water content of sample 1 15 can be measured more stably.
  • current, water content, and water mobility are measured using local magnetic resonance signals acquired by a small detection coil.
  • FIG. 16 is a diagram showing a schematic configuration of the measurement apparatus of the present embodiment.
  • the basic configuration of the apparatus shown in FIG. 16 is the same as that of the measuring apparatus shown in FIG. 13.
  • the solvent information calculation unit 3 0 9 of the calculation unit 1 30 is further replaced with the mobility calculation unit 1 3 3.
  • the movement amount calculation unit 1 3 4 is provided in the calculation unit 1 3 0.
  • the apparatus shown in FIG. 16 has a gradient magnetic field application unit (a pair of G coils 1 5 1) that applies a gradient magnetic field to sample 1 15. And a pair of G coils 1 5 1 Further prepare.
  • the pair of G coils 15 1 is a gradient magnetic field application coil disposed away from the small RF coil 1 1 4. As shown in FIG. 17, the pair of G coils 15 1 is arranged so that a gradient magnetic field can be applied to the sample 1 15. Two G coils 15 1 are arranged with respect to one small RF coil 1 1 4, and are arranged opposite to each other with the small RF coil 1 1 4 interposed therebetween.
  • the G coil 15 1 has a half-moon shape as shown in FIG. 17 illustrates a case where a plurality of small RF coils 1 1 4 are provided in one sample 1 1 5, and a pair of G coils 1 5 1 is arranged for each small RF coil 1 1 4. Yes.
  • the G coil 15 1 is arranged in parallel to the surface of the sample 1 15.
  • the G coil 15 1 is arranged above the small RF coil 1 14. As a result, a gradient magnetic field having a magnetic field gradient in the y-axis direction can be formed on the central axis of the small RF coil 1 14.
  • a shield shield (not shown) is provided between 14 and the other G coil 15 1. This shielding shield prevents noise from the G coil 15 1 from affecting the small RF coil 1 14.
  • the shielding shield has a thickness that prevents the passage of noise and allows a magnetic field to pass.
  • the mobility calculation unit 1 3 3 is based on the NMR signals acquired by the small RF coils 1 1 4 obtained corresponding to different gradient magnetic fields. Calculate the mobility of the protonic solvent (water).
  • the movement amount calculation unit 1 3 4 is based on the water amount calculated by the moisture amount calculation unit 1 3 2 and the self-diffusion coefficient calculated by the mobility calculation unit 1 3 3. Calculate the travel distance. For example, the movement amount calculation unit 1 3 4 calculates the movement amount of water molecules. A parameter storage unit in which parameters for storing are stored, and a movement amount calculation unit that reads a calculation formula stored in the parameter storage unit and calculates a movement amount of water molecules.
  • the parameter storage unit stores a calculation formula for calculating the amount of movement of water molecules from the self-diffusion coefficient and the amount of water for each sample 1 15 type. Based on this calculation formula, the movement amount calculation unit can calculate the movement amount.
  • the mode switching control unit 169 in the control unit 307 is a first measurement mode in which the current of the sample 1 15 is measured, and the water content in the sample 1 15 is measured. Switch between the second measurement mode and the third measurement mode that measures the mobility of water in sample 1 1 5.
  • an excitation oscillating magnetic field is applied to the small RF coil 1 1 4 force sample 1 15 and a nuclear magnetic resonance signal corresponding to the excitation oscillating magnetic field and gradient magnetic field is acquired.
  • the mobility calculation unit 133 calculates the mobility of a specific portion of the sample 1 15 based on the information of the nuclear magnetic resonance signals obtained corresponding to different gradient magnetic fields.
  • the small RF coil 11 14 applies the excitation oscillating magnetic field in a pulse sequence including the following (a) to (d).
  • the gradient magnetic field applied in (b) may be zero.
  • (a) 90 ° pulse is in the first phase and (c) 1 80 ° pulse is in the second phase that is 90 ° shifted from the first phase.
  • FIG. 18 is a flowchart showing an example of a procedure for measuring current and mobility.
  • FIG. 18 is an example of using the gradient magnetic field applied in (b) of the above pulse sequence as an outlet, and includes the following steps.
  • Step 3 0 1 (Step 1 0 2): Applying a static magnetic field to the sample 1 1 5, Step 3 0 3 (Step 2 0 2): Pulse with zero gradient field and the above (a) to (d) In sequence, apply the excitation oscillating magnetic field and acquire the corresponding echo signal (first step)
  • Step 30 Difference between the frequency of the echo signal and the frequency of the oscillating magnetic field for excitation using the real and imaginary parts of the echo signal corresponding to the pulse (d) obtained in step 303 To calculate,
  • Step 3 07 From the difference obtained in Step 3 0 5, find the current at a specific part of the sample (second step)
  • Step 2 04 The gradient magnetic field is set to a non-zero predetermined magnitude, and the excitation oscillating magnetic field is applied in the pulse sequence including the above (a) to (d), and the corresponding echo signal is acquired (third step) ),
  • Step 2 06 From the peak intensities of the NMR signals obtained in Step 2 0 2 and Step 2 0 4, the self-diffusion coefficient D of water at a specific location in Sample 1 1 5 Seeking (fourth step)
  • Step 3 0 9 (Step 2 0 8): After that, the result is output.
  • the gradient magnetic field in step 2 0 2 is set to zero.
  • the current and self-diffusion coefficient D can be obtained by a series of measurements.
  • the current and moisture content of sample 115 can be obtained by a single measurement by the procedure described above with reference to Fig. 15 for example.
  • the water content and the mobility of water molecules can also be measured. For this reason, the state in the membrane at the time of power generation or water electrolysis operation can be grasped in more detail.
  • the "current distribution" in the fuel cell, the "water content” in the polymer electrolyte membrane, and the "mobility of water molecules” can be Or it can measure alternately. Measurement that integrates these is useful as an integrated monitoring method for fuel cells, and is effective as a device that provides new measurement amounts to development sites aimed at improving fuel cell performance and expands the application range of NMR sensors. It is.
  • the amount of current flowing in a sample such as a fuel cell can be converted from the frequency shift amount of the NMR signal received by the small detection coil. At that time, even with one small detection coil, if the current flowing through the sample is uniform, the amount of current can be easily converted from the amount of frequency shift.
  • multiple small detection coils can be incorporated into a single unit cell, and NMR signals can be acquired from the polymer film, and the current distribution can be converted from the frequency shift amount ⁇ of the NMR signal that occurs when current flows. It becomes. This makes it possible to effectively measure the current distribution with multiple coils.
  • the “local water content”, “local mobility of water molecules” and “current distribution” of the polymer membrane for fuel cells are measured locally using the same device and the same sensor. By alternately measuring the three quantities every few seconds, both values can be obtained at almost the same time (same device, same position, almost same time).
  • non-invasive measurement using electromagnetic waves is performed simply by attaching to the surface of the polymer film.
  • the amount of water, the mobility of water molecules, and the current distribution can be measured while generating power from the fuel cell.
  • the state of the polymer membrane can be grasped from various information such as “water content”, “mobility of water molecules”, and “current distribution”, and the power generation state or water electrolysis state in the fuel cell Monitoring can be performed to control power generation efficiency to the highest level.
  • the “Magnet ⁇ Magnetic field gradient coil integrated rod-shaped local measurement sensor” allows the sensor to be easily installed in the fuel cell, and only the RF detection coil section needs to be in the measurement area. Measurements can be made without obstruction (Fig. 19).
  • the frequency shift amount is It can also be calculated from a simple FID (Free Induct on Decay).
  • F ID Free Induct on Decay
  • the current calculation unit 303 may obtain a real part and an imaginary part of the F I D signal to obtain a current from a phase change amount. In this way, measurement can be performed with only 90 ° pulses, so the repetition time of the excitation pulse can be shortened, and the amount of phase change can be determined in a short time. In addition, since the F ID signal itself can be observed whether it is smaller or larger than 90 °, adjustment can be simplified compared to the spin echo method.
  • the current when measuring the amount of water or the mobility of water together with the current, the current can be calculated from the FID signal corresponding to (a) by using the pulse sequence including (a) above. .
  • Example 1 In this example, water was poured into a narrow gap, a copper plate was placed in close contact with the sample, and the change in frequency of the NMR signal when current was passed through the copper plate was measured.
  • FIG. 20 is a diagram for explaining the outline of this embodiment.
  • an electric current j was passed through a copper plate to form a magnetic field! "!
  • an experiment was conducted to measure the frequency change of the N MR signal from the water placed next to it. It was.
  • the copper plate simulates a fuel cell
  • the water simulates a polymer membrane in the fuel cell.
  • the sample is a sealed container with two cover glass with dimensions of 18 mm x 18 mm and a thickness of 0.1 2 mm, with a gap of 0.5 mm, and 2.5 mm mo I in it. / encapsulating C u S0 4 aqueous solution L.
  • Figure 21 shows the sample used in this example.
  • the interval between the 90-degree excitation pulse and the 180-degree excitation pulse was set to 5 ms, and the echo signal was measured with an echo time of 1 Oms.
  • a gradient magnetic field was applied for 1 ms before and after the 180 ° excitation pulse so that the NMR signals immediately after the 90 ° and 180 ° excitation pulses did not interfere with the echo signal.
  • Figure 22 shows the echo signal when the current j is zero.
  • the NMR signal was detected by the phase-sensitive detection method, and two signals were acquired, the real part and the imaginary part.
  • the real part and the imaginary part are indicated by "R e a I” and "I mag", respectively.
  • phase difference ⁇ [rad] of the NMR signal was calculated using ta n- 1 (R e / I mg) based on the real part R e and the imaginary part I mg.
  • the phase reference is the reference wave from the oscillator of the NMR device, and this frequency is set in advance to the resonance frequency of the NMR signal. It was.
  • the phase difference between the reference wave (phase ⁇ 0 ) that does not change in time and the measured NMR signal is ⁇ .
  • the relationship between the real part, the imaginary part, and the phase difference ⁇ is shown in FIG.
  • FIG. 23 shows the phase difference ⁇ calculated from the tan- 1 (Re / I mg) in FIG. However, in this figure, 1 of the time when the echo signal is observed
  • the echo signal measured when the current j is 0.8 OA is shown in FIG. 24, and the phase difference ⁇ calculated based on this is shown in FIG.
  • the pulse sequence is the same as when the current j is closed.
  • the real part and imaginary part of the NMR signal oscillate, and the frequency is shifted from the reference wave.
  • the real part comes first and then the imaginary part vibrates.
  • FIG. 25 shows the phase difference ⁇ calculated by obtaining tan- 1 (Re / I mg) in FIG.
  • the phase difference ⁇ increases (upward straight line), and you can see how the phase of the N MR signal advances from the reference wave over time.
  • Fig. 26 shows the echo signal measured at 10.8 OA when current j is applied in the reverse direction
  • Fig. 27 shows the phase difference ⁇ calculated based on this.
  • Figure 27 is calculated by obtaining ta n- 1 (Re / I mg) in Figure 26.
  • the phase difference ⁇ is shown.
  • the phase difference ⁇ decreases with the passage of time (straight-down straight line), and the phase of the NMR signal delays with the passage of time from the reference wave. it can.
  • Fig. 28 shows the results of experiments conducted in increments of 0.20A with the current j flowing through the copper plate varied from 0.80A to + 0.80A.
  • the phase difference ⁇ of the NMR signal that changes in 1 ms is defined as “frequency shift amount of N MR signal ⁇ [rad / ms]”.
  • This “frequency shift amount ⁇ [rad / ms]” corresponds to “slope of phase difference ⁇ ” in FIGS. 23, 25, and 27, and the graph of phase difference ⁇ is linearly approximated by the method of least squares. It was calculated from the gradient.
  • N The factor that increases or decreases the MR signal frequency over time is not only the current j flowing through the copper plate, but also the increase or decrease in magnetic field strength in the case of permanent magnets. If the temperature of the permanent magnet increases or decreases, the magnetic field strength also increases or decreases in inverse proportion to it. For this reason, a time-stable magnetic field is required for frequency shift measurement.
  • the permanent magnet has a large heat capacity, and the frequency change due to a sudden temperature change
  • the phase difference can be measured with higher accuracy.
  • this method can also be used by subtracting the frequency of the reference wave even if it is set to a state slightly deviated from the true resonance frequency of the N MR signal. There is a feature that can offset the “deviation”.
  • 1 ⁇ 1 ⁇ 8 is a £ 1 ⁇ 1 (Polymer Electrolyte Membrane) with electrodes joined.
  • the ME A used here was manufactured by attaching Pt and Ir on the anode side and Pt on the cathode side electrolessly on a polymer electrolyte membrane manufactured by Asahi Glass.
  • the dimensions of M E A are 17 mm x 15 mm square, 500; U m thick.
  • ME A was standardized, pulled up from the ion-exchanged water just before the experiment, and wiped off the water moderately. From the T 2 (CPMG) relaxation time constant of ME A just before water electrolysis operation, the water content of MEA was about 10 [H 2 0 / S 0 3 -H +].
  • ME A was sandwiched between a poly force single-point cell with a small surface coil and a Pt electrode, and power was applied.
  • the voltage applied between the two electrodes was 2 to 3.5 V, and the current j was 0.10 to 0.3 OA.
  • ME A decomposes the water contained in the PEM and releases hydrogen and oxygen.
  • water was not supplied during water electrolysis (no humidification condition).
  • the cell temperature was 24 ° C.
  • Fig. 29 shows the echo signal measured when 0.30 A of current j is passed through the MEA
  • Fig. 30 shows the phase difference ⁇ calculated based on this echo signal. From FIG. 29 and FIG. 30, it can be seen that the echo signal acquired from ME A progresses although the phase difference is slight with time.
  • Figure 32 shows an overview of the entire system.
  • the water to be measured is sandwiched between two glass covers (15 mm x 15 mm) and sealed to a thickness of 0.5 mm.
  • two copper plates are placed in close contact, and their dimensions are 19 mm x 9 mm.
  • Each copper plate is connected to a constant current power source, and can individually control the current amount, 1 2 (j 2 in the figure).
  • FIGS. 33 (a) to 33 (c) The coils and devices used in this example are shown in FIGS. 33 (a) to 33 (c).
  • FIG. 33 (a) is a diagram showing a small surface coil used in this example.
  • FIG. 33 (b) is a diagram showing a pair of polycarbonate holders used in this example.
  • FIG. 33 (c) is a diagram showing an RF coil holder used in this example.
  • FIG. 34 is a diagram showing the static magnetic field H Q applied to the sample by the permanent magnet and the direction of the magnetic field created by passing a current through the copper plate. If current I 2 flows in the positive direction (in the direction of the arrow in the figure) through the copper plate, static magnetic field H. Is applied to the sample, and the frequency shift ⁇ increases.
  • N MR signal was obtained by changing current and I 2 independently, and the dependence of frequency shift amount ⁇ on currents I and I 2 was measured.
  • the frequency shift amount of coil 2 is about 3.92 rad / ms larger than that in Fig. 35, but it is inversely proportional to current I, its slope is negative, and its slope is You can see that they are almost the same.
  • the shape of the copper plate is symmetric between 1 and 2, and the positions of coils 1 and 2 are the same distance from the central axis of the poly force cell. Therefore, the measurement results (Fig. 3 8) are shown in Fig. 3 Matches the result of swapping 1 and 2.
  • the above (i) is a method using the fact that the frequency shift amount strongly depends on the current amount near the coil from the above measurement result.
  • the total current value for example, I and + I 2 are required, but this can be easily measured. According to this method, a more accurate calculation of the current distribution can be expected.
  • Example 3 the current value was measured when two small coils were placed using the spin echo method. Specifically, the current was obtained from the amount of phase change of the echo signal (equivalent to the frequency shift ⁇ ) at a certain time interval.
  • the NMR signal is not limited to the echo signal, and the current can be measured by the same method even from the F ID signal.
  • Fig. 43 shows the phase of the FID obtained by calculating arctan (Real / Imaginary).
  • time 1 FID up to 5ms It can be seen that the phase can be calculated without much dispersion.
  • the FID in Fig. 45 has a waveform that decays with a shorter T 2 * relaxation time constant.
  • the difference in calculating the amount of phase change observed with FID and echo signals is as follows. That is, in the FID, the T 2 * relaxation time constant increases or decreases, and the observation time of the N MR signal whose phase can be calculated is increased or decreased by the amount of current. On the other hand, in the case of an echo signal, it is almost constant. It is. From the viewpoint of analyzing the obtained N MR signal, it is preferable that “the time during which the N MR signal can be analyzed so that the phase can be analyzed is constant” like the echo signal.
  • Figure 48 shows the FID waveform obtained when the current was increased to 0.8 OA.
  • the FID in Fig. 48 has a waveform that decays with a shorter T 2 * relaxation time constant compared to Figs. 42 and 45.
  • phase change at a certain time interval calculated from the F I D waveform, that is, ⁇ was 8.01 r a d / ms.
  • the observed FID shows that the current flows in the opposite direction and the direction of the magnetic field is reversed, so that the progression of the Real and I magnary waveforms is reversed compared to Fig. 45. I understand.
  • Fig. 52 shows the phase of the FID obtained by calculating arctan (Real / Imaginary). In this figure, the phase It can be seen that the retreat (progress with a negative slope).
  • the frequency shift amount in coil 1 in Fig. 35 is compared with that in Fig. 54, the straight lines of both are almost the same slope, and the value of the intersection with the vertical axis is also the same. I understand. As a result, the frequency shift amount with respect to the current is almost the same regardless of whether it is an FID or an echo signal, and the relational expression (calibration equation) between the current and the frequency shift amount is created by either method. If this is done, it can be used to calculate the current value from the frequency shift amount regardless of the method.
  • FIG. 55 is a diagram showing directions of a static magnetic field H 0 and a magnetic field H i described later.
  • the conductor is placed in the static magnetic field and flows from left to right in Fig. 55.
  • This current creates a magnetic field around the conductor.
  • the magnetic field applied to the sample is the sum of the static magnetic field Ho by the magnet and Hi formed by the current.
  • the frequency ⁇ [H z] of the N MR signal is proportional to the magnetic field strength H [ga u s s], as shown in the following equation.
  • ⁇ is the nuclear gyromagnetic ratio [H z / gauss], and in the case of 1 H hydrogen nucleus, it is 4260 H z / gauss.
  • the magnetic field Hi that increases or decreases when the current I flows through the conductor is obtained as the frequency shift amount ⁇ [H z].
  • FIG. 56 is a perspective view showing the positions of the copper plate, RF detection coil, and water sample in the magnetic field analysis of this example.
  • the magnetic field Hi created by the current I flowing through the conductor can be calculated based on Bio-Savart's law. Specifically, when the conductor is placed in a vacuum (magnetic permeability is 4 ⁇ X 1 0- 7 N / A 2), conductors position (x p, y p, z p
  • FIG. 57 is a diagram showing a coordinate system in the following equation.
  • the copper plate is divided into small elements, and it is assumed that the current divided equally flows through the elements. Specifically, there is a small element at point Q, a current flows for each element, and these multiple elements at point Q create a magnetic field at point P.
  • the magnetic field at point P was calculated by integrating all elements at point Q using the above equation.
  • the longitudinal direction of the copper plate (L direction in Fig. 57) was equally divided into 64, and the width direction (W direction in Fig. 57) was equally divided into 32.
  • FIG. 58 shows the analysis results.
  • Figure 58 shows the results when the current I is changed to _ 1 [A], -0.5 [A], 0 [A], 0.5 [A], and 1 [A].
  • the unit of the frequency shift amount ⁇ on the vertical axis is [rad / ms].
  • H z was multiplied by 27 ⁇ to convert to rad, and s was divided by 1 000 to convert to ms.
  • FIG. 59 is a diagram showing the relationship between the frequency shift amount ⁇ and the current I at the position of the water sample measured by the RF detection coil. From Fig. 59, it can be seen that the relationship between the two is directly proportional. It can also be seen that the current I can be easily converted from the frequency shift ⁇ .
  • Example 5-1 the amount of frequency shift was measured using a small surface coil and compared with the result of Example 5-1.
  • a copper plate, an RF detection coil, and a water sample (pure water) were placed in the positional relationship shown in Fig. 55, and the relationship between the current I flowing through the copper plate and the frequency shift amount ⁇ was measured. .
  • FIG. 60 is a diagram showing a small surface coil produced in this example.
  • N MR measurement system a base system manufactured by MRL Technology was used with high sensitivity.
  • the magnet used was a modified Halbach magnetic circuit with a magnetic field strength of 1. O Tes Ia and an air gap of 45 mm manufactured by Neomax.
  • FIG. 61 is a diagram showing the relationship between the current I [A] flowing through the copper plate and the measured frequency shift amount ⁇ [rad / ms].
  • Fig. 61 shows the magnetic field of Fig. 59. The relationship (solid line) obtained from the analysis results is also shown.
  • Example 6_2 prior to actual measurement in Example 6_2 described later, the magnetic field was analyzed when ME A was modeled with a one-dimensional equivalent circuit.
  • FIG. 62 is a cross-sectional view showing a schematic configuration of MEA used in the analysis in this example.
  • the MEA is an assembly of an electrode and a polymer electrolyte membrane, and the polymer electrolyte membrane (PEM) is sandwiched between upper and lower current collectors.
  • PEM polymer electrolyte membrane
  • the dimensions of ME A are 23 mm X 2 Omm square, thickness 356
  • the current collector was a one-pong mesh with a thickness of 30 Om.
  • Pt and Ir were electrolessly attached to one surface of the polymer electrolyte membrane to form an anode-side catalyst layer.
  • Pt was unelectrolyzed as a catalyst on the other surface of the polymer electrolyte membrane to form a force sword side catalyst layer.
  • the RF detection coil is embedded in the lower side (electrolyte membrane side) of the bonbon mesh and is in contact with PEM with catalyst inside ME A. Also The RF detection coil is sandwiched between the upper force mesh and the lower force mesh. The magnetic field created by these two conductors is shown in Fig. 5.
  • the current distribution flowing through the anode-side collector electrode and the cathode-side collector electrode is a distribution that decreases or increases linearly with respect to the X direction. This was assumed to be an equivalent circuit of ME A.
  • FIG. 64 is a diagram showing the frequency shift amount ⁇ analyzed with respect to the positions of the cross section a_a and the cross section b_b shown in FIG.
  • the horizontal axis represents the frequency shift (rad / ms)
  • the vertical axis represents the position z (mm) in the thickness direction of the PEM.
  • the frequency shift amount at the position where the RF detection coil ⁇ is placed is about 6 rad / ms from the upper side to the lower side. Decreases from about 1 rads.
  • the current I 1 (at s e c t i o n a _ a) force flowing in the upper collector electrode is larger than the current I 2 (at s e c t i o n a _ a) flowing in the lower collector electrode (I ⁇ >
  • Example 6-1 the frequency shift distribution when the ME A shown in FIG. 62 was water electrolyzed was measured and compared with the analysis results of Example 6-1.
  • Pt and Ir are electrolessly attached to one side (anode side) of a polymer electrolyte membrane manufactured by Asahi Glass Co., and Pt is electrolessly attached to the other side (force sword side). Produced.
  • the dimensions of ME A are 23 mm x 2 Omm square and thickness 356 m.
  • the obtained MEA was sandwiched by a force of one meter mesh (manufactured by Japan Gore-Tex) with a thickness of 300 m.
  • the 0.03 mm thick Pt electrode foil was also sandwiched.
  • the current from the power source was changed from the conductor to the Pt electrode foil, Force —Bonmesh, Pt electrode foil, and lead wire flow in this order to return to the power source.
  • the applied voltage during water electrolysis was about 3 V, and the current density was 0.26 A / cm 2 .
  • the temperature of ME A during water electrolysis operation was about room temperature. In this embodiment, no steam is supplied to ME A.
  • a N MR measuring surface coil having an inner diameter of 0.6 mm was used. A small hole was made in the bonbon mesh, and the lead part of the coil was passed there to fix the coil on the bonnet mesh.
  • FIG. 65 is a diagram showing a coil used in this example.
  • the small surface coil was made by winding a copper wire with a polyurethane film diameter of 40 m in a spiral shape with an inner diameter of 0.6 mm five times in a flat shape. This coil was manufactured by Star Engineering.
  • FIG. 66 is a view showing a coil placed on a carbon mesh in this embodiment.
  • the amount of frequency shift in coil A is indicated by a white square (mouth), and the amount of frequency shift in coil B is indicated by a white triangle ( ⁇ ).
  • the measurement area in the depth direction of the coil is about one-fifth of the coil diameter, so the area has a width of about 0.1 mm from the PEM surface. This width is shown as a bar in Figure 67.
  • the frequency shift amount could be measured at multiple points even in the system corresponding to the operating fuel cell.
  • this method it is possible to know the local distribution of the current in the MEA of the fuel cell.
  • the current measurement method based on NMR was applied to ME A (Membrane Electrode Assembly) in which electrodes were joined to PEM, and the frequency shift of the NMR signal was measured by PEM during water electrolysis operation.
  • ME A Membrane Electrode Assembly
  • magnetic field analysis was also performed and compared with the measurement results.
  • the magnetic field analysis in ME A was performed, and the frequency shift amount of the N MR signal that increases or decreases due to the magnetic field was calculated.
  • the magnetic field analysis was performed with an electrode arrangement in which current was supplied to the end of the MEA and the current flowed from the opposite end of the force sword.
  • the protons pass through the PEM uniformly, have the same resistance value in the entire area of the PEM, and further, current flows through the electrodes on the anode side and cathode side of the PEM. It is assumed that a closed circuit is constructed.
  • the current I, flowing through the anode electrode decreases linearly with position X, and conversely, the current I 2 flowing through the force source electrode increases linearly with position X.
  • the current flowing through the anode and cathode electrodes forms a magnetic field inside the PEM.
  • This magnetic field strength distribution Hi was analyzed using Bio-Savart's law, and the frequency shift amount of the N MR signal was calculated from the magnetic field strength.
  • the origin of position X is the center of PEM.
  • Figure 69 shows the z-direction distribution of frequency shifts analyzed at the four sensor positions.
  • Position (x _7. 5mm) at the anode side of the current I, of the sensor _ A but larger than the current I 2 forces cathode side, as a result, the frequency shift amount has a positive value.
  • the force sword side The current I2 of the current is also large on the far node side, and the frequency shift amount in that case is a negative value opposite to that of the sensor A.
  • Sensors _B and C located in the middle of both are the frequency shift amount in the middle.
  • the ME A used in the water electrolysis operation was manufactured by electrolessly attaching Pt and Ir on the anode side of the polymer electrolyte membrane and Pt on the force sword side.
  • the dimensions of ME A are 23 mm x 20 mm square and 356 m thick.
  • this MEA is sandwiched as a GDL (Gas Diffusion Layer) by a 300m thick force pon mesh, and current is supplied from the stabilizing power source through the Pt electrode.
  • GDL Gas Diffusion Layer
  • a structure is used in which pressure is applied from both sides of the positive and negative electrodes using a cushioning material so that GDL and MEA are in uniform contact.
  • Applied voltage current density of about 3 V was 0. 26A / cm 2.
  • the surface coil for NMR measurement was a copper wire with an inner diameter of 0.6 mm and a wire diameter of 0.04 mm, and four coils were placed at 5 mm intervals between the MEA and the force pon mesh. Polyurethane coating is applied to the surface coil wire, and the carbon mesh is in an insulated state.
  • the measurement area in the coil depth direction is about one-fifth of the coil diameter, and is a disk-like area with a depth of about 0.1 mm from the PEM surface.
  • the representative value of the frequency of the NMR signal measured by this coil is considered to be the center of this disk, and its depth is assumed to be 0.05 mm.
  • Figure 71 shows the measured frequency shift amount and the X-direction distribution of the frequency shift amount obtained from the analysis.
  • the horizontal axis in this figure is position X, which corresponds to the position of each sensor.
  • the vertical axis represents the frequency shift amount.
  • the country (black square) force ⁇ , sensor _A ⁇ D of the position X of the frequency shift measured value is shown.
  • the solid line shows the amount of frequency shift calculated using the analysis method shown in the (MEA equivalent circuit) section and considering the thickness of GDL. In this case, the frequency shift
  • a small surface coil is inserted between the GDL and PEM in the fuel cell, and the current value when the fuel cell is operated in water electrolysis is measured from the frequency shift amount of N MR. The water content was also measured. If water electrolysis is performed without supplying water, the PEM will gradually dry and the current will gradually decrease. In this example, the spatial distribution of current and water content in this case was measured in time series.
  • FIG. 72 is a diagram showing the arrangement of MEAs and small coils in the present embodiment. As shown in Fig. 72, small coils are inserted between PEM and GDL, and the number is 3 on the anode side and 1 on the force side.
  • the direction of the static magnetic field HQ is opposite to that of Example 7.
  • the cell direction was reversed so that the direction of the static magnetic field was reversed.
  • the sign of the frequency shift amount measured in Example 7 is reversed.
  • the absolute value of the frequency shift amount for sensor A is the same, but the sign is negative.
  • MEA was immersed in distilled water until immediately before the experiment, and the surface water was wiped off with Kimwipe just before the experiment to obtain an appropriate water content.
  • the water content of ME A when treated in this way is about 1 0 [H 2 0 / S 0 3 ⁇ H +]. Even in this experiment, the water content is considered to be about 10 [H 2 0 / S0 3 -H +].
  • FIG. 73 is a diagram showing the measurement timings of PGS E and C PMG in this example when the time when water electrolysis operation was started (voltage was applied) was set to zero.
  • the echo signal obtained before applying a voltage to ME A is required. Therefore, in this example, as shown in FIG. 73, a series of measurements was started before applying a voltage (“O ff” in the figure) and used as a reference echo signal. At this time, the temperature of the magnet increases or decreases with time, the static magnetic field strength increases or decreases, and the frequency of the N MR signal may change. In this measurement, the amount of change in frequency was about 1 OOH z in 1 hour. The frequency that changes in the experiment time of 150 seconds, which was conducted this time, is about a few Hz, which is sufficiently smaller than the frequency shift amount (about 1 kHz) that increases or decreases with the current. For this reason, it is assumed that the reference frequency is the same as the value (echo signal waveform) obtained before applying the voltage. Can do. Here, we calculated the frequency shift amount based on this assumption.
  • an NMR detection coil is installed in a place that is not affected by the magnetic field formed by the current flowing through the MEA.
  • the time change of the static magnetic field strength can be offset.
  • FIG. 74 is a diagram showing the time change of the current flowing through M EA.
  • FIG. 75 is a diagram showing the time change of the voltage applied to MEA.
  • the current density at this time was 0.25 A / cm 2 .
  • the voltage applied to ME A reached 3.4 V, and at the same time, the current decreased to about 0.8 A. After that, the current gradually dropped to about 0.5 A.
  • PGSE measurement is performed six times at a time interval of about 30 seconds, which corresponds to the time indicated by the arrows in FIG. 74 (PGS E # 1 to # 6). In this way, the frequency shift was measured while the current was decreasing. On the other hand, C PMG was measured 5 times during P G S E.
  • FIG. 76 is a diagram showing the frequency shift amounts obtained by the three sensors A, C, and D from the measurements in PGS E # 1 and # 4 (FIG. 74).
  • the horizontal axis represents the sensor position x and the vertical axis represents the frequency shift amount, as in FIG. 71 shown in Example 7 above.
  • the plot for country shows the frequency shift measured with PGSE # 1 (current is 1.2 A).
  • the ⁇ (black triangle) plot shows the frequency shift measured with PGS E # 4 (current is 0.6 A).
  • the solid line and the alternate long and short dash line in the figure are analysis values calculated by the same magnetic field analysis as in Example 7. [0350] The result shown in Fig.
  • Fig. 71 differs from the result of Fig. 71 in the sign of the frequency shift amount, because the static magnetic field direction is reversed.
  • the direction of the current flowing through the GDL is the same, but if the direction of the static magnetic field is reversed, the sign of the frequency shift amount is reversed. If the measurer knows the direction of the static magnetic field, it will not be a problem.
  • Fig. 77 is a diagram showing the time change of the frequency shift amount obtained by the measurement from PGS E # 1 to # 6.
  • the measured and analyzed values of the three sensors A, C, and D are also shown. From this figure, the measured values and the analyzed values are in good agreement. Therefore, even when the current flowing through the MEA decreases transiently from 1.2A to 0.6A, it can be inferred that the current flows uniformly through the MEA.
  • the time change of the echo signal intensity acquired on the anode side is shown in Fig. 78, and the echo acquired on the force sword side (sensor 1E) is shown in Fig. 78.
  • Figure 79 shows the signal strength over time.
  • the echo signal intensity is the intensity of the echo signal that was observed second using the C PMG method, and the average signal intensity when five C PMG measurements were performed.
  • the values on the vertical axis are normalized by the average signal strength of 4 sets measured before applying current.
  • Fig. 8 1 is a diagram for explaining a phenomenon occurring in P E M during water electrolysis operation of M EA.
  • Fig. 81 water moves to the power lead side by electroosmotic flow in the PEM, and further water is decomposed by electrolysis on the anode side, so the water content on the anode side decreases. To do.
  • the force sword side the water content once increases due to electro-osmotic flow.
  • the water content of the entire P EM will decrease due to electrolysis, and eventually the water content will gradually decrease with time even on the power source side.
  • the frequency shift measurement by the PGSE method and the moisture content measurement by the CPMG method are shown alternately.
  • the frequency shift measurement and the moisture content measurement are common. You may carry out simultaneously by a pulse sequence.
  • the moisture content is measured by the frequency shift measuring device is shown, but the mobility of water can also be measured by the frequency shift measuring device.

Abstract

L'invention concerne un dispositif de mesure (300) qui mesure localement un courant à une position particulière d'un échantillon (115) à l'aide du procédé RMN. Le dispositif de mesure comprend : un aimant (113) qui applique un champ électrostatique à l'échantillon (115) ; une bobine RF de petite taille (114) plus petite que l'échantillon (115), laquelle bobine applique un champ magnétique oscillant pour l'excitation de l'échantillon (115) et acquiert un signal RMN généré à la position particulière de l'échantillon (115) ; et une unité de calcul de courant (303) qui calcule une différence entre la fréquence du signal RMN acquis par la bobine RF de petite taille (114) et la fréquence du champ magnétique oscillant d'excitation, et calcule le courant à la position particulière de l'échantillon (115) à partir de la différence.
PCT/JP2007/001055 2006-09-29 2007-09-28 Dispositif de mesure et procédé de mesure utilisant un procédé de résonance magnétique nucléaire WO2008041361A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008537406A JP5170686B2 (ja) 2006-09-29 2007-09-28 核磁気共鳴法を用いた測定装置および測定方法

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2006-268976 2006-09-29
JP2006268976 2006-09-29
JP2007126039 2007-05-10
JP2007-126039 2007-05-10
JP2007200654 2007-08-01
JP2007-200654 2007-08-01

Publications (1)

Publication Number Publication Date
WO2008041361A1 true WO2008041361A1 (fr) 2008-04-10

Family

ID=39268227

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/001055 WO2008041361A1 (fr) 2006-09-29 2007-09-28 Dispositif de mesure et procédé de mesure utilisant un procédé de résonance magnétique nucléaire

Country Status (2)

Country Link
JP (1) JP5170686B2 (fr)
WO (1) WO2008041361A1 (fr)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009302040A (ja) * 2008-05-13 2009-12-24 Keio Gijuku 燃料電池システム
JP2011002283A (ja) * 2009-06-17 2011-01-06 Keio Gijuku 測定装置および燃料電池システム
JP2013519073A (ja) * 2010-02-01 2013-05-23 クリア−カット メディカル エル・ティー・ディー 体外サンプルの断端評価
JP2013536940A (ja) * 2010-08-31 2013-09-26 メッツォ オートメーション オイ 固形物及びスラリーの含水量を測定する低磁場ｎｍｒ装置
JP2014098716A (ja) * 2014-02-19 2014-05-29 Keio Gijuku 測定装置および燃料電池システム
JP2019018148A (ja) * 2017-07-14 2019-02-07 三菱重工業株式会社 分離膜装置の汚損度評価方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60190846A (ja) * 1984-03-10 1985-09-28 Jeol Ltd 核磁気共鳴装置
JP2007051990A (ja) * 2005-08-19 2007-03-01 Keio Gijuku 核磁気共鳴法を用いて膜の透過特性を測定する測定装置および測定方法
JP2007121037A (ja) * 2005-10-26 2007-05-17 Keio Gijuku 磁気共鳴法を用いて試料中のプロトン性溶媒の挙動の分布を測定する測定装置、測定方法およびプログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60190846A (ja) * 1984-03-10 1985-09-28 Jeol Ltd 核磁気共鳴装置
JP2007051990A (ja) * 2005-08-19 2007-03-01 Keio Gijuku 核磁気共鳴法を用いて膜の透過特性を測定する測定装置および測定方法
JP2007121037A (ja) * 2005-10-26 2007-05-17 Keio Gijuku 磁気共鳴法を用いて試料中のプロトン性溶媒の挙動の分布を測定する測定装置、測定方法およびプログラム

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
BOGDAN BUHAI ET AL.: "Dissimilar Electro-Osmotic Flow and Ionic Current Recirculation Patterns in Porous Media Detected by NMR Mapping Experiments", PHYSICAL REVIEW LETTERS, vol. 96, no. 17, 2006, pages 174501-1 - 174501-4, XP003022111 *
HATADA H. ET AL.: "FEM ni Motozuita MRI Satsuzoji no Nonai RF Denjikai Bunpu", JOURNAL OF MAGNETICS SOCIETY OF JAPAN, vol. 29, no. 3, 1 March 2005 (2005-03-01), pages 364 - 367, XP003022110 *
OGAWA K. ET AL.: "Kogata Hyomen Coil ni yoru Nenryo Denchiyo Kotai Kobunshi Denkaishitsu Makunai no Suibunshi no Kyokusho Jiko Kakusan Keisu", NETSUKOGAKU CONFERENCE KOEN RONBUNSHU, no. 05-17, 2 November 2005 (2005-11-02), pages 237 - 238, XP003022109 *
OGAWA K. ET AL.: "Multi NMR Sensor ni yoru Kotai Kobunshi Denkaishitsu Makunai Gansuiryo Bunpu no Real Time Keisoku (Real-time measurement of water-content distribution in PEM under water electrolysis condition using multi-NMR sensor)", DAI 10 KAI NMR MICROIMAGING STUDY GROUP KOEN YOSHISHU, July 2006 (2006-07-01), pages 32 - 35, XP003022107 *
OGAWA K. ET AL.: "Multi NMR Sensor ni yoru Suidenkai Untenji no PEM-nai Gansuiryo no Bunpu Keisoku (Distribution measurement of water-content in PEM under water electrolysis condition using multi-NMR sensors)", NETSUKOGAKU CONFERENCE KOEN RONBUNSHU, no. 06-2, 23 November 2006 (2006-11-23), pages 245 - 246, XP003022108 *
SCOTT G.C. ET AL.: "Sensitivity of Magnetic Resonance Current-Density Imaging", JOURNAL OF MAGNETIC RESONANCE, vol. 97, 1992, pages 235 - 254, XP000275059 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009302040A (ja) * 2008-05-13 2009-12-24 Keio Gijuku 燃料電池システム
JP2011002283A (ja) * 2009-06-17 2011-01-06 Keio Gijuku 測定装置および燃料電池システム
JP2013519073A (ja) * 2010-02-01 2013-05-23 クリア−カット メディカル エル・ティー・ディー 体外サンプルの断端評価
JP2013536940A (ja) * 2010-08-31 2013-09-26 メッツォ オートメーション オイ 固形物及びスラリーの含水量を測定する低磁場ｎｍｒ装置
JP2014098716A (ja) * 2014-02-19 2014-05-29 Keio Gijuku 測定装置および燃料電池システム
JP2019018148A (ja) * 2017-07-14 2019-02-07 三菱重工業株式会社 分離膜装置の汚損度評価方法

Also Published As

Publication number Publication date
JPWO2008041361A1 (ja) 2010-02-04
JP5170686B2 (ja) 2013-03-27

Similar Documents

Publication Publication Date Title
JP4997380B2 (ja) 試料中のプロトン性溶媒の易動性を局所的に測定する方法、試料中のプロトン性溶媒の易動性を局所的に測定する装置
Siaw et al. A versatile and modular quasi optics-based 200 GHz dual dynamic nuclear polarization and electron paramagnetic resonance instrument
JP4798350B2 (ja) 磁気共鳴法を用いて試料中のプロトン性溶媒の挙動の分布を測定する測定装置、測定方法およびプログラム
JP5170686B2 (ja) 核磁気共鳴法を用いた測定装置および測定方法
Xu et al. Construction and applications of an atomic magnetic gradiometer based on nonlinear magneto-optical rotation
JP4849623B2 (ja) 試料中のプロトン性溶媒量を局所的に測定する方法、装置
US20130088232A1 (en) Solid-State NMR Spectrometer, Sample Holder Therefor, and Method of Solid-State NMR Spectroscopy
Liu et al. Pulsed-field nuclear magnetic resonance: Status and prospects
JP5337569B2 (ja) 燃料電池システム
JP4997620B2 (ja) 核磁気共鳴法を用いて膜の透過特性を測定する測定装置および測定方法
Tsushima et al. Magnetic resonance imaging of water in operating polymer electrolyte membrane fuel cells
Koppe et al. Frequency-swept ultra-wideline magic-angle spinning NMR spectroscopy
JP5212972B2 (ja) 計測装置および計測方法
Aguilera et al. The parallel-plate resonator: an RF probe for MR and MRI studies over a wide frequency range
JP5337413B2 (ja) 燃料電池用測定装置および燃料電池システム
JP6917251B2 (ja) 測定解析装置、燃料電池システム及び測定解析方法
JP2016145723A (ja) セルモデル、計測システム及び同時計測方法
JP5257994B2 (ja) 測定装置および測定方法
US20240094316A1 (en) Hyperpolarization micro-magnetic resonance imager and three-dimensional imaging of a biological composition with cellular resolution
CN217085237U (zh) 一种利用流动的液体来精确测量磁场的仪器
JP5513783B2 (ja) 測定装置および燃料電池システム
JP2007080613A (ja) 燃料電池内電流分布計測用スケール冶具及び燃料電池内電流分布計測用固定装置
JP5046203B2 (ja) 測定装置およびこれを備える燃料電池、ならびに測定方法
JP5046203B6 (ja) 測定装置およびこれを備える燃料電池、ならびに測定方法
JP2014098716A (ja) 測定装置および燃料電池システム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 07827834

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2008537406

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 07827834

Country of ref document: EP

Kind code of ref document: A1