A High Duty Cycle Transmitter Unit for Steady-State Surface NMR Instruments

Gaikwad, Nikhil B.; Liu, Lichao; Griffiths, Matthew P.; Grombacher, Denys; Larsen, Jakob Juul

doi:https://doi.org/10.5194/gi-2023-5

Preprints

https://doi.org/10.5194/gi-2023-5

Preprints

29 Jun 2023

| 29 Jun 2023

Status: a revised version of this preprint is currently under review for the journal GI.

A High Duty Cycle Transmitter Unit for Steady-State Surface NMR Instruments

Nikhil B. Gaikwad, Lichao Liu, Matthew P. Griffiths, Denys Grombacher, and Jakob Juul Larsen

Abstract. Groundwater measurements using surface nuclear magnetic resonance (NMR) have been notoriously challenged by a poor signal-to-noise ratio (SNR), but a new steady-state methodology based on long, high duty cycle, phase-locked pulse trains have demonstrated huge SNR increases. The hardware requirements for transmitters for steady-state surface NMR are significantly increased compared to transmitters for standard surface NMR use, due to the need for very high pulse-to-pulse stability over long survey times and the increased thermal load caused by a much higher duty cycle. Further, the increased SNR leads to increased production rates, which necessitates lightweight equipment, that can be carried easily between many field sites during surveys. Here we demonstrate a novel steady-state surface NMR transmitter with a maximum 93 A peak current. The stability of the transmitter is evaluated on 10 minutes, 10 % duty cycle, pulse trains containing pulses of either 5 ms, 10 ms, 20 ms, or 40 ms duration, and low or high current. We observe less than 150 ns pulse-to-pulse timing jitter and amplitude variations below 0.4 % between pulses for all pulse durations and currents. During tests, we observed no temperature effects on the timing and current stability. We have designed a customized heatsink, which reduces the transmitter weight by 30 % and size by 16 % without compromising safe thermal operating conditions. We evaluate the capacitor bank size and current stability and demonstrate that a 10 mF capacitor bank is an appropriate trade-off with insignificant current drooping in measurements. The extensive analysis and verification demonstrate that the transmitter generates highly stable pulse trains resulting in high SNR signals.

Received: 04 May 2023 – Discussion started: 29 Jun 2023

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.

Download & links

Nikhil B. Gaikwad, Lichao Liu, Matthew P. Griffiths, Denys Grombacher, and Jakob Juul Larsen

Status: final response (author comments only)

RC1:
'Comment on gi-2023-5', Anonymous Referee #1, 06 Sep 2024
Abstract
The authors present the development of a transmitter used for steady-state surface Nuclear Magnetic Resonance (NMR) applications in groundwater exploration. The new design and particularly the method using high-duty cycle enables to reduce the SNR as well as the time of data acquisition.

One of the main contribution is the design of a customized heatsink that allows the instrument to keep its performances while reducing its weight and size. Authors present a thermal simulation of the system in order to design their heatsink and a validation on field.

General comments
The authors are first thanked for their contribution. In particular, the article is well written and comprehensive with many details that significantly aid in the reader's understanding. However, some parts need clarification, notably in the results section where it is not always clear if the results are from the simulation part or the experiment one. Furthermore, graphics are not always easy to read and the data representation could be improve to further support the authors arguments, see the specific comments below.
Under what environmental conditions were the field experiments conducted? Specifically, do you expect your customized heatsink to perform as effectively as the standard one in extreme conditions, such as during a heatwave?

Specific comments
⚠️ notations on p.7: $T_j$ and $T_p$ are confusing. Use capital `T` for temperature (T_j) and `t` (t_p) for time (even if it is a time period)

⚠️ the way to reference the figures is not homogeneous. Sometimes it is entitled "Fig." while other times it is "Figure". Authors need to check the manuscript. Also, the references are not actual active links.

On Section 5, line 228, you state that "The temperature essentially follows the trend obtained in the simulation with poor and passive cooling". What about the other cooling methods? Could you provide a comparison between the simulation and the actual measurements for both the standard Heatsink and the customized one? I believe that it would be interesting to add some simulations results on Figure 6 to support your arguments.

In Table 2, it is not clear if the values are measured ones or simulated ones. Could you explain how you obtained all critical design parameters? In particular, the values on this table are difficult to link with the results on Figure 5.

In Table 2, it would be important, for reproducibility purposes, to specify what "Extreme Test Conditions" are.

Figure 8 is hard to read, and particularly Figure 8a. The negative time must be removed. As previously, adding the comparison between the measurement and simulation for 10mF (i.e. plotting the simulation with the measurements) would back up your statements.

Figure 10 is referenced before Figure 9

Figure 9b is only showing the first intensity (~70A) compared to Figure 9a (~70A & ~20A).
Citation: https://doi.org/10.5194/gi-2023-5-RC1
- AC1: 'Reply on RC1', Jakob Juul Larsen, 13 Nov 2024
  
  Dear Reviewer 1
  
  Thanks for the kind review and constructive comments for the manuscript.
  
  Main comment:
  
  Under what environmental conditions were the field experiments conducted? Specifically, do you expect your customized heatsink to perform as effectively as the standard one in extreme conditions, such as during a heatwave?
  
  The experimental data presented in the paper was acquired outside of the laboratory under calm conditions and an ambient temperature of about 15^oC. We will add a comment about this to the revised paper.
  However, this is not the full answer to this question. Much more importantly, since this manuscript was submitted, our instrument with the customized heatsink has been successfully used for extensive field work in, e.g., Senegal. Here, more than 100 soundings have been conducted at temperatures above 30^oC and more than 50 soundings have been conducted at temperatures above 40^oC. We can therefore now safely conclude that the performance of the customized heatsink is good. We will add a comment on our field experiences to the concluding remarks of the manuscript.
  Specific comments
  
  1. notations on p.7: $T_j$ and $T_p$ are confusing. Use capital `T` for temperature (T_j) and `t` (t_p) for time (even if it is a time period).
  
  We will revise the manuscript according to these instructions.
  
  2. The way to reference the figures is not homogeneous. Sometimes it is entitled "Fig." while other times it is "Figure". Authors need to check the manuscript. Also, the references are not actual active links.
  The journal instructions states that: The abbreviation "Fig." should be used when it appears in running text and should be followed by a number unless it comes at the beginning of a sentence, e.g.: "The results are depicted in Fig. 5. Figure 9 reveals that..." We have checked the manuscript, and we will correct one instance where we depart from the requested formatting. Regarding the actual links in the manuscript, we have used the LaTeX bibliography style from the editor. We will of course ensure that all references are properly formatted in the revised manuscript.
  
  3. On Section 5, line 228, you state that "The temperature essentially follows the trend obtained in the simulation with poor passive cooling". What about the other cooling methods? Could you provide a comparison between the simulation and the actual measurements for both the standard Heatsink and the customized one? I believe that it would be interesting to add some simulations results on Figure 6 to support your arguments.
  It is a good suggestion, but unfortunately, we did not record useable temperature data for the standard heatsink before it was replaced with the customized heatsink. We will expand the text and add more details on the similarities and dissimilarities between the simulations and the experiments.
  
  4. In Table 2, it is not clear if the values are measured ones or simulated ones. Could you explain how you obtained all critical design parameters? In particular, the values on this table are difficult to link with the results on Figure 5.
  These are simulated values. We will update the caption, so this is immediately apparent. Second, the wording used in the text, i.e., “critical design parameters” will be changed. As such, the only design parameter is the thickness of the heatsink. The other values are then either computed directly, e.g., the mass, or found from the simulations, e.g., the temperature at the junction.
  
  5. In Table 2, it would be important, for reproducibility purposes, to specify what "Extreme Test Conditions" are.
  
  The extreme test conditions correspond to continuous pulsing at 93 A for 3000 s. In reality, we never pulse continuously at the highest current for more than 60 s. The text will be revised to make this point clear.
  
  6. Figure 8 is hard to read, and particularly Figure 8a. The negative time must be removed. As previously, adding the comparison between the measurement and simulation for 10mF (i.e. plotting the simulation with the measurements) would back up your statements.
  The negative time will be removed and the text surrounding Figure 8a will be improved so it better guides the reader in understanding the figure. Please note that Figures 8a and 8b are plotted on the same time axis, so the time to stabilize the current can be directly compared between the simulations and the measurement. However, the simulations and measurements are not plotted on the same panel as the initial and final currents do not match perfectly. Importantly, we can’t expect them to match perfectly as, e.g., the DC supply for the capacitor is assumed an ideal source – which it is not.
  
  7. Figure 10 is referenced before Figure 9.
  This is a simply latex copy-paste error, and it will be fixed in the revision.
  
  8. Figure 9b is only showing the first intensity (~70A) compared to Figure 9a (~70A & ~20A).
  This is a deliberate choice made to avoid having an excessive number of panels on the figure. The variations in the amplitude of the ~20 A pulses are similar to the variations in amplitude of the ~70 A pulses, so we only show the high amplitude results. We will add a comment about similar variations in the low amplitude pulses to the revised manuscript.
  
  Citation: https://doi.org/10.5194/gi-2023-5-AC1
RC2:
'Comment on gi-2023-5', Anonymous Referee #2, 23 Oct 2024

Review of "A High Duty Cycle Transmitter Unit for Steady-State Surface NMR Instruments"
by

Gaikwad, Nikhil B., Liu, Lichao, Griffiths, Matthew P., Grombacher, Denys and Larsen, Jakob Juul
First of all, I would like to admit that I (a Geophysicist with experience in SNMR) have rather limited expertise in electrical engineering, and therefore cannot judge the hard details of the circuitry (i.e. Fig.3) or the implementation of the Simulink simulations. However, from the explanation the authors make and the results they present, I tend to say that I understood enough to give a proper review.
General comments:

The authors present their novel development of a customized heat sink that decreases the weight and size of their SNMR device and increases therefore its mobility in the field. They present the concept behind the measurement approach and show with simulations and field tests that the newly developed heat sink has no negative influence on the quality of the measured data.
The manuscript is generally well structured with a slight tendency to jump a bit back and forth in the text regarding figure descriptions. Also, I think the number of figures is a bit high and maybe there is a possibility to merge some figures into one or prepare them in a more compact manner. For me as a non-native English speaker the quality of the manuscript (spelling) is high with some minor typos here and there.
I think that all my comments can be addressed by a minor revision after which I recommend this manuscript for publication.
Specific comments:

Do you really need the TEM part in the introduction?
L.33: To me, relaxation time is only indirectly controlled by porosity. It is controlled by the pore size and mineral components of the pore walls. The amplitude is more closely connected to porosity because of: NMR amplitude -> (mobile) water content -> porosity
L.57: Isn't the 5s argument not a bit stretched? This strongly depends on the T1 of your subsurface and I can rarely remember ever seeing T1 above 500ms in the field. So maybe it is more fair to link the experiment time to the T1 time of the subsurface?
L.60-68: I totally agree on your point about the new steady-state approach and the massive SNR enhancement. Yes, this helps us to get better data to estimate water content (and therefore porosity) but still lacks the possibility to infer proper relaxation times, right? Or at least comparable to the standard FID measurements. So without the relaxation time, we only get half of the picture. Don't get me wrong, the half we are getting with steady-state is very good but a short hint/comment on the relaxation time issue would be nice (and objectively fair) I guess.
L.90: Maybe it is worthy to explain the term "drooping"? I have no electrical engineering background and I could basically infer from the context what it means. So maybe to an inexperienced reader a short explanation would be quite helpful?
Fig.1: I don't know what the policy of the journal is but I personally like to have the figure explained to me in the text also and not just in the figure caption. I can live with it as it is but maybe this is a point to reconsider?
L.125-L.130: You use the letter "T" for temperature (T_j) and time (T_p). This is bad practice. Consider maybe using also the letter \tau for time related variables (with the obvious exception of T1 and T2) like you do in eq. 2.
L.148: As you are already referring the Fig.8b (also a bit of bad practice), the caption says (c) instead of (b)
Paragraph 5.1: I guess it is clear but you could explicitly mention that Fig. 5 shows only simulations.
L.235: Technically, you could merge Fig.7 into Fig.4 so that the reader already there gets an impression of the decreased dimensions of the customized heatsink.
L.252: Description of Fig.8b: Can you comment on the variability of the pulse amplitude after the initial drooping. So from about 3s onward. Why does the intra-pulse drooping in the Apsu is much more pronounced compared to the simulations for 10mF (Fig.8a)?
Paragraph 5.3: Maybe Fig.10 should come before Fig.9 as your explanation also starts with Fig.10?
L.316: I like the beginning of last sentence ("We foresee") but would rather take it with the wink of an eye.
Technical corrections:

L.29 - Earth's magnetic field (also in L.161)

L.39 - are of small amplitude or have a small amplitude

L.43 - have been very

L.63 - standard one

L.159 - contain instead of containing

L.175 - vary from

L.316 - and has been

Citation: https://doi.org/10.5194/gi-2023-5-RC2
- AC2: 'Reply on RC2', Jakob Juul Larsen, 13 Nov 2024
  
  Dear Reviewer 2
  
  Thanks for the kind review and constructive comments for the manuscript.
  Specific comments:
  Do you really need the TEM part in the introduction?
  
  Yes and no. This paper is certainly not about TEM, but we have more and more field work where the acquisition of both TEM and NMR data allows us to draw much firmer conclusions than if data was collected with either TEM or NMR alone. We have therefore decided to keep the few lines on TEM to address this fact.
  
  L.33: To me, relaxation time is only indirectly controlled by porosity. It is controlled by the pore size and mineral components of the pore walls. The amplitude is more closely connected to porosity because of: NMR amplitude -> (mobile) water content -> porosity
  
  The text will be revised to “…relaxation time of the NMR signals is controlled by the host material, e.g., pore size and mineral composition,…”
  
  L.57: Isn't the 5s argument not a bit stretched? This strongly depends on the T1 of your subsurface and I can rarely remember ever seeing T1 above 500ms in the field. So maybe it is more fair to link the experiment time to the T1 time of the subsurface?
  
  The wait time is limited by the need for re-magnetization, but obviously this time is not known at the time of measurement. Further, the wait time can also be affected in some instruments by the need to recharge capacitors. The sentence will be revised into “…a wait time of typically several seconds…”
  
  L.60-68: I totally agree on your point about the new steady-state approach and the massive SNR enhancement. Yes, this helps us to get better data to estimate water content (and therefore porosity) but still lacks the possibility to infer proper relaxation times, right? Or at least comparable to the standard FID measurements. So without the relaxation time, we only get half of the picture. Don't get me wrong, the half we are getting with steady-state is very good but a short hint/comment on the relaxation time issue would be nice (and objectively fair) I guess.
  
  We would like to stress that steady-state surface NMR is sensitive to amplitude as well as relaxation time. See for example the inversion results in this paper: M. P. Griffiths, D. Grombacher, L. Liu, M. Vang, and J. J. Larsen, Forward modelling steady-state free-precession in surface NMR, IEEE Transactions on Geoscience and Remote Sensing, DOI 10.1109/TGRS.2022.3221624, 2022.
  
  L.90: Maybe it is worthy to explain the term "drooping"? I have no electrical engineering background and I could basically infer from the context what it means. So maybe to an inexperienced reader a short explanation would be quite helpful?
  
  This is a good point. We will modify the text to “…consistent current drooping (the current amplitude decreases during the pulse) within ..”
  
  Fig.1: I don't know what the policy of the journal is but I personally like to have the figure explained to me in the text also and not just in the figure caption. I can live with it as it is but maybe this is a point to reconsider?
  
  We will move the majority of the caption text to section 2, “Challenges in steady-state surface NMR”.
  
  L.125-L.130: You use the letter "T" for temperature (T_j) and time (T_p). This is bad practice. Consider maybe using also the letter \tau for time related variables (with the obvious exception of T1 and T2) like you do in eq. 2.
  
  This was also pointed out by reviewer 1 with an almost identical suggestion. We will update the notation to according to the reviewer 1 suggestion.
  
  L.148: As you are already referring the Fig.8b (also a bit of bad practice), the caption says (c) instead of (b)
  
  This will be fixed in the revision.
  
  Paragraph 5.1: I guess it is clear but you could explicitly mention that Fig. 5 shows only simulations.
  
  This will be explicitly mentioned in the revision. Please also see our response to reviewer 1 on this.
  
  L.235: Technically, you could merge Fig.7 into Fig.4 so that the reader already there gets an impression of the decreased dimensions of the customized heatsink.
  
  This is indeed a possibility, but we would prefer to keep them as separate figure as they have different scope. Figure 4 is a picture of the experimental setup, whereas Figure 7 shows the “before and after” result.
  
  L.252: Description of Fig.8b: Can you comment on the variability of the pulse amplitude after the initial drooping. So from about 3s onward. Why does the intra-pulse drooping in the Apsu is much more pronounced compared to the simulations for 10mF (Fig.8a)?
  
  This is a sampling artefact. In the experiments we use pulses with a 2127 Hz Larmor frequency. The duration of the pulses and the time between pulses are defined by the Larmor period. As this doesn’t match well with the sampling rate used to measure coil current data for this figure, we are essentially sampling the current at different times during each (drooping) pulse. This gives rise to the observed variations of about 1 A after 3 s. The right way to assess the stability of the pulses after the initial drooping is the analysis presented in Figures 9 and 10.
  
  The most obvious reason for the discrepancy between simulations and experiments is that the simulations are based on an ideal DC supply feeding the capacitor bank, but the DC supply used in practice is not ideal.
  
  We will add comments on both these questions during the revision.
  
  Paragraph 5.3: Maybe Fig.10 should come before Fig.9 as your explanation also starts with Fig.10?
  
  This is a simply latex copy-paste error, and it will be fixed in the revision.
  
  L.316: I like the beginning of last sentence ("We foresee") but would rather take it with the wink of an eye.
  
  It is indeed a speculative comment. Time will tell if we are right.
  Technical corrections:
  L.29 - Earth's magnetic field (also in L.161)
  L.39 - are of small amplitude or have a small amplitude
  L.43 - have been very
  L.63 - standard one
  L.159 - contain instead of containing
  L.175 - vary from
  L.316 - and has been
  These will be fixed in the revision process.
  
  Citation: https://doi.org/10.5194/gi-2023-5-AC2

Nikhil B. Gaikwad, Lichao Liu, Matthew P. Griffiths, Denys Grombacher, and Jakob Juul Larsen

Viewed

Total article views: 885 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
695	108	82	885	37	36

HTML: 695
PDF: 108
XML: 82
Total: 885
BibTeX: 37
EndNote: 36

Views and downloads (calculated since 29 Jun 2023)

Month	HTML	PDF	XML	Total
Jun 2023	13	6	3	22
Jul 2023	152	21	7	180
Aug 2023	103	7	4	114
Sep 2023	94	0	94
Oct 2023	78	6	2	86
Nov 2023	29	0	29
Dec 2023	47	2	2	51
Jan 2024	33	3	3	39
Feb 2024	6	4	0	10
Mar 2024	10	5	2	17
Apr 2024	13	2	2	17
May 2024	7	5	6	18
Jun 2024	11	2	2	15
Jul 2024	18	5	3	26
Aug 2024	11	6	1	18
Sep 2024	20	3	41	64
Oct 2024	29	19	1	49
Nov 2024	17	10	3	30
Dec 2024	4	2	0	6

Cumulative views and downloads (calculated since 29 Jun 2023)

Month	HTML	PDF	XML	Total
Jun 2023	13	6	3	22
Jul 2023	152	21	7	180
Aug 2023	103	7	4	114
Sep 2023	94	0	94
Oct 2023	78	6	2	86
Nov 2023	29	0	29
Dec 2023	47	2	2	51
Jan 2024	33	3	3	39
Feb 2024	6	4	0	10
Mar 2024	10	5	2	17
Apr 2024	13	2	2	17
May 2024	7	5	6	18
Jun 2024	11	2	2	15
Jul 2024	18	5	3	26
Aug 2024	11	6	1	18
Sep 2024	20	3	41	64
Oct 2024	29	19	1	49
Nov 2024	17	10	3	30
Dec 2024	4	2	0	6

Viewed (geographical distribution)

Total article views: 858 (including HTML, PDF, and XML) Thereof 858 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 13 Dec 2024

Short summary

The work presents simulations, modelling, and experimental verification of a novel steady-state surface NMR transmitter used for the non-invasive exploration of groundwater. The paper focuses on three main aspects of high current transmitter instrumentation, i.e., thermal management, current drooping, and pulse stability. This work will interest readers in geoscientific instrument prototyping for groundwater exploration using portable geoscientific instrument.


Total:	0
HTML:	0
PDF:	0
XML:	0